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ABSTRACT 

Using  the  pseudopotential  formalism,  calculations 
are  made  for  the  residual  resistivity  of  alloys,  devia¬ 
tions  from  Matthiessen ’ s  rule  and  the  resistivity  due 
to  dislocations  in  metals. 

Residual  Resistivity  of  Alloys:  It  is  shown  that 
the  failure  of  Harrison  (19 66)  to  obtain  a  good  agreement 
with  experimental  residual  resistivities  is  primarily  due 
to  the  fact  that  his  calculation  of  the  form  factors  of 
an  impurity  ion  takes  a  rather  limited  account  of  the 
effect  on  them  of  the  change  in  the  surroundings  of  the 
ion  compared  to  the  case  of  the  ion  in  the  pure  metal. 
Taking  a  proper  account  of  this  effect  and  on  considering 
a  few  further  improvements  (nonlocal  screening,  effective 
mass  etc.)  a  satisfactory  agreement  with  experiment  is 
obtained . 

Deviations  from  Matthiessen ' s ' Rule :  Attention 
is  drawn  to  an  interesting  cause  of  the  deviations  A 
from  Matthiessen ' s  rule  which  arise  from  the  interference 
between  the  scattering  due  to  lattice  vibrations  of  the 
host  atoms  and  the  excess  potential  due  to  the  impurities. 
A  semiquant itative  comparison  with  experiment  based  on 
the  impurity  form  factors  introduced  above,  shows  that 
for  simple  metals  at  high  temperatures,  the  deviations 
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resulting  from  such  an  interference,  make  a  significant 
contribution  to  A. 

Resistivity  due  to  Dislocations:  Treating  the 
dislocation  as  a  row  of  weak  scatterers  with  a  spheri¬ 
cally  symmetric  potential,  an  exact  solution  of  the 
Boltzmann  equation  in  the  presence  of  a  set  of  parallel 
dislocations  is  obtained.  This  is  used  for  making  a 
systematic  calculation  of  the  resistivity  in  the 
Harrison  model  of  the  dislocation.  It  is  also  pointed 
out  that  the  usual  method  of  calculating  the  resistivity 
due  to  random  dislocations ' from  the ~ resist ivity  of 
parallel  dislocations  would  be ' violated' under  certain 
conditions  and  an  alternative  method  is  explored. 
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INTRODUCTION 


This  thesis  is  concerned  with  a  study,  based  on  the 
pseudopotential  formalism,  of  (a)  the  residual  resistivity 
of  alloys  and  the  deviations  from  Matthiessen ' s  rule,  and 
(b)  the  resistivity  due  to  dislocations. 

The  first  attempts  to  explain  quantitatively  the 
electrical  resistivity  phenomenon  came  after  the  advent 
of  quantum  mechanics,  particularly  through  the  work  of 
Bloch,  Nordheim,  Mott  and  Bardeen. 

The  essential  step  in  the  estimation  of  the  electri- 
cal  resistivity  is  the  calculation  of  P(k,k’),  the 
probability,  per  unit  time,  that  an  electron  in  a  state 
characterized  by  the  wave  vector  k  is  scattered  to  another 
state  of  wave  vector  k’ .  This  scattering  arises  from  the 
aperiodicity  of  the  lattice  potential  caused  by  disorders 
relating  to  the  structure  (thermal  vibrations,  lattice 
distortions  etc.)  and/or  to  the  composition  (presence  of 
impurities).  The  f orement ioned  works  are  based  on  obtaining 
P(S,lc 1  )  by  calculating  the  perturbing  potential,  the 
difference  in  the  potential  in  the  perturbed  and  the  perfect 
lattice,  due  to  the  various  mechanisms.  This  aspect  of  the 
subject  has  been  reviewed  in  many  excellent  texts  [Mott  and 
Jones  (1936),  Wilson  (1953),  Ziman  (I960)]. 

i 

In  the  last  decade  an  alternative  approach,  known  as 
the  effective  potential  theory,  has  been  developed  by 
Phillips  and  Kleinman  (1959),  Austin,  Heine  and  Sham  (1962), 
Harrison  (1963),  Ziman  (1964),  and  Abarenkov  and  Heine(1965). 
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The  advantage  of  this  approach  is  that  since  the  effective 
potential  is  weak,  one  can  work  directly  with  the  total 
potential,  instead  of  with  the  perturbing  potential.  The 
scattering  can  then  be  written  as  a  product  of  atom  form 
factors  and  the  structure  factor  (or  the  partial  structure 
factors  in  the  case  of  an  alloy).  The  former  depend  only 
on  the  individual  (pseudo-)  atoms  and  the  latter  only  on 
their  arrangement  in  the  crystal.  Such  a  factorization  is 
a  consequence  of  the  weak  scattering  formalism  and  has  been 
known  for  a  long  time  in  X-ray  and  neutron  scattering 
theories.  Its  application  to  electron  scattering  and 
electrical  resistivity  problems  was  first  made byKrishnan 
and  Bhatia  (19^5)  and  Bhatia  and  Krishnan  (19^8).  The 
effective  potential  formalism  enables  one  to  calculate  the 
atom  form  factors  and  hence  the  electrical  resistivities 
quantitatively  [Ziman  (1961),  Faber  and  Ziman  (1965)  and 
Harrison  (19633  1966  )3- 

The  effective  potential  theory  has  been  formulated  in 
three  different  ways.  The  earliest  approach  is  known  as 
the  pseudopotential  method  which  is  the  one  used  in  this 
thesis.  The  other  two  are  the  quasi-  and  the  model  potentials. 
The  three  methods  differ  from  each  other  only  in  the  calcula- 

I 

tion  of  the  form  factors.  The  calculation  of  the  properties 
of  metals  and  alloys  proceeds  exactly  the  same  way  in  all  the 
three  methods. 

The  thesis  is  divided  into  three  parts.  Part  I  deals 
with  the  background  material.  Part  II  discusses  the  residual 
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resistivity  of  dilute  alloys  and  deviations  from 
Matthiessen ’ s  rule.  Part  III  is  devoted  to  the  resistivity 
due  to  dislocations. 

In  Chapter  1  we  collect  together  some  formulae  for 
calculating  the  electrical  resistivity  if  P(k,k’)  is  known. 

In  Chapter  2  we  review  the  pseudopotential  technique  which 
forms  the  basis  of  our  calculations  of  P(k,kf). 

In  Chapter  3  we  show  that  the  environment  of  an,  ion, 
which  changes  on  alloying,  has  a  much  greater  effect  on  the 
form  factors  of  the  ion  than  had  been  considered  in  the 
calculations  of  Harrison  (1966).  On  calculating  the  impurity 
form  factors  by  taking  a  proper  account  of  this  fact,  the 
agreement  with  experimental  residual  resistivities  is  found 
to  be  much  better  than  obtained  by  Harrison. 

In  Chapter  4  we  point  out  an  interesting  cause  of  the 
deviation  (A)  from  Matthiessen ' s  rule  which  arises  from  the 
interference  between  the  scattering  due  to  lattice  vibrations 
of  the  host  atoms  and  the  excess  potential  due  to  the 
impurities.  Such  interference  terms  are  missed  out  in  most 
previous  treatments  of  the  problem  [Koshino  (I960),  Taylor 
(1962,  1964)  and  by  Damon,  Mathur  and  Klemens  ( 1 9 6 8 ) J  (and 
references  given  there)  because  of  the  usual  practice  of 
adding  the  probabilities  of  scattering,  rather  than  the 
scattering  amplitudes,  due  to  the  two  mechanisms.  We  shall 
see  that  these  interference  terms  play  an  important  part  in 
determining  both  the  qualitative  and  quantative  behaviour 
of  A,  particularly  at  high  temperatures. 


' 
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In  Chapter  5  we  give  a  brief  survey  of'  the  experimental 
situation  and  the  previous  theories  of  the  resistivity  due 
to  dislocations.  Chapter  6  and  7  give  the  results  of  our 
investigation  of  this  problem. 


PART  I 


THE  BACKGROUND  MATERIAL 
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Chapter  1 

METHODS  OP  SOLVING  THE  BOLTZMANN  EQUATION 
AND  THE  CALCULATION  OP  ELECTRICAL  RESISTIVITY 

§1.1  Elementary  solution  of  the  Boltzmann  Equation 

The  standard  treatment  of  electrical  resistivity 
entails  the  solution  of  the  Boltzmann  transport  equation, 
which  in  the  steady  state,  is  given  by 


3  t 


J  field 


3fr 


9  t  , ,  .  . 

*-  J  collision 


=  0 


(1.1) 


where  fj  is  the  distribution  function  of  the  electrons  i.e. 
1  -> 

- 5-  ffdk  is  the  number  of  electrons  of  a  given  spin  per 

(2tt)  J  k 

unit  volume  of  wave  vector  space  dk. 

In  (1.1)  the  first  term  arises  from  the  drift  of  the 
electrons  in  the  presence  of  the  electric  field  e,  and  is 
given  by 


field 


-*  3  ft 

ee  _ k 

n  '  35 


(1.2a) 


For  weak  fields  this  may  be  approximated  to 


r9 

_ k 

9  t 


field 


-* 

ee 

fi 


9  f 


o 

1 


9  k 


(1.2b) 


where  f2  is 
k 

function . 


the 


equilibrium 


(Fermi-Dirac ) 


distribution 


. 
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The  second  term  in  (1.1)  represents" the' change  in 
f  due  to  collisions.  If  P(K,Kf)  denotes  the  probability. 


per  unit  time,  that  an  electron  in  state  k  makes  a  tran¬ 


sition  to  state  k’ ,  then 


r3f 


k 


3  t 


J  coll 


is  given  by 


r3f£ 


3  t 


J  coJL.1 


=  -/{fj(i-f£t  )p(k,5»  )-f+,  (l-f^)p(J'  35)}d£.» .  (1.3) 


k 


k 


k 


In  general,  the  Boltzmann  equation  cannot  be  solved 
exactly.  However  if  (1)  the  scattering  is  elastic  i.e.. 
Eg  =  Eg,  and  hence  P(k,k’)  =  P(k’,S),  (2)  P(k,k' )  depends 
only  on  j  k— k '  |  ,  and  (3)  Fermi  surface  is  spherical,  then 

a  simple  solution  for  f  can  be  readily  obtained.  To  see 
this  we  first  note  that  for  weak  fields,  f  can  be  taken, 
quite  generally,  to  be  of  the  form 


Following  the  normal  practice  we  shall  take  the  crystal 

volume  ft=l  and  omit  the  factors  and  —  in  the  kf-  and 

8u^  4  TT 

k-integrat ions .  Note  that  in  the  k-integrat ion  we  have  a 
factor  2  arising  from  the  spin  multiplicity;  it  does  not 
appear  in  the  k ’ -integration  because  scattering  takes  place 
only  between  some  spin  states. 
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o 

fi*-  =  f^  —  ( 
k  k 


ee 

fi 


o 

3  f  g 

grad5  E)  3Ej  * 


(1.4) 


where  xg  has  the  dimensions  of  time  and  is,  as  yet,  an 
unknown  function  of  5.  Now  for  elastic  scattering  (1.3) 
reduces  to 


3  t 


J  coll 


-  /  (fj-  fj, )p(5,S’ )dKf 


(1.5) 


and  hence  substituting  (1.4)  into  (1.5)  and  using  (1.1) 
and  (1.2b),  one  obtains  for  x+  the  integral  equation 


->  ■+ .  .  ,  -> 


e -grad^.E 

1  =  7(xt>-  x£,  - - )P(K,K’)dk’  . 

e.gradj;  E 


(1.6) 


For  a  spherical  Fermi  surface  (1.6)  reduces  to 


1  = 


/  (TjJ- 


T£- 


sdL)p(J  £.)d£. 

£  .  k 


(1.7) 


If  one  now  takes  PGc,S')  =  P(|K-ic',|)  so  that  P(ic,E' ) 


->  -> 


depends  only  on  the  angle  0  between  the  incident  and 


-*■ 


-> 


the  scattered  directions  (and  not  separately  on  k  and  kT), 

* 

then  one  obtains  for  xj  the  expression 


1/x^  =  /  (1  -  cos0 ) P ( | k-5 T | ) d5 ’  , 


(1.8) 


%  y  . 

Note  that  P(k,kf  )  is  proportional  to  6(E+-  E+f  )  so  that 


the  integral  in  (1.7)  is  essentially  a  two  dimensional  one 
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which  is  independent  of  the  direction  of  Ic. 

Using  (1.8)  the  resistivity  may  be  readily  obtained 
as  below. 

The  current  density  J  is  given  by 


J  =  e  /  Vf  ff  d£ 
k  k 


(1.9) 


Substituting  (1.4)  and  restoring  the  factor 

had  been  left  understood,  one  obtains 

o 


J 


2 

e 

4  7T  3h 


3  f 


k 


3  E 


k 


dE,  dS. 
k  k 


1/4tt3 


which 


A  — 

where  k  is  a  unit  vector  in  the  direction  of  k  and  dS, 

k 

is  a  surface  element  on  the  energy  surface  in  k-space  with 

energy  E^.  Taking  the  components  of  J  along  arbitrarily 

chosen  Cartesian  axes,  one  gets 

o 


J. 

i 


47T3fl2 


dE,  3  f ,  /s  ^ 

f  Tk  d\T  3“E”  dEk-/'(^e-k)k1dSk 
k 


2 

e  n 


o 


3  TT  ft 


0  dE,  3  f, 

/  t  k2  — —  — —  he 
2^2  J  L  k  dk  3  E,  k 


k 


(1.10) 


Remembering  that  for  any  smoothly  varying  function  g(E) 


f  g(E) 


3  f 


o 


3  E 


dE 


•g(Ep) 


for  knT  <<  E„,  one  gets 
Jd  r 


J. 

1 


2 

e  e . 


i 


?  2*2 
3  tt  n 


dE,-, 

_ k 

dk 
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(1.11) 


(1.12) 
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Comparing  this  with  the  definition  of  the  conductivity 


tensor 


J  .  —  Z  G  •  •  £  . 

1  j  1J  J 


one  finds  that  for  this  case 


a .  . 

ij 


2  2 
3tt  fi 


o  dE,  ■ 

i  2  k 

t,  k  -t; — 
k  dk 


F 


(1.13) 


1  =  J  5 


(1.14) 


=  0 


(1.15) 


Thus  here  the  conductivity  tensor  is  represented  by  a 
single  number  g^  which  is  the  same  for  i=l,2,3  and  we 
may  write  it  as  g,  the  conductivity.  Its  reciprocal  p  is 
the  resistivity. 

Notice  that  a  in  (1.14)  is  proportional  to  the  square 

For 


dEk  dE^ 

of  (jjk- )p  -  the  other  factor  of  arises  from 


metals  where  the  Fermi  surface  can  be  assumed  to  be 
spherical,  it  is  usually  possible  also  to  assume  that  the 
E-k  relation  is  parabolic  so  that  for  these  metals  o  is 
proportional  to  (m*/m)~’  |^or  p  «  (m*/m)  J,  where  m*  is  the 
effective  mass  for  the  electrons  in  the  metal  and  m  the 
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free-electron  mass.  For  this  case  (1.14)  can  be 
written  in  the  convenient  form 


o  = 


2 


ne 


m* 


(1.16) 


where  n  is  the  density  of  electrons.  Although  n  appears 
explicitly  in  (1.16),  it  may  be  noted  that,  as  will  be 
clear  from  (1.14),  the  conductivity  (  and  other  transport 
properties)  are  determined  entirely  by  the  electrons  at 
the  Fermi  surface.  The  low-lying  electrons  are,  of  course, 
important  but  only  in  building  up  the  Fermi  level. 


x  2 

It  may  be  mentioned  here  that  the  appearance  of  (m*/m) 

in  the  expression  for  a  has  been  a  subject  of  some  doubt 
since  the  work  of  Edwards  (1962)  and  Faber  (1966)  who  gave 
a  formalism  for  calculating  a  without  using  the  Boltzmann 
equation  (according  to  these  authors  a  should  be  indepen¬ 
dent  of  (m#/m));  For  a  recent  review  on  this  point  see 
Mott  (1967);  see  also  §3-3  of  the  thesis. 


3  v  ■  i  i  ■  i  -  ' : 

Q-ssuqs  n  i  ■  rod  r IA  .Sf.o^osle  *  o  s  isb  sd :  o.  n 
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§L .  2  Discussion  of  the  General  Case  and 
Variational  Principle 


In  the  general  case  when  P(k,k* ) depends  on  the 
directions  of  k  and  k'  separately  or  (and)  the  Fermi 
surface  is  not  spherical,  xg  would  depend  on  the  direction 
of  S  and  the  integral  equation  (1.6)  or  (1.7)  has  to  be 
solved  for  each  specific  case.  An  alternative  but  equi¬ 
valent  method  of  solving  the  Boltzmann  equation  for  these 
cases  consists  of  in  expanding  distribution  function  in 
spherical  harmonics,  which  we  shall  use  in  Chapter  6  on 
dislocations.  In  general,  however,  one  has  to  resort  to 
approximate  methods,  in  particular,  the  variational 
technique  which  is  applicable  for  both  elastic  and  in¬ 
elastic  scattering. 

Before  describing  this  technique,  we  remark  that 

the  conductivity  tensor  a. .  will  now,  in  general,  not  be 

— ■  J 

a  scalar.  However,  it  follows  from  Onsager’s  reciprocity 

relations  that  a..  =  a...  Now  for  a  symmetrical  tensor 

ij  j  l 

of  rank  2,  it  is  always  possible  to  choose  a  set  of  axes, 
called  the  principal  axes,  with  respect  to  which  the 
non-diagonal  components  of  the  tensor  are  zero.  With  respect 
to  these  axes,  the  non-diagonal  components  of  the  resisti¬ 
vity  tensor  p . . 

p..  =  j  (cofactor  of  a..) 

1 J  A  -L  J 


arfcf  no  aflneqebrM)*  nsriw  sbb'd  lens  nos  n  ^ 

Ians’*  erf*  (bnfi)  no  ^Isctsnfiqse  f  bn&  o  anc.:.  •  -  r.  b| 

,d  o.  a£l1'  ( V  -  X  >  no  (d.J)  noitfsups  Isng  octal  bn-J  ^  ‘  ° 

■ 

c  ri  r:w  t  .?x  '  '  "  1 
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& 

where  A  Is  the  determinant  of  \o..\,  will  also  be  zero. 

Clearly  along  the  principal  axes,  J||e  and  for  such 
cases  the  variational  expression  for  the  resistivity  is 
given  by  [see,  for  example,  Ziman  (i960)] 


P  = 


2lTT  ;  (<h-  <(>k')2  d£' 

B  ,  >4^ 

1  evs  n  3-e7  dk  1 

k 


where  Q(k,k')  is  defined  by 


(1.17) 


Q(k,5’  )  [=Q(£'  ,£)]=  (l-f|,  )P(k,5'  )  , 


(1.18) 


and  <()g  is  a  trial  function 

bution  function  fg  with  the 
o 

tion  fj  by  the  relation 

o 


o 


3  f 


*5 


k 


3  Eic 


connecting  the  perturbed  distri- 
equilibrium  distribution  func- 


(1.19) 


The  variational  principle  states  that  <j>g  must  be  chosen  to 
make  p  of  (1.17)  a  minimum. 

Thus  the  problem  of  finding  p  boils  down  to  (i) 
selecting  a  trial  function  with  adjustable  parameter/s. 


For  cubic  symmetry,  the  three  principal  (diagonal)  compo¬ 
nents  of  a. .  or  p. .  are  equal  and  a  and  p  are  scalar, 
ij  iJ 


13 


(11)  determining  the  parameters  by  minimizing  (1.17),  and 
(iii)  substituting  the,  now  known,  parameters  back  into 
(1.17)  to  get  p.  Reasonable  results  can  often  be  obtained 
by  using  the  simplest  trial  function 

4>j*  =  a  £  .  e  ,  (1.20) 

where  e  is  a  unit  vector  in  the  direction  of  the  electric 
field  and  a  is  a  constant  [the  value  of  which  does  not 
matter  since  it  cancels  out  from  (1.17)]. 


,19  3*10  s  .  .Q  i  •  .  1 
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Chapter  2 

PSEUDOPOTENTIALS  IN  METALS  AND 
THE  CONCEPT  OF  WEAKLY  SCATTERING  NEUTRAL  PSEUDO-ATOMS 

One  of  the  most  interesting  things  about  the  valence 
electrons  in  metals  is  the  extent  to  which  their  observed 
properties  can  be  explained  in  terms  of  models  where 
either  the  ions  are  ignored  or  their  role  vastly  simplified, 
except  for  maintaining  charge  neutrality.  Another  interest¬ 
ing  thing  is  that  it  is  not  readily  apparent  why  such  models 
should  work  at  all:  The  nearly  free  electron  approximation, 
for  example,  requires  the  atomic  potential  to  be  weak 
compared  with  the  bandwidth,  and  this  certainly  is  not  the 
case.  It  was  pointed  out  by  Phillips  and  Kleinmann  (1959) 
(to  be  abbreviated  PK)  that  the  reason  for  this  is  that  in 
the  region  of  the  ion  core  where  the  potential  is  the 
deepest,  there  is  an  almost  complete  cancellation  between 
the  large  negative  potential  energy  and  the  large  positive 
kinetic  energy  associated  with  the  rapid  oscillations  of 
the  wave  function  at  the  core.  The  observable  effects  are 
produced  by  a  weak  net  potential  which  has  come  to  be 
called  the  pseudopotential .  Not  only  did  this  provide  an 
immediate  understanding  of  the  nearly-free  behaviour  of 
the  electrons,  it  set  the  stage  for  the  development  of  the 
theory  of  metals  from  an  entirely  new  point  of  view. 

Previously,  whenever  one  wanted  to  calculate  the 
field  due  to  some  arrangement  of  ions  -  for  example,  in 


•  B  flJ3  5  1  b"tor,  ' 
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the  neighbourhood  of  a  dislocation  In  an  otherwise  perfect 
lattice  -  one  was  confronted  with  the  problem  of  allowing 
for  the  redistribution  of  charge  in  the  electron  gas  and 
for  the  field  that  this  produces.  The  construction  of  a 
self-consistent  field  in  such  circumstances  was  a  major 
difficulty.  However,  if  the  effective  interaction 
between  the  ions  and  the  electrons  is  weak  we  can  use 
the  approximation  of  linear  screening  according  to  which 
the  perturbations  of  the  electron  gas  due  to  different 
ions  can  be  added  independently  provided  the  screening 
cloud  is  not  too  large  a  fraction  of  the  total  electron 
density.  This  now  leads  to  the  result  that  in  any 
arrangement  of  the  ions,  the  ion  and  its  screening  charge 
may  be  considered  as  a  single  entity.  The  entity  is 
electrically  neutral  and  is  called  a  pseudo-atom.  After 
we  have  built  up,  say,  a  dislocation  out  of  these  neutral 
pseudo-atoms,  there  can  be  no  further  charge  shift;  all 
the  conduction  electrons  have  been  automatically  accounted 
for  and  their  potential  is  already  contained  in  the 
potential  of  the  pseudo-atoms.  Furthermore,  the  field 
due  to  each  pseudo-atom  being  the  (linearly)  screened 
pseudopotential  of  the  real  ion,  the  scattering  due  to 
each  pseudo-atom  is  weak  and  in  any  problem  the  information 
relating  to  the  arrangement  of  the  ions  enters  only  through 
a  structure  factor  (or  structure  factors  if  we  are  dealing 
with  an  alloy).  Thus,  we  can  hope  to  apply  the  theory, 
as  has  indeed  been  done,  to  a  study  of  the  general  proper- 
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ties  rather  than  simply  the  properties  of  perfect, 
undistorted  crystals. 

It  should  be  emphasized  that  the  concept  of  the 

weakly  scattering  pseudo-atom, despite  its  simplicity 

and  apparent  naivety,  is  not  a  phenomenological  model. 

It  is  based  upon  the  true  physical  system  and  is  almost 

as  precise  in  detail  as  the  best  calculations  that  have 

been  made.  However,  it  may  be  mentioned  that  the 

pseudopotential  formalism  is  based  on  the  small-core 

approximation,  i.e.  overlap  of  core  wave  functions  of 

* 

adjacent  ions  is  negligible  .  Because  of  this  the  method 
can  be  applied  only  to  the  alkali  and  poly  valent  metals, 
or  simple  metals  as  they  are  called.  Its  application  to 
noble  metals  is  questionable  and  to  transition  metals, 
it  is  quite  inappropriate. 

In  this  chapter  we  wish  to  describe  the  pseudo¬ 
potential  formalism,  the  concept  of  the  pseudo-atom  and 
the  calculation  of  the  scattering  from  a  single  pseudo¬ 
atom.  Since  the  first  two  points  relate  to  the  general 
properties  of  pseudopotentials,  the  discussion  of  these 
points  would  apply  to  pure  metals  as  well  as  alloys. 


This  remark  equally  well  to  the  other  two  formulations 


of  the  effective  potential  theory. 
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However,  the  explicit  calculation  of  the  scattering  from 
a  pseudo-atom  will  here  be  discussed  only  for  a  pure 
metal,  the  modification  of  the  scattering  on  alloying 
to  be  taken  up  later. 

§2.1  The  Pseudopotential  Method  -  Formal  Work 

§2.1.1  The  P-K  Formulation 

The  formulation  of  Phillips  and  Kleinmann  employs 
the  orthogonalized  plane  wave  (OPW)  method  which  is  based 
on  expanding  the  (true)  wave  function  ijjg  in  terms  of 
plane  waves  orthogonalized  to  the  core  states.  Thus 

f'K  =  l  aq  (*b  Xg+5  >  (2-D 

q 

where  the  orthogonalized  plane  waves  x£  are  given  by 

XJ  =  | 5  >  -  Za> <  a|k  > 

a 

E  (1  -  P) |£  >  ,  (2.2) 

the  projection  operator  P  being  given  by 

P  =  £ | a><a |  .  (2.3) 

a 

Here  |  j£>  are  plane  waves  normalized  to  unity  over  the 
crystal  volume  ft,  |a>  are  the  core  states,  the  term 
P|Io  in  x  j  is  for  making  the  latter  orthogonal  to  |ot>. 


■;  <  '  '  7 


■  t. 
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Now  ipj*  satisfies  the  equation 

(T  +  V)  ijjg  =  Eg  ipg  (2.4) 

2  2 

where  T  is  the  kinetic  energy,  -h  V  /2m,  V  is  the  true 
potential  and  Eg  the  eigenvalue  of  the  Hamiltonian  for 
the  eigenfunction  .  On  substituting  (2.1)  and  (2.2) 
in  (2.4),  using  (T+V)|a>  =  E  |a>  where  E  is  the  energy 
of  the  core  state  |a>  ,  and  rearranging  we  get 

(T  +  W)  (f)g  =  Eg  <j>g  (2.5) 

where 

(f)g  =  Zag(S:)|5  +  q>  (2.6) 

q 

and 


W6k  =  V(j)g  +  E(E^  -  E^  )  |  ot><ot  | 

a 


k  '  "  v  ^k  y  i  i  Yk' 


(2.7) 


Thus  the  true  wave  function  ij/g  and  the  true  potential  V 
may  be  replaced  by  <f>g  and  W  respectively  without  affecting 
the  energy  eigenvalue.  <j>g  is  called  the  pseudo  wave 
function  and  W  the  pseudopotential. 

We  notice  the  following  points 
(a)  From  (2.1),  (2.2)  and  (2.6)  we  get 


i4g  =  <J>£  -  Z  |  a><a  j  (p£> 


k 


a 


k' 


(2.8) 


so  that  <j>g  and  ^g  are  the  same  outside  the  region  of  the 

Using  this  and  (2.5)  it  can  be  proved  (Austin,  Heine 


core . 


« 

' 

.  ;  I '  *  l ;  I 


19 


and  Sham  (1962))  that  the  equivalence  of  (ijjj,  V)  problem 
to  the  (<J>£,  W)  problem  is  carried  over  to  the  determina¬ 
tion  of  the  scattering  probability. 

(b)  The  second  term  in  (2.7)  is  positive.  For  this 
reason  it  is  called  the  repulsive  potential  and  has  the 
effect  of  largely  cancelling  V  in  the  core  region  so 
that  we  can  use  perturbation  theory  to  calculate  the 
scattering  and  the  energy  of  the  eigenstates.  Experience 
shows  that  for  scattering  which  is  what  concerns  us  in 
this  thesis,  the  results  obtained  by  the  lowest  order 
perturbation  theory  are  reasonably  accurate. 

(c)  The  repulsive  part  V\F  of  W  is  an  integral  operator 
or  a  nonlocal  potential: 

WR  cf)(r)  =  /WR(r  ,r  '  )<J>  (r  '  )dr  ? 

WR(r,r’)  =  Z  (E^  -  E  )\p  (r  ’  )\p  (r  ) 

5  a  k  a  ra  ra 

where  ib  is  the  wave  function  for  the  core  state  |a>. 
a  1 

Because  of  this  nonlocal  character  of  W  ,  its  matrix- 
elements  between  two  states  |k>  and|k+q>  will  depend  on 
both  k  and  q,  rather  than  on  q  alone.  Thus  we  have  lost 
some  of  the  simplicity  of  the  original  problem  where 
V(r)  is  a  local  potential.  However,  the  nonlocal  nature 

p 

of  W  only  causes  some  calculat ional  complexity  and  no 
fundamental  difficulty. 


33  ■ 
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§  2 . 1.2.  The  Nonuniqueness  of  PK  pseudo  wave  function  and 

the  AHS  Theorem 

t 

It  was  pointed  out  by  Cohen  and  Heine  (196l)  that 
the  PK  pseudopotential  does  not  generate  a  unique  pseudo 
wave  function.  <j)g;  is  arbitrary  to  the  extent  of  the 
addition  of  a  linear  combination  of  the  core  states. 

If  one  adds  Z  a  |a>to  (J)f,  where  a  is  arbitary,  it  will 
remain  a  solution  of  (2.5)  with  the  same  eigenvalue. 
Further,  one  notes  from  (2.5)  that  the  true  wave  function 
also  remains  unchanged  as  the  addition  to  the  first  term 
is  neutralized  by  a  subtraction  in  the  orthogonalizat ion 
term . 

The  reason  for  this  nonuniqueness  lies  in  the  OPW’s 
Xg:  being  overcomplete,  hence  nonorthogonal .  The  plane 
waves  | £>  are  complete  and  independent,  but  the  effect 
of  orthogonalizat ion  is  to  Introduce  linear  relations 
among  the  Xjg* 

The  arbitrariness  in  4>g;,  allows  us  to  impose  a 

p 

constraint  on  it.  This  will  affect  W  and  hence  W. 

Thus  depending  upon  the  type  of  the  constraint  applied, 
we  will  get  various  forms  for  W. 

It  was  shown  by  Austin,  Heine  and  Sham  (1962)  that 
the  arbitrariness  of  W  is  mpre  deep-rooted.  In  general, 
the  pseudopotential  may  not  be  connected  with  the  OPW’s 
at  all.  They  proved  that  the  most  general  pseudopotential 


is  given  by 


v  *  z  . 
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W(}>g  =  V(r)<f>g(r)  +  Z  ipa  (r )  (f  (r’,a)  ,<j>g(rO  )  ,  (2.9) 

a 

where  f(r’,a)  Is  an  arbitrary  function  of  position  and 
core  index  a.  This  result  is  known  as  the  AHS  theorem. 

In  writing  (2.9)  we  have  used  the  inner-product  notation 

( u (r ),v(r))=/u*(r)v(r)dr.  (2.10) 

The  relation  between  the  true  wave  function  and  the 
pseudo  wave  function  in  the  AHS  formalism  may  be  deduced 
from  (2.9)  and  (2.5)  making  use  of  the  orthogonality  of 
\p+  and  \p  .  It  is  given  by 

(J)V:  =  ipT*  +  £  b  ip  ,  (2.11) 

Yk  rk  a  a, 

a 

where  b  is  determined  from  the  set  of  equations 
a 

b  (5)  (E+  -  E  )  =  Z  f(n,a)a  (k)<  a|n>5  (2.12) 

oc  oc  ^  ri 

where  |n>'s  are  a  complete  set  of  orthonormal  states 
(which  may  be  chosen  as  plane  waves  if  necessary)  and 
a  (k)  and  f(n,a)  are  defined  by  the  expansions 

<J>£  =  £  a  (£)  |  n>  (2.13) 

n 

f(r,a)  =  2  f*(n>a) <  n|a>|n>  •  (2.14) 

n 

Before  we  can  obtain  b  (k)  from  (2.12)  we  will  naturally 

a 

have  to  eliminate  a  (k)  using  (2.11)  and  (2.13). 


!<  i  ■  tn)* 
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The  detailed  evaluation  of  b  is  of  no  interest 

a 

to  us  but  we  may  note  that  when  f  (n,a)  is  independent 

of  n  and  f(a)  =  (Er>:  -  E  ) ,  as  is  the  case  for  the  PK 

k  a 

pseudopotential,  b^  =  <a|<J>£>,  which  is  what  we  obtained 
earlier  in  (2.8).  It  may  be  noted  that  this  does  not 
imply  any  particular  value  of  b  .  It  really  expresses 
only  the  orthogonality  of  and  and  leaves  b^ 
arbitrary,  as  discussed  in  the  beginning  of  §2.1.2. 

§2.1.3  The  Selection  of  a  Pseudopotential 

The  nonuniqueness  of  the  pseudopotential,  discussed 
in  §2.1.2,  not  only  does  not  give  rise  to  any  inconsis¬ 
tency  in  principle,  in  practice  also  it  is  not  entirely 
a  limitation.  On  the  contrary  it  can  be  very  helpful 
because  it  allows  us  to  choose  our  pseudopotential  and 
by  making  a  judicious  choice,  one  may  hope  to  be  able  to 
use  simple  computational  techniques  for  estimating  the 
observed  properties. 

Cohen  and  Heine  ( 1 9 6 1 )  suggested  that  we  choose  a 
pseudopotential  which  leads  to  a  pseudo  wave  function 
with  the  most  rapidly  convergent  plane  wave  expansion. 
Such  a  pseudo  wave  function  might  be  replaced  without 
much  error  by  a  single  plane  wave  -  beyond  doubt  the  most 
convenient  representation  of  (p .  Cohen  and  Heine  further 
noted  that  if  we  require  our  pseudopotential  to  give  the 
smoothest  possible  pseudo  wave  function,  we  might  just 
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get  the  <f>  we  are  looking  for:  it  is  reasonable  to  expect 
that  the  smoothest  <f>  will  have  an  optimum  convergence  of 
its  expansion  in  plane  waves. 

Noting  that  the  criterion  of  smoothness  is  equi¬ 
valent  to  minimizing 


/  |  V<|>  |  ^ dr//  |  (j)  |  ^dr 


(2.15) 


Cohen  and  Heine  found  (see  also  Harrison  (1966)  §8.1)  that 
the  corresponding  W  is  given  by 


W<f>+  =  (1 


(2.16) 


(2.16)  is  not  a  convenient  expression  for  W  because  if  it 
is  solved  explicitly  for  W,  one  gets  a  nonlinear  W;  hence 

(2.16)  is  not  of  AHS  form.  Cohen  and  Heine  suggested 
that  as  an  approximation,  the  second  term  on  the  R.H.S. 
of  (2.16)  (which  gives  rise  to  the  nonlinearity)  be 
dropped.  It  was  noticed  by  Austin  (see  AHS)  that  this  is 
not  an  approximation  at  all,  but  is  still  a  vlid  form 
for  W  satisfying  (2.5).  This  pseudopotential  is  called 
the  Austin  pseudopotential. 

Harrison  proceeding  somewhat  differently  about  the 
second  term  in  (2.l6)  obtained  the  so-called  Harrison 
pseudopotential.  He  replaced  it  by  a  linear  term  such 


that  the  resulting  pseudopotential  is  equivalent  to  the 
pseudopotential  in  (2.16)  up  to  the  second  order  pertur- 


bation  theory.  This  leads  to  a  pseudopotential  given 

by 

W  =  V  +  E  I  a><a  I  ( -fi2V2/2m  +  <S|W|5>  -  E  )  . 
a  a 

Operating  it  on  the  zero  order  state  |  k> 

W|k>  =  V  |  £>  +  Z  I  a>  (fi2k2/2m  +  <k|W|k>  -  E  )  <a  I  £>  .  (2.17) 

i  i  i  ii  a  1 

a 

On  eliminating  <S|W|5>  from  the  right  hand  side  we 
obtain  for  the  matrix  element  of  W  between  two  plane 
waves 


•<k  +  q  I  W  I  k> 

=  <  S+q|v|£>  +  E(fi2k2/2m  +<k|v|k  -E  )  <k+q  |  a><a  |  £> 

a 

+  Z(-h2k2/2m+ <k|  V|k>-E  )<k|a><a|k>  <k+-9.Lp  I  k>  .  (2.18) 

a  a  1-<K|P|S> 

The  smoothness  criterion  given  by  Cohen  and  Heine 
has  been  questioned  by  Abarenkov  and  Heine  (1965)  and 
by  Pendry  (1968).  The  issue  raised  by  Abarenkov  and 
Heine  is  whether  we  should  prefer  the  smoothest  wave 
function  or  the  smoothest  pseudopotential.  They 
suggested  the  use  of  a  criterion  which  is  a  compromise 
between  the  two  but  they  found  that  the  calculations 
became  very  heavy.  Pendry  ( 1968 )  has  suggested  that 
a  more  useful  pseudopotential  may  be  obtained  by 
requiring  the  pseudopotential  to  lead  to  the  most  rapid 
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convergence  of  the  Born  series.  The  usefulness  of  this 
procedure  has  still  to  be  demonstrated. 


§2.2  Pseudopotential  and  the  Use  of 
Perturbation  Theory 

In  §2.1  we  found  that  the  general  pseudopotential 


W<f)+=  V(r)<f)jj(r)  +  £a^a(r)  (f  (r,a)  ,<f>£(r) )  , 


with  <j>g  normalized  according  to 


<f)U  -  +  2  (ip  ip 

Yk  rk  a  a  k  ra 


reproduces  the  true  valence  energies  and  the  true  scatter¬ 
ing.  Now,  if  we  were  to  attempt  an  exact  solution,  we 
may  not  find  the  ( cp ,  W)-problem  to  be  much  simpler  than 
the  original  problem  of  ip  and  V.  However,  the  value  of 
the  pseudopotential  transformation  lies  in  that  in  the 
transformed  problem,  with  a  suitable  choice  of  W,  we  can 
expect  to  obtain  a  reasonably  good  estimate  of  the  exact 
results  for  most  problems  by  the  use  of  the  perturbation 
theory.  It  is  this  feature  which  enables  one  to  explore 
a  vast  array  of  the  properties  of  metals,  rather  than  just 
the  energy-band  calculations. 

The  perturbation  theory  is  essentially  used  in  two 
places:  (l),to  obtain  the  screening  of  the  bare  pseudo¬ 

potential  and  (2),  to  calculate  the  scattering  from  the 
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screened  pseudopotential.  In  both  of  these  we  shall 
restrict  ourselves  to  the  lowest  order  of  perturbation. 

§2.3  Separation  of  the  Pseudopotential 
and  the  Concept  of  Pseudo-atoms 

So  far  we  have  been  talking  of  the  total  pseudo¬ 
potential,  the  pseudopotential  due  to  all  the  ions  and 
the  conduction  electrons.  But  there  is  an  extremely 
important  feature,  known  as  the  separation  of  the  pseudo¬ 
potential  without  which  the  formalism  would  have  been  of 
little  practical  use.  This  consists  in  expressing  the 
total  pseudopotential  as  a  sum  of  the  potentials  centred 
around  individual  ionic  sites.  The  potential  associated 
with  each  site  includes  the  effect  of  the  ion  and  the 
screening  charge  and  is  ascribed  to  a  neutral  entity  called 
a  pseudo-atom.  The  separation  of  the  pseudopotential  is 
achieved  in  two  steps  -  the  separation  of  the  bare  pseudo¬ 
potential  and  that  of  the  potential  due  to  screening.  The 
former  is  based  on  the  small  core  approximation  and  the 
latter,  on  the  approximation  of  linear  screening  (i.e.  on 
the  approximation  of  estimating  the  screening  by  the  first 
order  of  perturbation) .  Now  we  shall  briefly  describe  these 
two  steps. 

§2.3*1  Separation  of  the  Bare  Pseudopotential 

In  its  most  general  form  the  bare  pseudopotential  is 


given  by 


c  ,  'ifief  L  *io  r  >i  ixo  36  -io  ,19^51 
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Wb|J>  =  Vb  |  k>  +  Z  f ( K,a) | a>  <a | k>  ,  (2.19) 

ct 

where  | k>  is  an  arbitrary  plane  wave  with  vector  k  and  the 
superscript  b  stands  for  bare. 

First  let  us  show  the  separability  of  Wb  for  a  pure 
metal.  The  discussion  will  be  found  to  apply  to  an  alloy 
with  minor  modifications,  which  will  be  pointed  out  at 
the  end  of  §2.3.1. 

If  the  core  states  are  sufficiently  localized  for 
the  cores  of  the  adjacent  ions  not  to  overlap,  we  can 
write  Vb  in  a  pure  metal  as 


Vh(r)  =  Z .  vb  (r  -  r . ) 

J  J 


(2.20) 


j— 

where  vD(r  -  r. )  is  the  bare  potential  due  to  an  ion  seated 
at  r .  The  assumption  of  no  overlap  of  the  cores  also 
gives 


|  a>  =  ,  (r  -  r  .  )  ,  (2.21) 

L  J 

where  the  index  a  has  been  separated  into  the  pair  of 
indices  j  and  t;  j  as  before  labels  the  ions  and  t 

I 

denotes  the  energy  and  angular-momentum  numbers  for  the 
core  states  of  the  ion.  If,  now,  f(K,a)  depends  on  a 
only  through  t,  we  can  write 

Wb|K>  =  Z  wb  (r  -  r.)|K>  ,  (2.22) 

j  J 

with  wb  given  by 
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b  \  i^.  b  / -*  \ 

w  (r  -  r.)  k>  =  v  (r  -  r.)  e 
J  1  J 


ik  .  r 


/& 


+  E  f(K,t)i[;1.(r-r,)/Tp*(r»-r.)(ei'?-?l/n35)dr»  . 

t  ^  J  ^  J 


(2.23) 


(2.22)  and  (2.23)  represent  the  separation  of  the  total 
bare  pseudopotential  into  pseudopotentials  associated 
with  individual  ions. 

Now  let  us  examine  the  reduction  of  f(<,a)  into 
f(K,t).  Actually,  because  of  the  arbitrariness  of  f(ic,a) 
we  could  always  choose  for  it  the  form  f(i?,t).  But  in 
view  of  the  criterion  of  the  smoothness  of  the  wave 
function  we  need  to  look  into  this  question  more  closely. 

The  assumption  of  no  overlap  of  the  cores  leads  to 
(2.21)  which  expresses  the  fact  that  the  core  states  are 
localized  at  particular  ions.  Thus  must  be  equal  to 
E,  .  Since,  if  the  boundary  effects  are  neglected,  all 

r  j 

the  sites  in  a  pure  metal  are  equivalent,  E^ .  must  be 

J 

independent  of  j .  So  we  have 


E  =  E,  =  E,  .  (2.24) 

a  t  .  t 
J 

(2.21)  and  (2.24)  immediately  give  that  within  the  small 
core  approximation  f(K,a)  reduces  to  f(ic,t). 

It  may  be  pointed  out  that  the  small  core  approxima¬ 
tion  is  important  not  only  in  achieving  the  separation  of 
Wb ,  but  also  in  estimating  E^  and  For  a  sufficiently 

localized  core,  we  may  neglect  the  variation  over  the  core 


■  ' 
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of  the  potentials  due  to  the  conduction  electrons  and 
the  adjacent  ions.  It  follows  that  the  core  wave  func¬ 
tions  are  the  same  as  in  an  isolated  ion  though  the 
energies  of  the  cores  states  will  be  different.  The 
difference  between  the  energy  of  a  core  state  in  a  metal 
and  the  energy  of  the  same  state  for  an  isolated  ion  can 
be  calculated  by  the  Wigner  Seitz  approximation. 

Therefore  if  we  know  the  wave  functions  and  the 
energies  for  an  isolated  ion  (which  are  usually  obtained 
by  Hartree-Fock  calculations),  we  can  determine  $  and 
in  a  metal. 

j_ 

It  is  clear  that  the  separability  of  Wu  will  be 
carried  over  to  the  alloys  also,  since  the  only  modifica¬ 
tion  we  will  have  to  make  is  change  vb,  wb  and  t  to  vb , 

J 

wb  and  t.  respectively.  It  may  be  noted  that  within  the 
J  J 

Wigner  Seitz  approximation  E-t  .  depends  upon  j  only  through 

J 

the  specie  of  the  ion  occupying  the  site  j  and  as  a  result 
of  the  small-core  approximation,  the  same  clearly  applies 
to  \p+-  •  Therefore  for  a  given  ion,  not  only  vu  but  the 

j 

whole  of  wb  is  independent  of  the  location  of  the  ion  in 


the  crystal. 


§2.3.2  Linear  Screening  and  Separation  of  the  Potential 
due  to  Screening 

Now  we  shall  consider  the  linear  screening  of  the 

T_ 

bare  potential  W  .  Here  we  shall  directly  deal  with  the 
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alloys  (rather  than  first  doing  for  the  pure  metal  and 
then  generalizing,  as  we  did  in  §2.3*1)  because  here  the 
alloys  do  not  lead  to  any  extra  complication. 

We  shall  use  the  self-consistent  approach  followed 
by  Bardeen  (1937),  reformulated  in  terms  of  a  dielectric 
function  by  Lindhard  (195*0.  It  has  been  shown  by 
Ehrenreich  and  Cohen  (1959)  that  the  method  is  equivalent 
in  many  respects  to  the  many-body  treatment  using  the 
random  phase  approximation. 

In  the  absence  of  the  ionic  potential  the  density 
of  conduction  electrons  will  be  uniform  and  their  wave 
functions  will  be  plane  waves.  In  the  linear  approxima¬ 
tion,  the  evaluation  of  the  deviation  from  the  uniform 
distribution  is  based  on  the  use  of  the  zero  order  wave 
functions  corrected  to  the  first  order  for  the  effect  of 

the  ions.  For  an  electron  of  wave  number  k,  this  wave 

<  : 

function  is 

<J>£  =  |k>  +  I>a^(5)  |k  +  q>  ,  (2.25) 

where  the  prime  over  the  summation  sign  indicates  the 
exclusion  of  q=0  and  a^-(k)  is  given  by  the  usual  first 
order  expression 


a+(k) 

q 


<k+q  |  W  |  k>  _  <  k+q  [  W  |  k> _ 

E+  -  e2  +  (fi2/2ni)(k2- |  k+q  |  2  ) 

k  k+q  '  1 


(2.26) 


The  electron  density  corresponding  to  6+  is  given  by 
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^5^5  -  ( 1  / )  |_1  +  E->-aj(k)e^k‘r+  a+(k)e  ik,r 


']• 


(2.27) 


Thus  the  deviation  from  the  uniform  distribution  associa- 


ted  with  an  electron  of  wave  number  k  is 


>  .  -+ 

(l/fi)Z->  a->(£)e1  ,r  +  a->-(k)e 

q  q  q 


(2.28) 


which  we  rewrite  as 


(l/fi)Z+  fa->(S)  +  a*+(k)1  elq,r 
q  L  q  -q  J 


(2.29) 


Summing  it  up  over  all  the  occupied  states,  we  find  that 
N(r),  the  total  deviation  from  the  uniform  density,  is 
given  by 


N(r)  =  Z->  N(q)  eiq,r  , 


(2.30) 


where 


N(q)  =  h-  /  d£  [a^($)  +  &*-»(£)] 

4tt3  q  _q 

1  .  dk<k+q | W 1 k> 

2ttj  (n  /2m)(k  -  |  k+q 

Note  that  the  screening  charge  density  depends  upon  W  (not 
just  on  Wb),  which  requires  the  knowledge  of  the  screening 
charge  itself.  So  we  have  to  solve  the  problem  self- 
consistent  ly  . 

We  note  that  W  is  given  by 


2 


(2.31) 


w  =  wb  +  ws 


5 


(2.32) 


v;d  n^vis 
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g 

where  W  is  the  contribution  due  to  screening.  Further, 

W  (r)  is  related  to  N(r)  by  the  Poisson’s  equation 

V2Ws(r)  =  -47re2N(r)  , 

S  y  g  ^ 

so  that  W  (q),  the  Fourier  transform  of  W  (r),  is  given 

by 

WS(q)  =  N(q)  .  (2.33) 

q 


Substituting  (2.32)  and  (2.33)  in  (2.31)  and  noting  that 

g  ^  y 

W  (q)  is  independent  of  k,  we  get 


q2Ws  ( q )  _  1  j.  dS  <S+q  | Wb | _ 

47Te2  2 7T ^  (fi2/2m)(k2- |  £+q  |  2  ) 


<i2wS^)[,e(q)-l].  (2.34) 

^Tre 


where  e(q)  is  called  the  dielectric  function  and  is  given 

by 


£  ( q )  =  1  + 


2me  kp  r 

272~ 

Trq  h 


4kp-q2 

1  +  ■— o— ; -  In 


8qk 


F 


2kp+q 

2kp-q 


(2.35) 


Rearranging  (2.34),  we  get 


WS(q)  = 


2e‘ 


TT 


d£  <j£+q  [ Wb | k> 
2q2e(q)  '  (fi2/2m)  (i£- |  k+q  |  2  ) 


/ 


(2.36) 


Since  WD  has  earlier  been  shown  to  be  separable,  we  have 


b  it*  -iq.r, 


<k+q|Wb|k>  =  Z .  <£+q | w3 | £>e  '  j 

J  J 


3 


(2.37) 
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so  that 


Ws(q)  =  E . w? ( q ) e 

J  J 


s^^^-iq.r 


(2.38) 


with 


w^q)  =  2__2 


2g2  d5  <S+q  |  w1?  |  k> 


/ 


tt  q^e(q)  (  fi2/2m)  (k2- |  k+q  |  2  ) 


(2.39) 


(2.38)  immediately  gives 


Ws  (r ) 


S  /  -*■  s 

K.w.(r  -  r.)  . 


(2.40) 


where 


s  /->x 
w  (r ) 

J 


„  s  /->x  lq  .  r 
E+w. (q)e  H 
q  J 


(2.41) 


(2..40)  is  the  result  we  had  sought  to  establish. 

We  can  carry  the  argument  a  bit  further.  From  (2.38) 
and  (2.33)  we  can  write 


-* 


N(q)  =  E  .  n  .  ( q )  e~iq  *  rj 
J  J 


(2.42) 


where  n.(q)  is  related  to  w.(q)  in  the  same  manner  as  N(q) 
J  J 

to  W  (q).  From  Fourier  inversion  of  (2.42)  we  get 


N(r)  =  E.n.(r  -  r.) 

J  J  J 


(2.43) 


with 


-¥ 


n  .  ( r ) 

J 


I5n.(q)elq-rj 


(2.44) 


Thus  we  have  shown  that  if  the  bare  pseudopotential  is 
separable,  the  use  of  linear  screening  permits  us  to 
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divide  the  total  screening  charge  into  distributions  of 
charge  associated  with  the  ion  sites.  The  ion  and  the 
screening  charge  may  now  be  regarded  as  a  single  entity 
giving  rise  to  a  potential  w  given  by 


b  s 
w  =  w  +  w 


(2.45) 


Since  the  screening  charge  equals  the  charge  it  screens, 
this  entity  must  be  neutral  and  is  called  a  pseudo-atom . 

In  the  above  discussion  we  have  omitted  two  impor¬ 
tant  points. 

First  of  all  we  note  that  we  have  taken  the  zero 
order  energy  of  the  states  of  wave  vector  k  and  |  k+q> 

2  2  2  i  -*  i  2 

to  be  fi  k  /2m  and  fi  |  k+q  |  /2m  respectively.  Evidently, 
the  free  electron  mass  m  must  be  replaced  by  an  effective 
mass  m* . 


Secondly,  we  notice  that  we  have  neglected  the 
effect  of  exchange  on  the  screening  charge.  It  has 


been  shown  by  Sham  (1965)  that  this  effect  can  be  included 
by  multiplying  the  right  hand  side  of  (2.35)  and  (2.39)  by 
Ql-f  (q)]  ,  and  by  changing  e(q)  to  e*(q)  given  by 


e*(q)  =  1  +  -P-  -p-  - 

tt  q  fi 


with 


(2.47) 


^The  effect  of  correlation  is  not  significant  [Harrison 
(1966),  p.  323)] 
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where  k  is  half  the  square  of  the  Thomas  Fermi  screening 

o 

length  wherefore 


k 


2 

s 


1 

2 


mk. 


2*2 
t r  n 


2kp 

7T 


2 

me 

fi 


(2.48) 


It  is  obvious  that  while  both  of  these  points  are 
important  for  an  explicit  calculation  of  the  screening, 
they  do  not  affect  our  conclusion  regarding  the  separabi¬ 
lity  of  the  pseudopotential. 


§2.4  First  Order  Estimate  of  the  Scattering 
from  the  Screened  Pseudopotential  in 
terms  of  the  Form  Factors  of  Individual 
Pseudo-atoms 


In  the  first  order  of  perturbation,  the  probability 
of  scattering  from  a  state  with  vector  S:  to  another  with 
wave  vector  5+q  is  determined  by  the  matrix  element  of  W 
between  the  zero  order  states  | k>  and  | k+q> .  As  a  result 
of  the  separability  of  W,  this  matrix  element  can  be 
written  as 


<k+q|W|k>  =  E.<  k+q|w.|k>e  lc^*rj  .  (2.49) 

J  J 

By  defining  plane  wave  states  |K)  and  |£+q)  which  are  the 
same  as  |K>  and  |k+q>  respectively  except  that  the  former 
are  normalized  to  unity  over  the  volume  per  ion,  we 
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can  rewrite  (2.49)  as 


.  -*■  ->■ 

<E+q|W|E>  =  ( 1/N)  £  .  (E+q  |  w  .  |  E)  e-1C1  ’  r  j  , 

J  J 


(2.50) 


where  N  is  the  total  number  of  ions  in  the  crystal.  The 

matrix  elements  (E+q|w.|E)  are  called  the  form  factors 

J 

of  the  pseudo-atom  at  the  site  j .  The  calculation  of 
these  form  factors  is  of  interest  because  they  can  be 
used  to  deduce  a  host  of  atomic  and  electronic  properties 
of  metals  and  alloys.  That  they  are  related  to  the 
transport  properties  is  evident  from  (2.50).  Their 
relationship  with  other  properties  like  crystal  structure, 
elastic  constants  etc.  may  be  looked  up  in  Harrison  (1966), 
Chap.  6. 


For  a  pure  metal,  w.  is  the  same  for  all  j  and 

J 

(2.50)  takes  a  particularly  convenient  form 


<k+q|W|k>  =  (E+q I w I E)S(q)  , 


(2.51) 


where 

S(q)  =  (1/N)  I.  .  (2.52) 

J 

In  (2.51)  we  have  a  factorization  of  <E+q|W|E>,  which  in 
our  first  order  treatment  represents  the  scattering 
amplitude,  into  a  factor  (k+q|w|k)  which  depends  only  on 
the  properties  of  the  atoms  and  another  factor  S(q), 
called  the  structure  factor,  which  depends  only  on  the 
arrangement  of  the  atoms.  Such  a  factorization  of  the 


t 
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scattering  amplitude  is  a  direct  consequence  of  the  use 
of  the  first  order  theory  which  is  equivalent  to  the 
assumption  of  the  Born  approximation.  Below  it  will  be 
seen  that  the  result  is  carried  over  to  the  more  general 
case  of  the  weak  scattering  approximation. 

It  may  be  pointed  out  that  the  relation  (2.49) 
which  we  obtained  by  the  use  of  first  order  perturbation 
theory,  is  somewhat  more  general  and  may  be  written  as 

F(£,  £+q)  =  Z.fUfc,  i<+q)e“lq-rj 

J  J  j  \  2  •  bo 


where  F  is  the  total  scattering  amplitude  and ~  f .  is  the 
scattering  amplitude  associated ' with  the  site  j,  calcula¬ 
ted  exactly,  rather  than  by  the  first  order  theory.  It 
can  be  shown  [see,  e.g.,  Glauber  (1961)]  that  as  long  as 
the  total  potential  in  the  crystal  can  be  associated  with 
individual  ionic  sites,  all  we  need  to  obtain  (2.53)  is 
to  assume  the  validity  of  the  weak  scattering  approximation 
i.e.  to  assume  that  the  scattering  from  each  site  is  small, 
small  enough  for  multiple  scattering  to  be  neglected.  It 
is  not  at  all  necessary  to  calculate  the  scattering  from 
the  individual  sites  by  the  Born  approximation  as  we  did 
in  (2.49)  . 

It  is  important  to  note  an  inconsistency  in  the 


application  of  (2.53)  to  the  scattering  of  electrons  in 
metals.  It  is  clear  from  §2.3.2  that  unless  the  Born 
approximation  is  valid,  the  potential  due  to  screening 
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cannot  be  associated  with  the  ionic  sites  at  all.  In 
view  of  this  it  is  not  clear,  even  in  principle,  how  to 
go  consistently  beyond  the  first  order  theory. +  If  one 
is  not  satisfied  with  the  results  of  the  first  order 
theory,  the  only  option  is  to  include  higher  order  terms 
selectively.  This  is  precisely  what  one  does  [see 
Meyer,  Nester  and  Young  (1967)],  when  one  associates 

i  . 

a  potential  with  each  ion  site  and  calculate  the  total 

scattering  from  an  exact  estimate  of  f.  using  (2.53)  (or 

J 

its  counterparts  which  might  include  multiple  scattering 
effects ) . 

§2.5  Calculation  of  the  Form  Factors  of  a 
Pseudo-atom  in  a  Pure  Metal 

Now  we  shall  show  how  to  calculate  the  form  factors 
of  a  pseudo-atom  in  a  pure  metal.  For  the  sake  of  concrete¬ 
ness,  we  shall  take  the  Harrison  pseudopotential  as  the 


Note  that  but  for  the  conduction  electrons  the  total 
scattering  can  in  principle  be  calculated  exactly  (in 
practice,  the  only  difficulty  would  be  that  we  would  need 
higher  particle  correlations  to  take  into  account  multiple 
scattering  effects). 


. 

■ 
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basis  of  the  calculation.  For  any  other  pseudopotential 
the  method  is  similar. 

Let  us  recall  from  §2.1  that  for  the  Harrison 
pseudopotential  <k+q|W|k>  is  given  by 

<£+q  |  W  |  8>  =<  S+q|V|S>+I  ( k^+ <5  |  V  |  5>-E  )  4t+q  |  a>  <a  |  £> 

0 1  0 L 

+Z  (k2+<$|V|i{>-E  )  <S  I  a>  <a  I  £>  <k+q  I  P  I  k>  (2.54) 

01  a  i-<£|p|S> 

where,  for  convenience,  we  have  changed  to  what  we  shall 

call  the  semiatomic  units,  obtained  by  putting  Ti  =  2m  = 

2 

he  =  1.  In  this  system  the  unit  of  length  is  the  Bohr 
radius,  just  as  in  the  atomic  units,  and  that  of  energy 
the  Rydberg,  which  is  half  the  atomic  unit  of  energy. 

As  in  (2.51)  j  making  use  of  the  separability  of  W,  we  get 

<S+q|W|J>  =  (S+q | w | 5 ) S ( q )  ,  (2.55) 

with 

($+q|w|£)  =  (ic+q|  v  |£)+It  (k2+<£|  V|£>-Et )  (£+q|  t>  <t  |8) 

+  E,  (k2+  <£|V|k>-E.)  (£  1 1>  <t  |S)  5  (2.56) 

t  1- ( k | p | k ) 

where  p  is  given  by 

p  =  E 1 1 1 >  <t  (2.57) 

Evaluation  of  (2.56)  is  done  in  three  steps 
I.  Evaluation  of  (k+q|vb|k). 


;  1  :  ...  •  ' 
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^  y  TD  y 

II.  Evaluation  of  (k+q|w  |k).  This  consists  of  the 
second  and  the  third  terms  of  (2.56)  (the  superscript 
R  on  w  denotes  "repulsive"). 

h  R 

III.  Evaluation  of  the  screening  of  v  and  w  .  Note 

t)  V) 

that  since  v  consists  of  v°  and  the  screening  of  v  , 
w  should  not  have  anything  to  do  with  the  screening. 

o 

Why,  then,  does  the  expression  (2.39)  for  w  (q)  involve 
wb  (which  is  the  sum  of  vb  and  wR)  and  not  vb  alone? 

The  reason  is  obvious  i.e.  v13  being  a  deep  potential, 
the  calculation  of  its  screening  is  difficult.  So  we 
fall  back  upon  the  assumption  that  the  exact  screening 
of  v13  is  reasonably  well  represented  by  the  first  order 
screening  of  w*3 .  In  fact,  this  assumption  and  a  similar 
assumption,  namely,  the  exact  scattering  from  v  is 
reasonably  well  given  by  the  first  order  scattering  from 
w,  constitute  the  whole  spirit  of  the  pseudopotential 
technique . 

Finally,  we  add  I,  II  and  III  to  get  (£+q|w|S). 

The  potential  vb  consists  of  four  parts:  vj^, 

"U.  1-  l~v 

v(2)a  V(3)  and  V(4)5  which  denote,  respectively,  the 
potential  due  to  the  valence  charge  of  the  ion,  the 
electrostatic  potential  due  to  the  core,  the  potential 
due  to  conduction-core  exchange,  and  the  potential  due 
to  the  conduction  electrons  (for  reasons  that  will 
become  clear  later,  v^ ^  is  also  called  the  potential 
due  to  the  orthogonalization  charge).  If  one  is  more 
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ambitious,  one  might  include  the  potential  due  to  the 
conductron-core  correlation  also,  but  it  is  usually 
small . 

and  v^2)  can  be  calculated  by  noting  that 
p(r),  the  electron-density  associated  with  an  ion,  is 
given  by 

p (r )  =  p  (r)  -  (Z  +  n ) 6 ( r )  , 

where  p  (r)  is  the  density  of  electrons  in  the  core  and 
c 

(Z  +  n)  is  the  charge  on  nucleus  of  the  ion,  Z  being  the 
valence  of  the  ion  and  n  the  total  number  of  the  core 
electrons.  The  negative  sign  for  the  second  term  is  due 
to  the  fact  that  the  charge  on  the  nucleus  is  of  a  sign 
opposite  to  that  on  an  electron.  Using  Poisson’s  equation, 
(2.58)  gives 

(5+q|v^  2  |  k)  =  ^  ^  Dpc  ( Q )  “  2  -  ri]  ,  (2.59) 

5  o  q 

where  v^  2  denotes  the  sum  of  and  v^2  ^  and  Pc(5)  is 

given  by 

p  ( q )  =  /  e_iq’r  p  (r)dr  .  (2.60) 

Note  that  the  factor  2  in  (2.59)  arises  from  our  use  of 

Rydberg  as  the  unit  of  energy. 

p  (r)  may  be  calculated  from  the  wave  functions  of 
0 

the  core  electrons  since  we  have 


rij  sn±3on  ,\£d  b-  -J-.-'u 
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pc(r,e,*)  =  2Enj£jm(l/r2)p2Jl(r)Y;m(e^)YJlm(0,*)  ,  (2.6l)  + 

where  (l/r)P  (r)  denotes  the  radial  part  of  the  core  wave 
function  characterized  by  the  quantum  numbers  n  and  £. 
Noting  that 

^  » 


we  get 


Pc(r)  =  21 


l2-f.+  l)  i_  P2  (r) 
n ,  £  4tt  r2  nJT1  '  5 


so  that 


pc(q)  =  /  U(r)dr  , 


where 


(2.62) 


(2.63) 


U(r )  =  I  2(2J,+1)  P2a(r)  .  (2.64) 

To  obtain  v^)>  one  uses  the  Slater  approximation 
according  to  which  the  exchange  potential  due  to  like 
particles  of  density  d  is  given  by 

-3e2(3d/87r)1/3  (C.G.S.  units)  .  (2.65) 


Therefore,  the  potential  due  to  conduction-core  exchange 
in  rydbergs  would  be 


t 


The  factor  2  in  (2.6l)  arises  from  spin  multiplicity. 
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-6 ( 3pc/8tt ) 1//^  , 
whence 

(£+q|v^3)  |£)=  -(1.2707/q)(4Tr/no)/sin  qr  [rU(r )]  1/3dr ,  (2.66) 

where  1.2707  is  value  of  6  (  3/32tt^  )  -L//^  . 

Now,  regarding  the  conduction  electron  potential 
vj^,  let  us  remind  ourselves  that  as  yet  we  are  not  consi¬ 
dering  the  part  of  the  conduction  electron  potential  arising 
from  the  screening.  This  potential  arises  from  a  small 
nonuniformity,  present  irrespective  of  the  screening,  in 
the  distribution  of  the  conduction  electrons.  If  the 
electron-density  were  uniform,  there  would  be  no 
But  as  a  result  of  the  orthogonalizat ion  of  the  conduction 
and  the  core  states,  this  is  not  the  case.  The  difference 
may  be  described  by  superposing  on  the  uniform  distribution 
an  extra  charge,  localized  in  the  core  of  each  ion  in  the 
crystal.  This  charge  is  called  the  orthogonalizat ion 
charge.  Since  the  orthogonalizat ion  effectively  modifies 
the  charge  on  an  ion  and  since  the  orthogonalizat ion 
charge  must  be  compensated  by  suitably  modifying  the 
uniform  charge  densith  which  must  now  be  different  from 
what  it  would  have  been  in  the  absence  of  the  orthogonali- 
zation,  the  situation  is  the  same  as  if  the  actual  valence 


'  (SitS£\£)3  ;■  (lev  3 1  VO'  .  aneriw 
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+ 

Z  is  replaced  by  an  effective  valence  Z  .  The  uniform 

* 

density(of  Z  electrons  in  one  atomic  volume)  does  not 

contribute  anything  to  (£+q | vb | k) ,  but  the  orthogonali- 

¥; 

zation  charge  of  (Z  -  Z)|e|  per  ion,  does.  This  latter 

contribution  is  represented  by  vb^  .  To  determine  vb^ 

* 

we  need  to  know  Z  and  the  manner  in  which  the  charge 
(Z  -  Z  )  | e |  is  distributed. 

It  can  be  shown  [Harrison  (1966),  pp .  271-2]  that 
the  orthogonalizat ion  charge  per  valence  is  given  by 

«[<k|P|k>(l-<k|P|k>)_1]>>  ,  (2.67) 

where  <<  >>  denotes  the  average  of  the  enclosed  expression 

over  all  the  occupied  states.  Since  (2.67)  must  be  equal 

* 

to  (Z  -  Z)/Z,  we  get 

Z*/Z  =  «[l-<£  |P|E>[]_1>>  .  (2.68)* 

t 

Note  that  as  a  result  of  orthogonalizat ion ,  in  (2.58) 
(Z+n)  would  become  {Z  +  [n+(Z-Z  )J  } ,  which  is  the  same 
as  before.  This  is  as  it  should  be  because  orthogonaliza- 
tion  cannot  be  expected  to  change  the  charge  on  the  nucleus 

%  % 

Note  from  here  that  Z  >  Z  which  is  expected  since  ortho- 
gonalization  leads  to  an  exclusion  of  the  conduction 
electrons  from  the  region  of  the  core. 


j  J  :v  s  'f 

- 

*qq  < (dd£i  )  nco.  igH  .  nworis  e  l  . .  o:  j  - 

93Bt;9vb.  as3o no  .  «  >>  s'isriw 

si  t  \( 


: 0X9  HE  Oo • 3 b £ 9 

e  1<  '<  iiC  t-  1 .  .  •  -3 


4  5 


Since  for  a  pure  metal  <k  | P | k>  reduces  to  (k|p|k) 
and  since  we  have 

eik-r=  z  ^  [4tt(  2)1+1 )]  ^  ihJ_(kr)Y{o(6,<t»)  ,  (2.69) 

we  get 

Z*/Z  =  «  [1  -  I^A^OO]-1  » 

=  tV{1~Z  n.lAnhk)]2>~1  ^khk 

2 

j  4tt  k  dk 
o 

=  3  /  dx  x2  U-En>Jl|AnJl(xkF)]2}"1  ,  (2.70) 

where  is  given  by 

i/  i/  a 

AR£(k)  =  (  4tt/^o  )  2(  2  A+l )  2  /  r  j ^(kr )Pn£(r )dr  .  (2.71) 


The  distribution  of  the  orthogonalizat ion  charge 
can  be  shown  to  be  given  by  [Harrison  (1966)] 

<<  ( 1-  <lo  |  P  1 1?>  ) _1  {fi3*  elk,rE  i|)*<k|a> 

1  1  ora  1 

— r  — lk  .  r  1  ~y 
+  2  e  Z  ib  k> 

aTa  1 


■Z  ij>0<k|a><  3|k>}>> 

a ,  3  a  3  1 


(2.72) 


To  simplify  the  calculation  we  take  the  orthogonalizat ion 

t 

charge  to  have  the  same  distribution  as  the  core  charge 

and  in  view  of  (2.72),  this  is  not  expected  to  be  far  from 

b  * 

truth.  With  this  approximation,  v^^  is  (Z  -  Z)/n  times 


' 
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the  potential  due  to  the  core  charge 
comparison  with  (2.59)  we  get 


(£+q  |  v^2| )  |  £)  =  - 


b 


4tt  2 

ft  2 
o  q 


Z*Z 

n 


Pc(q) 


Therefore,  by 


(2.73) 


The  negative  sign  in  (2.73)  arises  from  a  difference  in 
the  sign  of  the  orthogonalization  and  the  core  charges. 

This  completes  the  evaluation  of  (S+q|v^|5).  But 
before  we  proceed  further,  we  note  a  point  which  will  be 
found  to  be  very  useful  later  on.  The  expressions  (2.59) > 
(2.66)  and  (2.73)  for  the  various  parts  of  (£+q|vb|k)  do 
not  contain  1 1.  Furthermore  they  depend  upon  q  only  through 
its  magnitude.  (The  first  is  a  consequence  of  the  fact 
that  v^  is  a  local  potential  and  the  second  arises  from 
v13  being  spherically  symmetric.) 

Now  let  us  see  how  to  determine  (£+q|w^|£).  We 
first  discuss  the  evaluation  of  (k+q,  1 1>  <t  |  k)  .  This  should 
be  straightforward  by  the  use  of  (2.69).  Allowance, 
however,  has  to  be  made  for  the  fact  that,  in  general, 

5  and  S+q  are  in  different  directions  and  since  the 
specification  of  (0,6)  in  (2.69)  is  in  polar  coordinates 
with  the  Z-axis  along  k,  the  expansion  of  ei(^k+cl)*r  wqn 
be  in  a  different  coordinate  system.  We  take  coordinates 
(61,61)  based  on  k  and  (02,62)  based  on  (k+q).  The  core 
wave  function  |t>  will  be  specified  in  coordinate  system 
of  k.  Thus  we  have 


.  o  :  -  !  'i 
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el^'r  =  I*  P*  (21+1)]*  i£jjl(kr)Y  o(01,01) 

Xj 

ei(k+q).r  _  2ji+i)]^  B  j  ^  (  |  K+q  |  r )  e2  ,  <J>  2  ) 

*t(?)  =  Vm(r)  =  (1/r)pnt(r)yei**l)  ' 


Noting  that  the  integrals  over  products  of  two  vanish 
unless  both  have  the  same  &-value  even  if  the  coordinate 
systems  are  rotated  with  respect  to  each  other,  we  get 

(5+q  1 1>  <t  | S) 

=  (  Wfto)(2£+l)  [/dr  r  j  £  (  |  S+q  |  r )  ( r ) /dftY^  q(  0  2  ,  <j>2  )  (  01 ,  )] 

x  [/dr  r  j  5/(kr)Pn£(r)-/'d^Y£m(e1^1)Y^(eia({)1)]  .  (2.74) 


The  final  angular  integration  gives  6mo,  which  means  we 
can  put  m=o  in  the  first  angular  integration  also.  The 
integral 


* 


fdSl  YJl0{e2’<|)2)Y-  (01  ) 


Zo 


l’Yl 


(2.75) 


may  be  evaluated  by  expressing  Y J.  0 ^  ,  (j> ^ )  in  terms  of 
' Y£m^  02  5  ^2 ^  the  formula 

Y£o(915<J)1)  =  Dm’°  (°  3  0  5  °)Y£r!ie25({)2)  a  (2.76) 

where  0  is  the  angle  between  the  axes  and  or,  in 

->■  ->■  ->* 

other  words,  it  is  the  angle  between  k  and  k+q.  By  the 
use  of  (2.76),  the  integral  (2.74)  becomes  D^Q(0,e,0)  which 


is  given  by 
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DqO(0,6 ,0)  =  Y£*(9,0)  =  P£(cose)  ,  (2.77) 

where  is  the  Legendre  Polynomial^ .  (2.74),  then, 

simplifies  to 


(5+q|t><£|^)  =  Ar£( |5+q| )An£(k)P^(cos6)  .  (2.78) 

To  complete  the  determination  of  w  ,  we  also  require 
the  value  of  <ic  |  V  |  lo-E^  .  For  this  we  have 


<k  V  k>-E  =  V  +  e  0  -A  n  , 
11  nl  1  n£  1  n£  5 


(2.79) 


where  V  is  the  average  potential  experienced  by  a  conduc¬ 
tion  electron  in  the  crystal,  e  is  the  core  energy  of  an 

n.  36 

(n, &)-electron  for  an  isolated  ion  and  A  „  is  the  extra 

potential  experienced  by  this  electron  when  the  ion  is  in 

the  crystal.  We  shall  treat  e  to  be  known  (it  can  be 

r*i  36 

determined  either  experimentally  or,  more  usually,  from 

a  Hartree-Fock  calculation) .  The  next  problem  is  to 

evaluate  V  and  A  „ .  Let  us  first  calculate  V. 

n& 

For  this  purpose  one  divides  the  crystal  into  Weigner 
Seitz  (W.S.)  cells.  Since  all  the  cells  are  equivalent, 

V  over  the  whole  crystal  is  the  same  as  V  over  any  one 


JL 

In  this  section,  some  symbols  have  been  used  for  more 
than  one  quantities  e.g.  note  the  use  of  P  in  P^  and  P^ 
This  was  done  only  in  order  to  stick  to  the  most  widely 
used  notation  in  each  case,  and  it  was  felt  that  it  will 
not  give  rise  to  any  confusion. 
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W.S.  cell.  In  this  cell  the  total  charge  distribution 
consists  of: 

(i)  a  nuclear  charge  (Z+n)|e| 

(ii)  a  core  charge  ne 

(iii)  an  orthogonalizat ion  charge  (Z  -Z)|e| 

(iv)  a  uniform  electron  density  of  total  charge  Z  e. 
Regarding  (ii)  we  recall  from  our  discussion  of  v13  that 
in  addition  to  an  electrostatic  potential,  it  also  gives 
rise  to  an  exchange  potential.  We  shall  denote  the 
electrostatic  and  the  exchange  constribut ions  as  (iia)  and 
(iib)  respectively. 

In  the  point-ion  approximation  (i),  (iia)  and  (iii) 
produce  the  same  potential  as  a  ( point ) charge  Z * | e |  .  The 
correction  to  (iia)  due  to  the  finite  size  of  the  core  is 

j_ 

just  v(2)^r^‘  The  Potential  due  to  (iii)  is  also  subject 
to  a  similar  correction  but  we  shall  neglect  this  since 
it  is  much  too  small. 

Now,  since  the  contribution  (iib)  is  clearly 
we  have 


3 


V  =  V^2)(q=o)  +  V  (*  3  )  ( q=0 )  +  Vrest  , 


(2.80) 


b 


where  the  first  two  terms  give  the  average  of  v^^r)  and 

v^o\(r),  and  V  ,  denotes  the  average  of  the  potential 
( j  j  re  St 

due  to  a  point  charge  Z*|e|  and  a  compensating  electronic 


charge  Z  e  distributed  uniformly  over  the  W.S.  cell. 
v^3)(q=°)  may  be  immediately  obtained  from  (2.66)  and  we  get 


SVSri  9.W 
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yl(3)(q-0)  =  ”^1*2707)  (^Tr/^o)/r  [rU(r)]  1//^>  dr.  (2.66)' 

tL 

For  evaluating  v^(q=o),  we  note  from  (2.59)  that 


v(2)(q=°)  -  “0  {"  iT  '  ^  [Pc(l)-n]} 

^  o  q 


_  Lt  4tt  2 


q-*o 


Q  2 
o  q 


[n  -  /  Sln-ar  u(r  )dr]  } 


Since 


/  U(r)dr  =  n  , 


we  immediately  obtain 


v(2)^q=0^  =  ^  ^2U(r)dr  .  (2.59)’ 

Assuming  the  W.S.  cell  to  be  a  sphere,  V  is  easily 

27  0  S  X/ 

shown  to  be  given  by 


V 


rest 


-0.6  Z*/r  . 

o 


(2.81) 


Here  rQ  is  the  atomic  sphere  radius  and  is  given  by 

r  =  (3ft  /^v)1/3  .  (2.82) 

o  o 

This  completes  the  evaluation  of  V.  Now  we  come  to 
the  calculation  of  A.  If  we  again  divide  the  crystal  into 
W.S.  cells,  the  only  contribution  to  the  potential  will 
come  from  the  cell  in  which  the  ion  under  study  is  situated. 
The  other  cells  being  neutral,  do  not  contribute.  In  the 
W.S.  cell  of  the  ion  we  again  have  the  various  charges  (i)  , 


<  n  *  Xu(t:)U  \ 
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(ii),  (111)  and  (Iv).  But  now  we  do  not  have  to  worry 
about  (1)  and  (ii)  because  their  effect  is  already 
contained  in  e  .  The  excess  potential  that  the  electron 
experiences  by  the  ion  being  in  the  crystal  compared  to 
the  case  of  the  isolated  ion  is  only  due  to  (iii)  and 
(iv).  Therefore,  An^m  (using  the  approximations  made 
earlier  in  this  section,  we  shall  show  presently  that 
depends  only  on  n  and  Z)  is  given  by 


n£m 


=  /[v°C(?)+vCC(r)]  [P^J,(r)/r2]Y*m(e,()))YI)m(8,<fi)dr  , 


Am 


Am 


(2.83) 


o  c  c  c 

where  v  and  v  are  respectively  the  potentials  due  to 
the  orthogonalizat ion  charge  and  the  conduction  charge. 

As  in  v10 ,  we  assume  the  orthogonalizat  ion  charge  to  be 
having  the  same  distribution  as  the  core  charge  and  as 
in  V,  we  take  the  W.S.  cell  to  be  a  sphere.  Then  we  have 
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(2.85) 


Since  v  (r)  varies  rather  slowly  with  r,  following 
c  c 

Harrison,  we  approximate  v  (r)  by  its  value  at  r=0. 

0  0 

However,  voc(r)  varies  much  more  rapidly  with  r  and  we 


J  -  . 
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cannot  make  the  same  approximation  in  this.  Thus  from 
(2.82),  we  have 


3Z 


n£ 


+  A 


o 


oc 

n£ 


(2.86) 


oc 


with  An^  ,  the  contribution  to  A^  from  the  orthogonali- 
zation  charge,  being  given  by 


An£  =  f  v°C^r^  Pn£('r')  dr  *  (2.87) 

(2. 86), with  (2.87)  and  (2.84)  determines  An^. 

This  finishes  the  evaluation  of  (k+q|w  |k).  It  is 
useful  to  note  that,  unlike  (S+q|vb|k),  it  depends  on  both 
S:  and  K+q  and  not  on  q  alone  (this  is  because,  as  noted 
in  §2.1,  wP  is  nonlocal).  However,  from  (2.56)  and(2.78) 
it  may  be  seen  that  this  dependence  reduces  to  the  dependence 
on  q,  k  and  |  S+q  | .  Thus  the  vector  triangle  of  k,  q  and 
5+q  may  be  rotated  in  space  without  affecting  the  value  of 
(k+q|w  |k).  This  in  fact  would  be  true,  see  Harrison  (1966) 
p.  20,  for  all  w  for  which  f(<,t)  depends  only  on  the 
magnitude  of  k  which  is,  of  course,  true  for  the  Harrison 
pseudopotential . 

Now  let  us  discuss  the  screening.  For  simplicity 
we  shall  discuss  Hartree  screening  and  will  take  the 
electron  mass  to  be  the  same  as  its  free  mass.  As  discussed 
in  §2.3.2,  the  expressions  thus  obtained  can  easily  be 
modified  to  include  the  effect  of  exchange  and  the  influence 
of  the  effective  mass. 
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Modifying  (2.39)  for  the  change  of  units  and 
remembering  that  now  our  plane  waves  are  normalized  over 
instead  of  ft,  we  get 


S  /->N 

w  (q) 


2  2  ,  , 
tt  q  e(q) 


rrdK(K+q 

b 

V 

|K)  ,  rdk(K+q 

i  R 
|  w 

k)7 

L  k2- 

k+q  | 

2  +  J  2  i 

k  - 

r*  ■+ 

k+q 

2 

(2.88) 


Both  the  integrals  in  (2.88)  contain  a  singularity  at 
k= | K+q | .  We  want  to  find  the  principal  values  of  these 
integrals.  The  first  term  does  not  present  any  difficulty 
since  (k+q | vb | £)  can  be  taken  outside  the  integration  sign 
and  the  rest  of  the  integral  can  be  evaluated  analytically. 
However,  the  same  does  not  apply  to  the  second  term  and 
we  have  to  evaluate  the  integral  numerically  taking  the 
singularity  into  account.  Let  us  postpone  our  worries  for 
the  time  being  and  just  denote  the  second  integral  by  I. 
This  gives 

wS(q)  =  vb(q)  +  -j-l -  I(q)  .  (2.89) 

£(qj  tt  q  £ ( q ) 

s 

Note  that  w  depends  only  on  the  magnitude  of  q  because 
from  the  foregoing  discussion  we  know  that 


(k+q|vb|k)  =  vb(q) 
and 

(k'|w^|k)  =  wR(q,k,k’)  , 

->-->■ 

where  k'  =  k  +  q. 
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Adding  the  matrix  element  of  wD  to  (2.88),  the  form  factor 


of  the  total  screened  pseudopotential  is  found  to  be 


it  q  e(q) 


Note  that  the  screened  form  factor  of  v  is  given  simply 
by  dividing  the  bare  form  factor  vu(q)  by  the  dielectric 


function  e(q).  This  is  a  standard  result  for  a  local 


R  R 

potential  and  will  not  apply  to  w  ,  since  w  is  nonlocal. 


However,  in  the  so-called  semilocal  approximation,  the 


nonlocal  nature  of  w  is  ignored  for  the  purposes  of  the 


calculation  of  the  screening.  In  that  case  we  get 


(2.91) 

semilocal  approximation. 


In  the  full  nonlocal  theory  we  cannot  use  (2.91)  and  must 
try  to  evaluate  the  principal  value  of  I(q)  and  determine 
the  form  factors  of  w  from  (2.90).  The  discussion  of  the 
evaluation  of  I(q)  is  given  in  Appendix  I. 

To  facilitate  the  writing  of  a  computer  program,  we 

summarize ' below  the  various  steps  in  the  evaluation  of  the 

pseudopotential  form  factors 

1.  Data 

ft  for  the  metal, 
o 

Z,  P  n(r)  and  e  n  for  the  ion.  (The  last  two  are 
5  n£  ni 

usually  obtained  from  the  Hartree-Fock  calculations.) 
Using  P  0,  obtain  U(r)  from  (2.64)  and  obtain  Z*/Z 

L  l  Xj 


from  (2.70)  and  (2.71) • 


2. 
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3.  (i) 

Using  (2.63),  obtain  p  (q)  from  U(r)  and 
substituting  for  Z*/Z,  obtain  v1?  p  h  from 

1  5  3  ^ 

(ii) 

b  4tT  2  r  ,  /  \  Vi  Z  *  — Z-I  , 

1,2,4  n  •  2  {  Z  n+PcU)[l  -  n  ]}, 

0  Q 

which  is  the  sum  of  (2.59)  and  (2.73). 

Using  (2.66),  obtain  v^(q)  from  U(r). 

(iii) 

Calculate  e*(q)  or  e(q)  from  (2.46)  or  (2.35) 

depending  on  whether  or  not  the  effect  of 

(iv) 

exchange  on  screening  is  to  be  included. 

Add  ^  ^(q)  and  v^(q),  and  divided  by  e*(q) 

or  e(q),  as  the  case  may  be,  to  get  v(q),  the 

screened  form  factor  of  v  . 

4.  (a) ( i ) 

Use  U(r)  in  (2.59)’  and  (2.66)’  to  get  vk(q=o) 
and  v^ (q=o ) . 

(ii) 

Substitute  these  in  (2.80)  to  get  V  with  the 

(b)  (i) 

help  of  (2.81)  and  (2.82) . 

Evaluate  voc(r)  from  (2.84). 

(ii) 

t 

0  c 

Substitute  in  (2.86)  to  get  A^  ,  and 

(iii) 

Use  (2.86)  to  get  A^. 

(c ) 

Obtain  |  V  |  lo>  — En^  from  (2.79). 

5.  Calculate  (k+q|t><t|k)  from  (2.71)  and  (2.78). 

6.  Substitute  <£  |  V  |  S>-En£  and  (k+q|txt|k)  in  the  last 
two  terms  of  (2.56)  to  get  w  . 

7.  Evaluate  I(q)  and  obtain  w(q,k,kT)  from  (2.90)  or 
use  (2.91)  if  one  does  not  mind  the  use  of  the  semi 
local  approximation. 
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It  may  be  mentioned  that  a  pseudopotential  calcula¬ 
tion  with  nonlocal  screening  takes  much  longer  as  compared 
to  a  similar  calculation  using  the  semilocal  approximation. 
The  actual  amount  of  extra  computer  time  depends  upon  the 
complexity  of  the  ion  under  consideration:  For  Li,  nonlocal 
screening  takes  3  times  longer  and  for  Zn  about  10  times 
longer.  Yet,  in  view  of  the  high  accuracy  demanded  of  the 
calculation  of  the  form  factors  to  be  used  for  obtaining 
the  resistivity  [see,  e.g..  Wiser  (1966)],  it  appears  to  be 
worthwhile  to  go  through  the  trouble  of  calculating  the 
screening  nonlocally. 

§2.6  Critique  of  the  Pseudopotential  Method 

Inspite  of  an  impressive  record  of  success  in  explain¬ 
ing  a  wide  range  of  the  properties  of  metals,  there  are 
certain  features  of  the  theory  of  pseudopotentials  which 
make  one  feel  uncomfortable  about  the  whole  thing.  One  is 
the  lack  of  a  suitable  criterion  for  selecting  a  pseudo¬ 
potential.  Should  one  try  to  have  the  smoothest  wave 
function  or  the  smoothest  potential  or  the  most  rapid 
convergence  of  the  Born  series?  Another  and  perhaps  by 
far  the  most  serious  one  is  that  while  with  a  suitably 
chosen  pseudopotential,  first  order  calculations  of  the 
scattering  may  be  reasonably  accurate,  it  is  not  clear  how 
to  proceed  to  higher  orders  consistently.  It  may  be  noted 


d.  I 
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that  if  we  knew  how  to  calculate  the  higher  orders,  the 
absence  of  a  criterion  for  choosing  would  not  have 
mattered  since  any  choice  would  give  the  same  final 
result  in  an  exact  calculation. 

But  then,  the  importance  of  the  method  of  pseudo¬ 
potentials  lies  not  in  that  it  leads  to  a  rigorous 
concept  but  in  that  it  is  a  very  fertile  approximation. 
Before  the  advent  of  pseudopotentials,  the  total  potential 
being  strong,  one  was  required  to  work  with  the  perturb¬ 
ing  potential,  and  for  every  new  situation,  a  new  postulate 
was  required  for  its  calculation.  Now  we  can  deal  directly 
with  the  total  potential,  the  scattering  due  to  which  is 
known  as  soon  as  we  have  calculated  the  pseudopotential 
of  the  ions  involved  and  have  determined  the  arrangement 
of  ions  by  some  means,  say  the  neutron  or  the  X-ray  scatter¬ 
ing.  The  technique  applies  no  matter  whether  we  want  the 
scattering  due  to  dislocations,  thermal  vibrations, 
impurities  and  what  have  you.  The  question  of  postulating 
the  perturbing  potentials  in  different  situations  just  does 


not  arise. 
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Chapter  3 

PSEUDOPOTENTIALS  IN  ALLOYS  AND  THE  CALCULATION 
OF  THE  RESIDUAL  RESISTIVITY* 

§3*1  Preliminary  Remarks 

As  we  stated  in  the  Introduction,  the  pseudopotential 
form  factors  of  an  ion  depend  on  its  environment.  In  this 
chapter  we  first  briefly  review  Harrison’s  (1966)  calcula¬ 
tions  for  the  form  factors  in  an  alloy.  Then  we  discuss 
the  improvements  that  are  necessary  and  calculate  the 
matrix  elements  afresh.  In  §3*3  we  use  the  matrix  elements 
to  calculate  the  residual  resistivity  of  a  dilute  alloy. 
We  shall  see  that  our  results  differ  substantially  from 
those  given  by  Harrison.  Further,  when  allowance  is  made 
for  the  effective  mass  factor  and  nonlocal  screening,  our 
calculated  results  are  in  satisfactory  agreement  with 
experiment . 

i 

Calculations  for  some  alloys,  not  considered  by 
Harrison,  are  also  presented. 

§3.2  Calculation  of  the  Form  Factors  in  an  Alloy 
§3.2.1  Harrison’s  Theory 

Harrison  suggested  that  the  form  factors  in  an  alloy 
can  be  calculated  by  applying  the  correction  for  the 

*  The  results  described  in  this  chapter  have  been  published 


in  Gupta  ( 1968 ) . 
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changes  in  the  dielectric  constant  e(q)  and  the  atomic 
volume  ft  .  Both  of  'these  are  simple  to  apply  and  we  get 
the  following  relation  between  the  form  factors  of  an  ion 
in  an  alloy  and  those  of  the  same  in  the  pure  metal: 


where  the  subscript  p.m.  stands  for  pure  metal. 


(3.1) 


If  one  wants  to  use  this  scheme  for  calculating  the 
form  factors  of  an  ion  for  a  number  of  different  conditions 
of  alloying  (achieved  either  by  changing  the  metal  with 
which  it  is  alloyed  or  by  changing  the  concentration),  the 
calculations  may  be  simplified  by  using  the  point-ion 
model  for  the  potential  of  the  ion.  In  this  model  the 
unscreened  pseudopotential  is  represented  by 


wD ( r )  = 


-  +  6  Sir) 


(3.2) 


where  Z  is  the  valence  of  the  ion  and  3  the  point-ion 
parameter  determined  by  adjusting  its  value  such  that 
the  screened  form  factors 


(  _  +  e)/noe(q)  (3.3) 

q 

agree  with  the  calculated  form  factors.  (As  in  chapter  2, 
the  units  we  are  using  here  are  the  semiatomic  units 
fi=2m=^e2  =  l )  .  Once  the  form  factors  for  the  pure  metal 

i  ’ 

have  been  calculated,  3  can  be  determined  and  the  form 
factors  for  any  situation  may,  then,  be  obtained  from  (3-3). 
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It  may  be  mentioned  that  the  results  obtained  from 
(3.3)  have  very  nearly  the  same  accuracy  as  those  obtained 
from  (3*1)  since  the  single  parameter  3  gives  a  rather 
good  fit  to  the  calculated  form  factors,  particularly  for 
the  set  of  form  factors  needed  in  resistivity  (i.e.,  q 
ranging  from  0  to  2kp  and  both  k  and  k*  equal  to  k^). 

Harrison  applied  this  theory  to  the  calculation  of 
the  form  factors  in  dilute  alloys.  In  such  alloys  the 
form  factors  of  the  host  may  be  taken  to  be  unchanged 
and  for  calculating  the  form  factors  of  the  impurity, 
[S10e(q)]alloy  may  be  taken  to  be  equal  to  [ftoe(q)]pure  host  • 
When  these  form  factors  were  used  for  determining  , 
agreement  with  experiment  was  found  to  be  poor. 

§3.2.2  A  Closer  Look  at  the  Form  Factors  of  an  Ion  in 
an  Alloy 

Compared  to  the  case  of  an  ion  in  the  pure  metal,  the 
ion  in  an  alloy  notices  three  important  changes  that  have 
a  bearing  on  its  form  factors: 

(a)  change  in  ft  , 

(b)  change  in  the  electron  density  surrounding  the  ion, 
and 

(c)  change  in  the  other  ions  among  which  the  ion  is 
placed . 

(a)  comes  into  the  picture  because  the  states  |k) 
and  |k+q)  characterizing  a  form  factor,  are  normalized 


.  .  .  )  •  ‘ i\  ,■  ai  s  :  *  •  flc 
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to  unity  over  ft  .  The  effect  of  this  change  is  easy  to 
take  into  account  and  it  is  clear  from  (3.1)  how  to  do 
it . 

The  effect  of  (b)  and  (c)  is  more  complex  in  nature. 

They  lead  to  a  shift  in  A,  the  zero  of  the  energy  in  the 

crystal  relative  to  the  free  ion,  to  a  change  in  and 

r 

to  a  change  in  the  orthogonalizat ion  charge.  Of  these, 
the  change  in  k^  is  brought  about  by  (b)  alone,  that  in 
the  orthogonalizat ion  charge  and  A  by  both  (b)  and  (c). 

All  these  three,  notably  A  and  kp ,  are  important  parameters 
in  a  pseudopotential  calculation. 

The  reason  why  A  enters  the  calculation  of  the  form 
factors  is  that  to  determine  the  scattering  from  an  ion 
in  the  crystal,  we  must  use  the  data  relevant  to  the  ion 
in  the  crystal.  Since  we  do  not  possess  these  data,  we 
deduce  them  from  the  corresponding  results  for  the  free 
ion.  As  discussed  in  chapter  2,  the  difference  between 
the  two  situations  may  be  described  to  a  good  approxima¬ 
tion  by  a  shift  A  in  the  core  energies  while  the  wave 

i 

function  may  be  taken  to  be  unchanged.  The  extent  of  the 
influence  of  A  on  the  form  factors  has  been  clearly  brought 
out  by  Lin  and  Phillips  (1965).  They  conclude  that  a  change 
in  a  by  0.1  ryd  changes  the  form  factors  by  0.01  ryd,  which 
is  the  accuracy  to  which  one  tries  to  determine  the  form 
factors.  (It  may  be  mentioned  that  the  shift  in  A  on  alloy¬ 
ing  is  usually  from  0.3  to  0.5  ryd.) 
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Regarding  how  the  form  factors  are  influenced  by 
kp,  one  point  is  obvious,  i.e.,  through  the  effect  of  kp 
on  the  potential  due  to  screening.  This,  of  course,  has 
been  taken  into  account  by  Harrison.  (It  hardly  needs 
to  be  mentioned  that  changing  e(q),  as  done  by  Harrison, 
takes  care  of  the  change  in  the  potential  due  to  screening 
only  in  the  semilocal  approximation. )  But  when  one  is 
interested  in  such  properties  (like  electrical  resistivity) 
as  depend  upon  the  matrix  elements  of  the  potential  with 
one  or  both  the  states  lying  on  the  Fermi  surface,  k^ 
enters  in  yet  another  way.  The  pseudopotential  being 
nonlocal,  the  matrix  elements  involve  not  only  q  but  also 
kp.  The  importance  of  taking  both  the  initial  and  the 
final  momentum  (and  not  just  the  transfer  of  momentum) 
into  account  for  the  calculation  of  pseudopotential  matrix 
elements  can  hardly  be  over-emphasized  in  view  of  the  work 
of  Animalu  (1965)  and  of  Animalu  and  Heine  (1965)  who 
obtained  considerably  improved  agreement  with  experiment 
over  the  earlier  treatments,  which  had  ignored  this  in  the 

i 

calculation  of  the  screening. 

Thus  when  A  and  kp  change,  as  they  do  on  alloying 
because  of  (b)  and  (c),  we  cannot  expect  the  new  form 
factors  to  be  related  to  the  previous  ones  in  a  simple 
manner.  A  proper  account  of  this  change  can  be  taken 
only  by  making  a  fresh  calculation  with  the  changed 
values  of  these  quantities.  The  effect  of  the  change 
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of  orthogonalization  charge  is  not  likely  to  be  very 
significant.  But  in  view  of  the  high  accuracy  in  the 
form  factors  demanded  for  the  resistivity  calculations 
[see  Wiser  (1966)],  it  is  advisable  to  include  this  also. 

It  may  be  parenthetically  pointed  out  that  the 

above  remarks  about  the  effects  of  A,  k^  etc.  on  the  form 

factors  are  by  no  means  restricted  to  the  pseudopotential 

approach.  They  apply  equally  well,  no  matter  which  of 

the  various  formulations  of  the  effective  potential  theory, 

+ 

the  pseudo-,  the  model  and  the  quasipotential,  one  chooses 
as  the  basis  of  the  calculation. 

f 

Below  we  shall  obtain  an  explicit  expression  for  the 
form  factors  of  an  ion  in  an  alloy  from  the  Harrison 
pseudopotential.  We  shall  consider  the  case  of  a  dilute 
alloy,  dilute  enough  for  the  impurity  to  be  regarded  as 
isolated.  As  mentioned  before,  in  this  situation  the  only 
problem  is  to  calculate  the  form  factors  of  the  impurity; 
those  of  the  host  may  be  taken  to  be  unchanged  by  alloying. 
Our  reason  for  using  the  Harrison  pseudopotential  as  the 
basis  of  the  calculations  is  that  by  comparing  our  results 
for  residual  resistivity  with  those  of  Harrison,  the 
importance  of  taking  into  account  the  changes  on  alloying 


^  The  dependence  of  the  model  potential  on  A  is,  perhaps, 
not  very  obvious  so  that  it  may  be  useful  to  mention  that 
it  arises  from  the  dependence  of  A^,  the  parameters  of  the 
model  potential,  on  A. 


r  {  3  3  J  >1  O :  t  1 1 9W  Y  [  1  r  p  9  q<3B  'l  tt‘ 
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in  A,  kp  etc.  in  the  calculation  of  the  form  factors, 
can  be  clearly  appreciated. 


§3*2.3  Impurity  Form  Factors  in  a  Dilute  Alloy  from 
the  Harrison  Pseudopotential 


The  matrix  element  of  Harrison's  optimized  and 
linearized  total  pseudopotential  W  is  given  (in  rydbergs) 
by  (2.54).  As  before  (§2.3),  on  writing  the  core  states 


.  ->  -> 


I a>  as  t.(r-r.).  where  t.  denotes  the  t-th  core  state  of 

the  ion  occupying  the  position  r.,  and  on  using  the  usual 

J 

assumptions  that  the  screening  is  linear  and  that  the  core 
energies  of  an  ion  depend  only  on  the  specie  of  the  ion 
(i.e.,  they  are  independent  of  the  location  of  the  ion  in 
the  crystal),  from  (2.54)  we  get 


->  -> 


<k+q|W|£>  =  N”1  E  .  e  iq,rj  (S+q|w.|£)  , 

J  J 


(3.4) 


with 


( k+q | w . | k )  =  ( K+q |  v  . | k ) 

<3  J 

+  Zt  [k2+<k|V|k>-Et  ]At  (|S+q|)At  (k)P£(cos0) 
^  j  j  J  j 

[k2  +  <5|V|5>-Et.]  [At  (k)]2 

+  {N  Ei  t  - J -  > 

j  1- <k  P  k> 


-> 


x  Zt  At  (|lc+q|)At  (k)P£(cos0)  , 
j  J  J 

where  A^  is  given  by  (2.70).  The  factor  P^(cos0) 
J 


(3.5) 

where 
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0  is  the  angle  between  k  and  k+q,  arises  (see  §2.5)  from 
the  projection  of  | t . >  on  two  plane  waves  |k)  and  |k+q) 
propagating  in  directions  which  are,  in  general,  different. 

We  notice  from  (3*5)  that,  in  general,  the  form 
factors  of  an  ion  in  an  alloy  depend  upon  the  nature  and 
concentration  of  the  constituents.  The  situation  simpli¬ 
fies  considerably  for  a  dilute  alloy.  For  this,  the  form 
factors  of  the  host  are  left  unchanged.  For  the  impurity 
we  get 

(k+q|w|k)  =  (k+q|v|k) 

+  Zt[k2+<5|V|k>-Et]  (  |  S+q  |  )  Aj.  ( k)  ( cos0  ) 

[k2  +  <S|v|K>-Et]  [At(k)]2 
Zt'  1  -  It„[At„(k)]2 

x  At ( I E+q I )At (k)P^( COS0 )  .  (3.6) 

The  method  of  evaluation  of  (3.6)  runs  parallel  to  the 
evaluation  of  (2.56)  for  the  pure  metal.  Here  we  shall 
go  over  the  various  steps  indicating  the  changes,  where 
necessary,  caused  by  alloying.  For  the  details  of  the 


Note  that  except  for  quantities  like  which  are 
obviously  for  the  host,  all  expressions  not  carrying  the 
subscript  host  are  to  be  evaluated  for  the  impurity  ion. 


I 


.  5 1  v  |  p+T )  =  '  I  p-»  - 

[(  o+A]  ra-<S|v|?  >+ 
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steps  involved,  one  may  refer  to  §2.5  (the  notation  used 
here  is  the  same  as  there). 


First  let  us  consider  the  part  (S+q|v|5)  of  (3.6). 


This  is  equal  to  the  sum  of  the  bare  potential  and  the 
potential  due  to  screening  of  the  sum  of  the  true  and 
the  repulsive  potentials.  The  bare  true  potential  and 


its  screening  add  up  to  give  vb(q)/c(q)  and  the  screening 

2  2 

of  the  repulsive  part  is  given  by  4  I(q)/Tr  q  e(q). 
Therefore,  we  get 


(5+q|v|£)  =  vb(q)/e(q)  +  4  I  ( q)/u2q2e  (q)  . 


(3.7) 


due  to  the  valence  charge  of  the  ion,  the  electrostatic 
potential  of  the  core,  the  potential  due  to  conduction- 
core  exchange  and  the  potential  due  to  the  orthogonaliza- 
tion  charge.  Since  the  Fourier  transforms  of  the  first 
three  depend  upon  alloying  only  through  ,  these  can  be 
evaluated  directly  using  (2.59)  and  (2.66).  Regarding 
vb^  ,  we  note  from  (2.67)  and  (2.72)  that  both  the  magnitude 
and  the  distribution  of  the  orthogonalizat ion  charge  per 
valence  charge  of  an  ion  depend  not  only  on  the  ion  itself 
but  on  (the  core  wave  functions  of)  all  the  ions  in  the 
crystal.  Therefore,  for  a  dilute  alloy  we  have 


V(4)(q)  =  (Z/Z 


host 


)[v(4)(q)]hOBt 


(3.8) 


.  '  li  £  •  ~  i  t 


(  p  '•  f,q  edd  r:  •  8U 
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For  the  evaluation  of  the  last  two  terms  of  (3.6) 
we  observe  that  the  average  potential  in  crystal 
<$|V|£>  will  be  the  same  as  in  the  absence  of  the 
impurity.  Therefore,  from  (2.80)  and  (2.8l)  we  get 

<£|V|5>  =  [-0.6Z*/ro  +  v(2)(°)  +  v(3) (°)]  host  *  (3-9) 

The  core  energies  of  the  impurity  ion  in  the  crystal  are 
given  by 


E,  E  E  n 
t  n  £ 


‘£njJ  +  An£  5 


where  A  0  is  given  by 
n  36 


AnZ=  3Z*/ro  +  /v°C(r)P 


For  obtaining  voc(r), 
for  obtaining  v^^(q) 


v 


oc 


<r>  =  <z/Zhost>£v 


oc 


(3.10) 

[see  eqns.  (2.85)  and  (2.86)] 

^(r)  dr  .  (3.11) 

we  use  the  same  argument  as  we  did 
above.  This  gives 

(r)]host  •  (3'12) 


In  view  of  the  above  changes,  the  outline  of  the 
complete  procedure  for  determining  the  impurity  form 
factors  in  a  dilute  alloy  would  be: 

T_  ‘U 

1.  Determine  v°  2  from  (2.59)  and  v^  from  (2.66)  using 
wave  functions  P  ^(r)  of  the  impurity  but  ftp  of  the 
host . 

2.  Calculate  [v°4 )  (q)]host  from  (2-^2)  and  thence  obtain 
v  (\  ) ( q )  for  the  impurity  from  (3.8). 


. 


1  )  W  £  /  '  L 
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3.  (i)  Calculate  f-0 . 6Z*/r  +v?-  >  (o)+v^,  v  (o  )1 .  . 

o  ( 2 )  ( 3 )  -  host 

and[voc(r)]host  . 

(ii)  Determine  C<£  |  V  |  io-Et7  from  (3.9)  to  (3.12). 

(iii)  Substitute  this  and  the  impurity  wave  functions 

in  (3*6)  to  obtain  w^(q,k,kf). 
b  R 

4.  Using  v  (q)  and  w  (q,k,k’)  obtained  above,  determine 
(ic+q|v|S)  from  (3.7),  where  I(q)  is  given  by  (2.88). 

4* 

Note  that  at  the  step  4  and,  in  most  cases  ,  also  at  3(iii), 

we  require  the  value  of  k^ .  The  k^  to  be  used  here  is 

r  r 

clearly  that  of  the  host. 

Further  details  of  the  above  steps  may  be  found  on 
pp.54-55  of  the  thesis.  For  the  sake  of  completeness,  a 
copy  of  the  matrix  elements  calculated  by  us  is  appended 
to  the  thesis. 

§3.3  Calculation  of  the  Residual  Resistivity 
and  Discussion  of  the  Numerical  Results 

Once  the  matrix  elements  are  known  the  calculation 
of  the  residual  resistivity  is  elementary.  The  probability 
of  scattering  P(i£,k')by  Nc  impurities  (N  total  number  of 


4- 

The  cases  for  which  3(iii)  involves  k^  have  been 
mentioned  before  i.e.  where,  as  in  the  calculation  of 
the  resistivity,  one  requires  the  value  of  (k+q|w|k)  with 
at  least  one  of  k  and  |£+q|  being  equal  to  kp. 
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atoms  in  the  crystal),  when  c  <<1,  is  given  by 
P(M')  =  |^[wA(q,k,k')  -  wB(q,k,k')]26(Ek-Ek,)  ,  (3.13) 


where  A  and  B  denote  the  host  and  the  impurity  atoms, 
and  q=k’-k.  Substituting  (3*13)  in  (1.8)  and  assuming 
E-k  relation  to  be  parabolic,  we  get 


k 


F 


^  ^OkFC  F  F  P 

- ~ —  f  (1-cose)  [w.  (q)  -  wfUq)]  sinedQ,  (3.14) 

2TThJ  A  b 


■Cl 

where  is  the  volume  per  atom  and  w  (q)  stands 
w(q,kp,kp).  From  (3.14)  and  (1.16),  we  get  for 
residual  resistivity  per  atomic  per  cent  of  the 
the  expression 

o  37T(m  ^  fio  2  f  F,  s  F /  N-i  2  3  n 

P „  =  - 9-  ",  -o  /  K(x)  ~  wR(x)J  x  dx  , 

400  e^hr  k^  o  A 
F 


for 

P°5  the 
impurity , 


(3.15) 


where  x=q/kp. 

Numerical  calculations  of  p°  were  undertaken  with 

r 

two  objects  in  view.  The  first  was  to  study  the  effect 
o  ' 

on  p^  due  only  to  the  modifications  of  the  impurity  form 
factors  discussed  in  §3.2.2  and  §3.2.3  because  while 
there  is  no  doubt  that  these  modifications  should  be 
taken  into  account,  it  is,  after  all,  not  very  meaningful 
to  go  through  all  the  trouble  involved  in  trying  to  include 
them  if  their  effect  is  going  to  be  small.  The  second 
object  was  to  see  if  by  incorporating  some  further  improve¬ 
ments,  it  is  possible  to  get  a  reasonable  agreement  with 
experimental  residual  resistivities. 


■ 


jj’ti  t  S£W 
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To  serve  the  f irst " obj ect ive ,  one  set  of  calcula¬ 
tions  was  performed  where,  while  the  modifications 
proposed  in  §3.2.2  and  §3.2.3  were  included,  the  treat¬ 
ment  was  otherwise  similar  to  that  of  Harrison  (local  screen¬ 
ing  and  m*=m) .  The  results  of  these  calculations 
are  given  in  column  A  of  the  "Present  *  calculations" 
in  Table  I.  The  alloys  treated  here  are  those  for  which 
Harrison  obtained  the  residual  resistivities  from  first 
principles.  The  difference  between* the ' result s  in  A 
and  those  of  Harrison  is  due  only  to  the  changes 
by  alloying  in  A,  k^  etc.  From  a  comparison  of  the 
results  it  is  clear  that  these  changes* exert  an  important 
influence  on  the  form  factors:  Taking  this  into  account 
in  the  calculation  of  the  form  factors  of  an  isolated 
impurity  can  change  the  residual  resistivity  by  as  much 
as  a  factor  of  4  or  more. 

After  the  importance  of  the  proposed  modifications 
was  established,  we  tried  to  improve ' the* set  A  from  two 
points  of  view:  first,  by  treating  the  screening 
nonlocally  including  the  effect  of  exchange  and  secondly. 


*  o 

Harrison  calculated  for  a  number  of  other  alloys 
also  by  adjusting  his  parameter  3  to  optimize  agreement 
with  experiment . 


V  £  t9  1-9'  3’  c  £  i  ’  ni 


' 
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Table  I.  Residual  resistivities  for  the  alloys  treated 
by  Harrison  from  first  principles.  All  the 
values  are  in  pft  cm. 


Alloy 

Host  purity 

Harri¬ 

son 

Present  calculations 

A  B  C  D 

Experi¬ 

ment 

A1 

Mg 

0.78 

0.17 

0.23 

0.23 

0.28 

0.33 

Zn 

0.1 

0.13 

0.14 

0.14 

0.17 

0.22 

Ca 

1.8 

3.1 

4,5 

4.1 

4.9 

>4.5 

Mg 

A1 

1.6 

0.63 

0.82 

0.83 

1.0 

2.1 

Li 

0.27 

0,33 

0.44 

0.58 

0.69 

0.75 

by  taking  into  account  the  influence  of  the' effective 
mass.  Column  B  in  Table  I  gives  the  results  incorporating 
the  first  improvement.  The  role  played  by  effective  mass 
has  two  aspects.  One  is  connected  with  its  influence 
through  the  screening  and  the  Fermi  energy.  The  calcula¬ 
tions  taking  this  into  account  are  given  in  C.  The  other 
aspect  relates  to  the  multiplication  of  the  resistivities 
obtained  by  using  the  f ree-electron  density  of  states, 

\L  p 

by  (m  /m)  .  The  set  D  represents  the  results  after  this 

multiplication.  The  effective  masses  used  for  obtaining 
the  sets  C  and  D  were  taken  from  the  work  of  Animalu  and 
Heine  ( 1965)  • 

Table  I  shows  that  while  the  single  most  important 
factor  to  which  the  poor  agreement  of  Harrison’s  results 
may  be  attributed  is  his  neglect  of  the  influence  of  the 
changes  in  A,  kp  etc.  for  the  impurity  ion,  the  agreement 


' 
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with  experiment  improves  as  successive  stages  of  refine¬ 
ments  are  incorporated.  It  is  gratifying ~ to" note  that 
except  for  the  case  of  A1  inMg,  the  set'D  reproduces 
the  experimental  resistivities  within  20%.  We  are  not 
clear  what  gives  rise  to  the  rather  large  discrepancy 
for  this  alloy.  It  is  possible  that  this  is  due  to  the 
occurrence  of  resonance  scattering  [[see,  e.g.,  Ziman 
(1964)J,  in  which  case  our  estimate  will  naturally  be 
much  lower  than  the  correct  value. 

We  may  draw  attention  tothe  fact' that  in  our 

treatment,  multiplication  by  the  *  somewhat ~ controversial 

%  2 

factor  (m  /m)  consistently  leads  to  a  better  agreement 
with  experiment;  see  discussion  on  p.io. 


Resistivities  of  dilute  alloys  of  Alkali  Metals 

j 

Recently  Dickey,  Meyer  and  Young  (1967)  have 
calculated  the  residual  resistivities  of  alloys  of 
monovalent  metals  containing  monovalent ■ impurities 
using  the  quasipotential  form  factors  obtained  from 
the  phase-shift  data  of  Meyer,  Nestor  and  Young  ( 1 967 ) 
(abbreviated  hereafter  as  MNY ) .  The  phase  shifts  used 
for  every  ion  are  relevant  to  the  appropriate  kp  so  that 
in  this  work  the  dependence  of  the  form  factors  on  kp 
is  automatically  taken  care  of.  Because  of  the  interest¬ 
ing  connection  between  the  influence  of  kp  on  the  quasi- 
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potential  form  factors  and  the  nonlocal  nature  of  the 
pseudopotential,  we  decided  to  perform  calculations  for 
alloys  made  from  the  combinations ' of  Li,  Na  and  K  to  see 
how  our  results  compare  with  those  of  Dickey  et  al. 


Table  II.  Residual  resistivity  for  the  alloys  obtained 
from  combinations  of  Li,  Na,  and  K.  All  the 
values  are  in  pfi  cm. 


Alloy 

Host  impurity 

Dickey 
et  al . 

Present 

B 

calculations 

C  D 

Experi¬ 

ment 

Li 

Na 

0.47 

0.02 

0.24 

0.43 

K 

5.69 

2.1 

0.79 

1.4 

Na 

Li 

0.33 

0.00 

■=j- 

0 

• 

0 

0 

• 

0 

K 

1.55 

2.0 

1.3 

1.6 

K 

Li 

0.50 

1.7 

3.7 

2.9 

Na 

0.67 

2.0 

3.4 

2.6 

The  results  of  these  calculations  are  given  in 
Table  II.  The  calculations  corresponding  to  column  A 
of  Table  I  have  been  omitted  here  and  those  given  in 
B,  C,  and  D  have  the  same  significance  as  before.  The 
effective  masses  used  in  C  and  D  are  taken  from  Ham  (1962). 
Since  Dickey  et  al.  have  not  considered  the  influence  of 
the  effective  mass,  their  results  should  be  compared  with 
our  set  B.  However,  our  best  estimate  is  represented  by 


the  set  D. 


' 
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From  Table  II  we  note  that' (a)  the  results  in  column 
B  do  not  compare  well  with  those  of  Dickey  et  al.  and  (b) 
except  in  the  case  of  two  alloys,  Na  in  Li  and  K  in  Na, 
our  best  estimates  also  differ  considerably  from  the 
results  of  Dickey  et  al. 

To  provide  at  least  a  partial  understanding  of  (a) 

we  would  like  to  draw  attention  to  certain ' result s 

obtained  by  Young,  Meyer  and  Kilby  ( 1967 ) .  They  found 

that  for  Li+  in  the  MNY  phase  shifts  at  the  k^  for 

* 

lithium  gave  a  rather  large  resistivity  for  the  liquid 

1 

metal  (93*3  y£2  cm  compared  to  the  experimental  value 
24.0).  Further,  they  found  that  the  agreement  is  much 
better  if  a  modified  potential  is  used  which  yields  the 
phase  shifts  0.434,  0.352,  0.015  and  0.001  instead  of  MNY 
phase  shifts  0.242,  0.417,  0.015  and  0.001.  On  using  the 
modified  phase  shifts  for  Li+,  we  get  0.06  and  4.8  yfi  cm 
for  p°  of  Na  in  Li  and  K  in  Li,  respectively  (note  that 
the  modified  phase  shifts  being  evaluated  at  the  k^  for 
lithium,  these  are  the  only  two  alloys  in  Table  II  for  which 

1 

they  can  be  used).  Keeping  in  mind  the  reliability  of  the 
second  decimal  place  in  these  calculations,  for  Na  in  Li 
the  new  value  may  be  seen  to  compare  well  with  the  value 


The  same  applies  to  Na  and  K,  but  to  a  lesser  extent. 
The  MNY  phase  shifts  give  17.0  and  27.5  yft  cm,  compared 
with  the  experimental  values  9.65  and  13.2. 


f. ;  =>c  -i  i J-  -•  !-,j-  - 
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in  column  B.  The  difference  remains  large  for  K  in  Li, 
though  the  gap  is  slightly  reduced. 

As  regards  (b),  experimental  work  is  necessary  to 
determine  which  results  are  closer  to  the' correct  value. 

It  may  be  mentioned  that  the ' agreement ' with  experiment 
obtained  by  Dickey  et  al . f or  the  alloy  Na  in  K  may  be 
fortuitous.  This  is  because  Na  and  K  atoms  have  a 
tendency  to  form  (Na2K)+  ion  and  this  would  affect  the 
resistivity . 

§3.4  Summary 

We  have  emphasized  the  role  of  the  surroundings  in 
the  form  factors  of  an  ion.  When  they  change,  as  they 
do  on  alloying,  Harrison's  recipe  is  too  inadequate  to 
take  into  account  the  effect  of  their  change  because  of 
the  accompanying  changes  in  the  Fermi  momentum,  the 
position  of  the  zero  of  the  energy  in  the"  crystal  relative 
to  that  for  the  free  ion,  and  the  orthogonalization  charge. 
These  can  be  taken  into  account  only  by  doing  the  whole 
calculation  all  over  again.  On  using  the  form  factors 
obtained  by  such  a  calculation  for  determining  the  residual 
resistivity,  we  are  led  to  a  considerably  better  agreement 
with  experiment  than  obtained  by  Harrison  except  for  the 
alloy  Al  in  Mg.  That  our  results  for  the  alloys  obtained 
from  combinations  of  Li,  Na,  and  K  differ  considerably 
from  those  of  Dickey  et  al.  should  not  be  disturbing,  in 


3  i  w  Jnoms  .s'ij}  ; 
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view  of  the  large  divergences  between  the  experimental 
values  for  the  resistivities  of  (pure)  liquid  Li,  Na, 
and  K  and  those  obtained  from  the  MNY  phase  shifts 
employed  by  Dickey  et  al .  It  will  be  of  interest  to 
have  the  experimental  results  for  these  alloys. 


■ 
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Chapter  4 


DEVIATIONS  FROM  MATTHIESSEN ’ S  RULE 


§4.1  Theoretical  Basis  of  Matthiessen ’ s  Rule 

and  Introductory  Remarks 

The  derivation  of  Matthiessen ' s  rule  employs  two 
assumptions  which  are : - 

A.  The  probability  P^  of  scattering  due  to  the  host 
lattice  and  P^  that  due  to  the  impurity  atoms,  considered 
at  rest,  are  additive  i.e.,  the  total  probability  of 
scattering  is  given  by 


P  =  P.  +  P 
i  r 

B.  The  trial  function  <J>g  in  (1.17)  for  the  alio 

taken  to  be  proportional  to  the  trial  function  <f>J 

k 

host  lattice.  In  other  words 


(4.1) 

y  is 
for  the 


=  a  =  b  >  (4.2) 

r 

where  <f>g  is  the  trial  function  when  P^  alone  is  present, 
and  a  and  b  are  constants  (independent  of  k).  (4.2)  is 
a  valid  assumption  if  relaxation  times  can  be  defined  for 
both  P^  and  P^  and  if  their  ratio  is  independent  of  the 
direction  of  it. 

The  total  resistance  p  of  a  dilute  alloy  is  usually 
written  as 


=  p .  +  p  +  A 
Ki  Kr 


P 


(4.3) 


' 

:  3"?  ;  .  fiw  sf  :  '•  •  . 

■ 

. 


si)',  a>  ti4-- 


fj&z  l  w 
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where  p.  and  p  are  the  ideal  and  the  residual  resisti- 
1  r 

vities,  A  is  the  deviation  from  Matthiessen ’ s  rule. 

These  deviations  would  occur  when  one  or  both  of  the 
assumptions  A  and  B  are  violated.  If  A  is  violated, 

A  can  have  either  sign  depending  upon  which  of  the  two 
sides  of  (4.1)  is  greater.  But  when  B  is  violated,  it 
can  be  seen  from  the  variational  principle  for  the 
resistivity  that  A  would  always  be  positive  [see,  e.g. 

Ziman  (i960),  pp.  285-286]. 

To  explain  the  observed  deviations  from  Matthiessen ’ s 

rule,  various  theories  have  been  proposed,  the  most  notable 

among  them  being  by  Dube  (1938),  Sondheimer  and  Wilson 

(1947)j  Koshino  (i960),  Kagan  and  Zhernov  (1966)  and 

Damon,  Mathur  and  Klemens  (1968).  The  first  two  are  based 

on  the  violation  of  B  and  the  other  three  on  that  of  A. 

* 

We  wish  to  point  out  here  a  cause  of  the  deviations 
from  Matthiessen ' s  rule  of  the  type  A  which  arises  from 
the  interference  between  the  scattering  due  to  the  thermal 
vibrations  of  the  host  atoms  and  the  excess  potential  due 
to  the  (thermally  vibrating)  impurities.  Such  interference 
terms  occur  because  it  is  the  matrix  elements,  rather  than 
the  probabilities  of  scattering,  which  are  additive,  and 
have  been  considered  in  different  contexts  before,  e.g.. 


The  results  described  in  this  Chapter  are  being 
published  in  Bhatia  and  Gupta  (1969). 


y  * 


. 
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in  Krishnan  and  Bhatia  (1944)  on  the  resistivity  of 
order-disorder  alloys,  in  neutron  scattering  formula¬ 
tions,  and  in  Faber  and  Ziman  (1965)  on  the  resistivity 
of  liquid  alloys.  Kagan  and  Zhernov  (1966)  have 
considered  the  interference  terms  in  connection  with 
the  deviations  from  Matthiessen ’ s  rule  also.  However, 
they  underestimate  significantly  the  importance  of 
these  terms  (see  footnote  on  p.  84  ,  and  discussion 
following  eqn.  (4.19));  moreover,  the  formulation 
given  here  is  different  from  theirs. 


§4.2  Calculation  of  the  Deviations 
from  Matthiessen ’ s  Rule 

i 

From  (1.17),  the  variational  expression  for  the 
resistivity,  we  have 


n 

p  "  j2 


(4.3) 


where 


11  =  2FT  /  )2  d^’  » 

B 


(4.4) 


and 

o 

9  f^ 

J  =  e  /  vg  4>g  9  E_^  dk  . 


(4.5) 


Recalling  that  Q( k,S ’ ) =f+( 1  -  fg t )P(£,£ 1 ) ,  and  using  time 
dependent  perturbation  theory,  we  have 


.  Oc.lB  •*!  •  '  "  *  3  '  HIO*  €-i  9b  i: 

. 
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Qdt.ic’  )=  ^r|<S'  I w  |ic> 1 2  f|(l-fS,)S(Ej-Ej>,+fiil))5 


where  fioo  is  the  energy  transferred  during  the  scattering 
Substituting  for  <k’  |w|k>  from  (2.50),  one  gets 


Q(  k,k'  )  = 


P  77  O  O 

fjd-fi,  )r(q,io) 


2 

fi  N 


k 


(4.6) 


where  f(q,u))  is  the  Van  Hove  scattering  correlation  func¬ 


tion  given  by 


f(q  ,co) 


N 


2ttN 


r  dt  -imt  v  f  \  f  \  -iq.r.(o)  iq.r.(t) 
J—  e  <  £  w.  (q)w.(q)e  H  l  e  ^  j  >, 

i,j=l  J 

(4.7) 


where  q=k'-k,  N  the  total  number  of  atoms  in  the  crystal, 

r.(t)  the  position  of  the  atom  j  at  time  t,  w.(q)  the 
J  J 

pseudopotential  matrix  element  of  the  atom  j  in  the  alloy 
and  <....>  denotes  ensemble  average.  Substituting  (4.6) 
in  (4.4)  and  restoring  the  factors  fiQ(the  crystal  volume ) , 

o  o 

1/4 it J  and  1/8tt  which  had  been  left  understood  before,  we 


get 


n  = 


ft 


o 


o 


o 


-> 


-¥  -* 


32ir5li2kBT 


/  ( .  )  ftfCl-ffct  )r(q,u)dk  dk' 


k 


k 


ft 


o 


32ir5f.2kBT 


/  f  t ) r((3 >“)dEjdEj i 


X 


dSS 


dS£, 


grad^E |  |grad+,E 


(4.8) 


(ft  =ft/N).  If  ho)  <<  Ep,  the  integrals  / dSg;  and  /dS+t  may 
be  taken  over  the  Fermi  surface  and  the  integrals  over 


. 
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Eg  and  Eg,  may  be  performed  independently  of  the  *  surface 

integrals.  To  evaluate  the  Eg  -  and  Eg,  -  integrals,  we 
o  o 

write  fg  and  fg,  explicitly,  define  eg=  Eg-  E^ ,  and  use 
Eg,  =  Eg  +  fio) .  This  gives 


n 


flo  dS£  dSS' 

32tt5Ti  F.S.  I  Srad£E  I  |gr-ad^,E 


(4.9) 


where  y(q)  is  given  by 


y(q)  = 


k  T 
B 


/  doa  r(q,m)  f 


de^- 

k 


(e  2,/kBT+l )  ( l+e“^  £K+flaJ^/kBT) 


"■Tico  /  if  T  ) 

The  integral  over  ef  gives  ¥iu)/(  1-e  7  B  ;  so  that, 


defining  $=fi/kgT,  we  get 


/ 


go) 

-Boo 


dw 


(4.10) 


Using  the  simplest  trial  function  (1.20),  and  assuming 
cubic  symmetry  and  spherical  Fermi  surface,  from  (4.3), 
(4.5)  and  (4.9),  we  get 


P  -  f  yav(x)x3dx  , 
o 


(4.11) 


where  x=q/k1-1 ,  ^=377  (m* ) 2  /(4h3e2k^),  and  y  (x)  is  the 

p  Or  civ 

average  of  y(q)  over  all  directions  of  q  in  the  crystal. 

(We  recall  that  k„  is  the  Fermi  wave  number  and  m*  the 

r 

effective  mass . ) 

Let  us  now  make  the  simplifying  assumptions  that 
(a)  the  alloy  is  of  the  substitutional  type  and  that  the 
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two  types  of  atoms  A  and  B  in  it  are  distributed  at 
random  and  (b)  A  and  B  atoms  have  roughly  the  same 
volume,  mass  and  force  constants,  so  that  the  effects 
of  local  strain  and  of  changes  in  vibration  frequency 
spectrum  can  be  ignored.  Then  in  (4.7)  we  can  replace 


wiwj  its  conf igurat ional  average  w^w^.  before  taking 
the  thermal  average.  Now 


w.w.  =  (w)2  +  6..  Tw2  -  (w)2l  . 

i  .1  i.i  L  J 


(4.12) 


If  c  is  the  concentration  of  B  atoms,  this  gives 


wiwj  =  [(l-c)wA  +  c  wfi]2+  61.  c(1-c)(wa-  wb)2  .  (4.13) 

Substituting  (4.13)  in  (4.7)  and  using  (4.10),  we  get 
for  y(q)  (which,  as  may  be  seen  from  (4.9)  or  (4.11), 
determines  the  resistivity),  the  expression 


Y(q)  = 


[(l-c)w.+c  wB]2  S ( q )  +  c(l-c)(w  -w  )2  S.  (q)  , 


(4.14) 


where 


S ( q )  =  [/ do)  3u)/(l-e  a))]/ 


-t5(A)N-i  r  dt  _.-iajtl  „-iq .  r  .  ( o )  iq . r  .  ( t ) 


2ttN 


'  V  -  -*  -> 

e  I  e" 

i  J 


(4.15) 

H>- 


o  Fr j  n  //n  -3ni\i  r  dt  -iwt  v  -iq.r.  (o)  iq.r .  (t). 

Sinc  U)=[/dw  geo/ (1-e  M  )J  e  <Ee^i  eMi  >. 

(4.16) 


Writing  S.  (q)  =  l+0(q),  we  get  from  (4.l4)  for  c  <<  1, 
i  in  o 


.  [  (5)  -  w]  a  +  s('  )  =  WjW 


. 

' 


. 
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y(q)  -  w2  S(q)  +  c(wA-  Wg)2  +  2w^ ( Wg-  wA)cS(q) 

+  c(wA-  wB)2e(q)  .  (4.17) 

The  first  term  in  (4.17)  is  responsible  for  the 

ideal  resistivity  of  the  host  metal  A,  the  second 

term  for  the  residual  resistivity  p  ,  and  the  remaining 

two  terms  give  rise  to  A.  One  verifies  that  at  high 

temperatures  (3m  «  1),  S.  (q)  =  l  and  hence  0(q)=O.  At 

n.  in  o 

low  temperatures  (T  <<  Debye  temperature  ©  )  9(q)  gives 
rise  to  a  term  in  the  resistivity  which  is  proportional 
to  T  Qsee  Taylor  (1964)J  and  can  dominate  over  the 
contribution  ( «  T  )  from  the  third  term  in  (4.17).  For 
T  >  O .  the  third  term  makes  the  dominant  contribution  to 
A.  This  is  the  interference  term,  referred  to  earlier  and 


Koshino  (I960),  referred  to  earlier,  ascribed  the 
deviations  from  Matthiessen ' s  rule ■ to ' a  term  similar  to 
the  last  term  in  (4.17),  but  with  0(q)  replaced  by 
essentially  S(q)  and  thus  wrongly  estimated  its  importance, 
specially  at  high  temperatures.  The  error  in  Koshino Ts 
work  has  been  pointed  out  by  Taylor  (1962,  64). 


- 

39Vi:i  (  p  )  6  (  O  9<d»a  >>  )  Z9?otBH>qm3t  wol 
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we  confine  further  discussion  to  it  only  .  Denoting  its 

contribution  to  A  by  A',  we  have  for  A1  and  p 

r 

2 

A’  =  2  c  /  wA(wg  -  wA)  Sav(x)x3dx  ,  (4.18) 

_  2 

p^  =  «$  c  /  ( wA  -  wg )  x J  dx  .  (4.19) 

First  we  note  that  because  of  Umklapp  collisions 
S  (x)  is  strongly  peaked  in  the  backward  directions 

ci  V 

(x  =!  2).  For  example,  whereas '  Sav(o)  (=nkBT3c,  where 
n=N/^  and  3  is  the  isothermal  compressibility) 

-  10  S  (2)  -  0.4  at  273°K  for  sodium;  at  40°K 
Sav(2)  -  0.04  £see  Greene  and  Kohn '  ( 1965  )J  .  Thus  for 
T  >  Q  ,  we  may  expect  from  ( 4 . 1 8 )  and  (4.19)  that 
A'/p  will,  in  general,  be  of  the  order  of  a  few  percent 
or  more,  which  is  of  the  correct  order ' of ' magnitude . 

Also  since  at  these  temperatures  Sav(x)  «  T,  we  have 
A’  a  T,  as  is  actually  observed. 


* 

It  may  be  of  interest  to  mention  that  the  theory  of 
Faber  and  Ziman  (1965)  for  the  resistivity  of  liquid  alloys, 
referred  to  earlier,  contains  just  this  type  of  interference 
term.  The  formulations  of  Krishnan  and  Bhatia  (1944)  and 
of  Kagan  and  Zhernov  (1966)  are  such  that  the  interference 
term  in  these  works  is  essentially  proportional  to  S(o). 


>-b  •'  )  \  = 
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The  interesting  point  about  the  expression  ( 4 . 1 8 ) 
is  the  combination  WA(WB~  w^)  in  which^the  pseudopoten¬ 
tial  matrix  elements  and  wg  occur  in  it.  This  is  in 
contrast  to  the  expressions  for  A  given  in  previous 
works  [Koshino  (I960),  Taylor  (1964),  and  Damon,  Mathur 
and  Klemens  (1968)J  where  w^  and  wg  essentially  occur 
in  the  combination  (w^-  Wg)^.  A  consequence  of  (4.18) 
is  that  A  can  be  of  either  sign.  Note  that  because  of 

o 

the  factor  S  (x)x  ,  the  sign  of  A  would  largely  depend 

on  the  sign  of  w^(wg-  w^)  in  the  neighbourhood  of  x=2. 

A  rough  calculation  using  the  pseudopotential  matrix 

elements  w^  ,  w^  obtained  by  Gupta  (1968)  for  Mg  in 

A1  and  vice  versa  show  that  A  should  be  positive  for  Mg 

in  A1  and  negative  for  the  complementary  alloy  (A1  in  Mg). 

Measurements  undertaken  by  Seth  (1967)  at  our  suggestion 

in  1965  show  such  a  difference  in  sign  for  the  complemen¬ 
ts 

tary  alloys  of  Al-Mg;  Seth  and  Woods  (private  communi¬ 
cation)  have  also  found  similar  results  for  Ag-Cd,  Ag-Mg 
and  Ag-Al  alloys. 

To  conclude,  we  have  shown  that  the  expression  (4.18) 
makes  a  significant  contribution  to  A  at  high  temperatures, 
although  for  numerical  agreement  effects  of  local  strain, 
deformation  of  vibration  spectrum  due  to  impurities,  and 


I  would  like  to  thank  Mr.  R.S.  Seth  and  Professor  S.B. 
Woods  for  letting  us  see  their  results  prior  to  publication. 
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the  anisotropy  of  the  Fermi  surface,  would,  in  general, 
have  to  be  taken  into  consideration. 


PART  III 


RESISTIVITY  DUE  TO  DISLOCATIONS 
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Chapter  5 


INTRODUCTION  AND  HISTORIAN,  BACKGROUND 

§5.1  A  Brief  Survey  of  the  Experimental  Work 

The  electrical  resistance  of  a  metal  is  found  to 
change,  usually  increase,  by  a  few  per  cent  when  the 
metal  is  plastically  deformed.  Change  is  also  observed 

4- 

in  other  electrical  quantities  like  thermoelectric  power 
Hall  coefficient,  m^gneboresistivity  etc.  The  change 
in  electrical  resistivity  was  first  observed  by  Bergmann 
(1889)  and  later  by  Coke  (1890).  More  systematic  inves¬ 
tigations  have  been  reported  by  Tammon  and  Dreyer  (1933), 
Dahl  (1936),  and  Smart,  Smit  and  Phillips  (19^1). 

The  association  of  this  extra  resistance,  6p  , 
with  dislocations  was  first  suggested  by  Koehler  (19^9). 
Various  studies  of  work  hardening  [see,  e.g.,  Cottrell 
(1953)J  undertaken  in  the  wake  of  the  discovery  in  193^ 


& 

In  some  cases  the  change  in  resistivity  on  plastic 
deformation  is  very  large.  Geiss  and  van  Liempt  (1927) 
observed  a  change  of  50%  in  tungsten  and  l8%  in  molybdnum 

4* 

This,  incidentally,  was  the  first  electrical  property 
for  which  a  change  on  plastic  deformation  of  the  metal 
was  observed  [Thomson  (1856)]. 
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of  the  connection  between  dislocations  and  plastic 

crystals  revealed  that  dislocations  get  multiplied 

during  cold  work.  In  well-annealed  crystals,  the 

discloation  density  is  of  the  order  of  10-  10°  cm 

11  12  —2 

and  on  cold  work  it  increases  to  10-10  cm  .  Thus 

the  extra  resistivity  6p  may  be  attributed  to  the 

increase  in  the  dislocation  density. 

However,  as  a  result  of  the  experiments  of 

Molenaar  and  Aarts  (1950)  and  of ' Broom  ( 1952 ) ,  it 

soon  became  clear  that  there  are  two  other  defects, 

* 

namely  vacancies  and  stacking  faults  which  contribute 

to  5p.  Let  us  denote  the  various  contributions  to  6p 

by  <$PD5  <5py  and  8p^,  where  the  subscript s- are  self- 

explanatory.  It  is  found  that  8p^  and  5py  are  roughly 

of  equal  magnitude  while  6p  is  usually  small  but  in 

s 

some  cases  (e.g„,  in  noble  .  jnatals  and  in  alloys  such 
Cu  +  Zn  which  have  low  stacking  fault  energies),  6p 

s 

is  about  10  -  20%  of  6p. 


* 

Stacking  faults  are  produced  by  plastic  deformation 
of  an  fee  crystal  because  dislocations  in  such  a  crystal 
may  dissociate  into  partial  dislocations  separated  by  a 
ribbon  of  stacking  faults. 


89 


Experimental  estimates  of  the  resistivity  of 
dislocations  are  usually  expressed  as  Sp^/N^,  which 
we  shall  denote  by  p^.  At  present,  it  is  not  possible 
to  estimate  pD  very  accurately;  the  results  are 
reliable,  at  best,  within  a  factor  of  2.  The  relia¬ 
bility  of  measurements  is  limited  primarily  by  the 
uncertainty  in  the  estimation  of  the  dislocation 
density.  The  uncertainty  in  the  separation  of  the 

measured  resistance  into  its  components  is  another 

* 

important  source  of  error. 

In  Table  III  we  give  the  values  of  pD  taken 
from  different  sources.  The  spread  between  the  various 
measurements  is  obvious  at  a  glance. 


Basinski  et  al .  (1963)  have  drawn  attention  to  yet 
another  source  of  error  in  the  estimation  of  p^, 
arising  from  the  possibility  of  a  rearrangement  of  the 
impurity  ions  when  the  crystal  is  plastically  deformed. 


d  nr  o  o  idl  oct  .  .  or  > 
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Table  III.  Values  of  in  10  ^  ohm-cm^  for  Cu  and 

A1  as  reported  by  different  workers. 


Cu 


PD 

Source 

PD 

Source 

2.3(1t.2°K)  + 

Blewitt  et  al . 
(1955) 

3-3  (80°K)+ 

Clarebrough 
et  al.  (1961) 

3 . 0( 80°K) 

Clarebrough  et 
al.  (1961) 

7  ±  2 ( 80°K ) 

Cotterill ( 1 9  6  3 ) 

1  ( 4 . 2°K ) 

Basinski  and 
Saimot 0 

(1967) 

6.4  ( 80°K ) 

Yoshida  et  al. 

(1963) 

1.3(4.2°K) 

1 . 6 ( 77°K) 

Rider  and  Poxon 

(1967) 

2 . 9±0 . 4  ( 80°K) 

1.8±0.1(4.2qK) 

Rider  and  Poxon 

(1966) 

A1 


+ 

This  indicates  the  temperature  at  which  the  measurement 
was  made. 


§5.2  Theory  of  the  Dislocation  Resistivity 


§5.2.1  Previous  Theories 

The  theory  of  the  resistivity  due  to  dislocations 
is  an  extremely  controversial  topic.  There  is  a  large 
discrepancy  between  the  experimental  and  the  theoretical 
values  of  -  the  latter  being  lower  by  a  factor  of  the 
order  of  10,  usually  more.  Even  though  the  experimental 
data  are  not  quite  reliable,  it  is  unlikely  that  they  are 
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systematically  In  error  by  such  a  large  factor.  It  is 
generally  agreed  that  the  theory  needs  to  be  Improved 
but  there  Is  general  disagreement' about  the  source  of 
the  discrepancy.  A  number  of  possibilities  have  been 
suggested  but  the  extent  of  their  influence  is  far  from 
clear . 

Three  different  approaches,  called  the  rigid  ion 
approach,  the  deformation  potential,  and  the  theories 
based  on  approximate  models  of  the  core  of  the  disloca¬ 
tion,  have  been  tried  for  the  calculation  of  the 
probability  of  scattering  (which  forms  the  basis  of  a 
resistivity  calculation)  in  the  presence  of  dislocations. 

The  earliest  of  these  was  the  rigid  ion  approach 
which  was  used  by  Koehler  (19^9)  in  his  pioneering  work 
on  the  calculation  of  for  edge  dislocations In  this 
approach  the  screening  charge  is  assumed  to  move  rigidly 
with  the  ions  as  they  are  displaced  from  their  positions 
in  the  perfect  lattice,  and  the  perturbing  potential  is 
given,  in  the  first  approximation,  by  the  product  of  the 
displacements  of  the  ions  and  the  gradient  of  the  ion 
potential.  The  difficulty  caused  by  the  displacements 
due  to  a  single  edge  dislocation  being  infinite  at  large 
distances,  was  avoided  by  considering  a  pair  of  edge 
dislocations  of  opposite  sign. 

Koehler’s  work  suffered  from  two  objections.  The 
first  was  that  to  solve  the  Boltzmann  equation,  Koehler 
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proceeded  by  expanding  the  distribution  function  as  a 
Fourier  series  and  retaining  only  the  first  few  terms; 
the  convergence  of  this  procedure  was  not  investigated. 

To  avoid  this,  Mackenzie  and  Sondheimer  (1950)  suggested 
a  procedure  equivalent  to  the  use  of  the  simplest  trial 
function  (1.20)  in  a  variational  calculation  of  the 
resistivity.  The  reliability  of  the  procedure  was  not 
clear  in  this  case  also,  but  it  is  expected  to  be  a  good 
approximation  when  the  resistivity  produced  by  disloca¬ 
tions  in  a  very  small  fraction  of  the  total  resistivity. 
Moreover,  the  method  is  much  simpler  than  that  employed 
by  Koehler  and  has  been  used  in  several  subsequent 
investigations  of  the  dislocation  resistivity.  The  second 
objection  against  Koehler’s  work  was  the  displacement  in 
the  presence  of  a  pair  of  edge  dislocations  of  opposite 
sign,  while  finite  everywhere,  is  discontinuous  across 
the  line  joining  the  dislocations;  and  this  invalidates 
the  procedure  for  evaluating  the  perturbing  potential. 

To  deal  with  this  difficulty,  Landauer  (1951)  introduced 
the  deformation  potential  approach  in  which  the  perturb- 
potential  depends  on  the  strain  (or  the  deformation) 
which,  unlike  the  displacements,  is  finite  and  continuous. 

In  the  deformation  potential  approach  the  screen¬ 
ing  charge  is  treated  in  a  basically  different  manner. 
Instead  of  requiring  this  charge  to  move  rigidly  with  the 
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ions,  here  one  permits  the  electrons  to  redistribute 
themselves  in  the  strained  lattice.  The  tendency  for 
the  electrons  to  redistribute  themselves  in  the  presence 
of  the- strain  arises  because "the 'Fermi "energy  Ep  depends 
upon  the  strain.  If  the  strain  is  inhomogeneous,  Ep 
varies  from  place  to  place  in  the  crystal  and  electrons 
flow  so  as  to  make  it  uniform.  The  perturbing  potential 
is  the  potential  caused  by  the  redistribution  of  charge. 

It  was  shown  by  Dexter  ( 1952 )'  [see  also  Harrison  (1958a)] 
that  this  method  of  evaluating  the  perturbing  potential 
is  justified  if  the  strains  vary  slowly. 

The  treatment  of  Landauer  was  based  on  the  nearly 
free  electron  (NPE)  approximation.  This  was  later 
generalized  by  Hunter  and  Nabarro  (1953).  In  the  NPE 
approximation  E^ ,  hence  the  perturbing  potential,  depends 
only  on  the  dilatational  strain.  However  the  general 
treatment ■ of  Hunter' and  Nabarro  showed  that  a  shear  strain 
also  gives  rise  to  a  perturbing  potential.  This  has  an 
important  bearing  on  the  resistivity  of  screw  dislocations 
whose  strain  field  is  one  of  pure  shear  in  the  approximation 
of  the  linear  elastic  theory. 

The  calculations  of  Hunter  and  Nabarro  gave  a  value 
-21  3 

6  x  10  fi-cm  for  of  edge  dislocations  in  Cu.  This  may 

be  compared  with  the  most  recent  experimental  estimate 

-19  3 

(see  Table  III)  of  1  x  10  £!-cm  .  The  experimental  value 
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of  pD  is,  of  course,  based  on  an  estimate  of  the  dis¬ 
location  density  which  includes- both  the  edge  and 
the  screw  dislocations,  but  this  only  makes  the 
discrepancy  more  serious  since  the  resistivity  due 

I 

■ 

to  screw  dislocations  is  expected  to  be  small 
[Harrison  (1958a)]. 

The  calculation  of  Hunter  and  Nabarro  has  been 
criticized  by  Carruthers  (1959)  and  by  Seeger  and 
Bross  (i960),  but  the  objections  do  not' appear  to  be 
well  founded  [see  Nabarro  (1968)  p.  630]  .  However, 
there  is  one  weakness  of  the  theory,  namely  the 
inadequacy  of  its  treatment  of  the  scattering  due 
to  the  core  of  the  dislocation,  which  had  been  known 
from  the  very  beginning  but  it  has  not  been  clear  what 
to  do  about  it.  The  effect  of  the  core  had  been  treated 
only  as  an  extrapolation  of  the  linear  strain  field. 

This  is,  of  course,  incorrect  because  of  the  severe 
distortion  of  the  lattice  at  the  core.  But  a  proper 
treatment  of  the  core  is  beset  with  two  difficulties  - 
evaluation  of  the  strain  in  the  region  of  the  core  and 
the  evaluation  of  the  perturbing  potential  due  to  this 
strain  (the  deformation  potential  would  not  be  adequate 
for  this  purpose  because  the  strain  in  the  core  is  not 
slowly  varying). 

To  provide  some  estimate  of  the  scattering  due  to 
the  core,  Harrison  (1958)  suggested  a  model,  called  the 
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hollow  cylinder  model,  for  the  core  of  the  dislocation. 

The  model  is  questionable  in  detail  but  has  the  essen¬ 
tial  features  of  the  core  i.e.,  the  interaction  is 
relatively  strong  and  is  restricted  to  a  region  of 
atomic  dimensions.  Considering  this  cylinder  to  be 

a  row  of  vacancies  and  ignoring  the  linear  strain  field 

-20  1 

altogether,  Harrison  found  for  Cu  to  be  5  x  10  fi-cm  . 

This  is  about  10  times  the  value  obtained  by  Hunter  and 

Nabarro  and  compares  favourably  with  the  experimental 

-19  1 

value  of  1  x  10  fl-cm  .  Calculations  have  also  been 
made  on  similar  lines  by  Basinski  et  al .  (1963)  and 
Ziman  (1964)  with  the  difference  that  their  estimate 
of  the  dislocation  resistivity  was  based  on  the  mean 
square  displacement  of  the  ions  in  the  presence  of  a 
dislocation.  Estimating  the  value  of  this  mean  square 
displacement,  from 'plausibility  arguments,  to  be  the 
Burgers  vector  of  the  dislocation  (-  interatomic  spacing), 
both  Basinski  et  al .  and  Ziman  obtained  reasonable  agree¬ 
ment  with  experiment . 

§5.2.2  The  Present  Work 

The  theories  discussed  above  are  inadequate  in 
one  or  both  of  the  two  important  aspects  of  a  resistivity 
calculation:  (1),  the  calculation  of  the  probability  of 

scattering  P(k,k* )  and  (2),  the  calculation  of  the  resis- 


t 
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tivity  using  this  probability  of  scattering.  Because 
of  their  inability  to  deal  properly  with  the  core,  the 
theories  of  Koehler,  Mackenzie  and  Sondheimer,  Landauer 
and  of  Hunter  and  Nabarro  are  mainly  inadequate  in  the 
first  aspect.  Their  calculation  of  the  resistivity  is 
based  on  approximate  methods,  the  reliability  of  which 
has  not  been  tested  before.  On  the  other  hand,  the 
theories  of  Harrison  (1958),  Basinski  et  al.  (1963)  and 
of  Ziman  (1964)  are  mainly  inadequate  in  the  second 
aspect.  Their  estimate  of  P(S,.S’),  though  still  not 
satisfactory,  is  considerably  improved,  but  their  cal¬ 
culation  of  the  resistivity  is  crude:  Even  the  aniso¬ 
tropy  of  the  dislocation  resistivity  has  not  been  taken 
into  account. 

The  present  work  was  started  in  an  effort  to 
construct  a  theory  satisfactory  in  both  the  aspects. 

The  calculation  of  P(k,k')  has  now  become  possible  (see 
Chapter  2)  by  the  development  of  the  pseudopotential 
technique.  In  this  method,  we  are  no  longer  bothered 
by  infinite  displacements  since  what  matters  is  the 
pair  correlation  function  or  the  relative  displacements. 
The  manner  in  which  the  strains  may  vary  is  also  not  a 
limitation  any  longer.  For  calculating  P(k,k'),  all 
we  have  to  know  is  the  form  factors  of  the  ions  and  the 


structure  factor. 
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However,  the  solution  of  the  second  difficulty 
is  not  quite  as  simple.  Difficulties  in  this  arise 
when  we  try  to  solve  the  Boltzmann  equation  in  the 
presence  of  such  a  highly  anisotropic  object  as  the 
dislocation.  One  is,  therefore,  faced  with  two 
alternatives  -  either  to  use  some  approximate  method 
for  solving  the  Boltzmann  equation  or  to  use  some 
simple  model  for  the  dislocation  so  that  the "Bolt zmann 
equation  may  be  solved  exactly:  We  decided  in  favour 
of  the  latter.  Treating  the  dislocation  to  be  a  row 
of  weak  scatterers  each  having  a  spherically  symmetric 
potential,  we  have  been  able  to  obtain  an  exact  analy¬ 
tical  solution  of  the  Boltzmann  equation  in  the 
presence  of  a  set  of  parallel  dislocations.  This 
decision  was  taken  for  several  reasons.  Firstly,  the 
preliminary  investigation  made  by  Harrison  shows  that, 
though  crude,  the  model  is  not  far  from  reality  and 
it  is  of  interest  to  obtain  a  reliable  estimate  of  the 
resistivity  on  its  basis.  Secondly,  if  we  can  extend 
our  solution  to  the  case  where  the  potential  due  to 
the  individual  scatterers  is  allowed  to  have  radial 
asymmetry  in  the  plane  perpendicular  to  the  row,  the 
model  would  simulate  the  scattering  due  to  an  exact 
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4* 

arrangement  of  the  ions  in  an  edge  dislocation  (see 
figure  1).  Thirdly,  an  exact ' solution ,  if *  only  for 
a  simple  model,  is  important  because  it  can" be  used 
as  a  guide  to  the  reliability  of ' the- various  appro¬ 
ximate  procedures;  or,  if  one  wants  to  use  the 
variational  technique,  for  making- a  more  educated 
guess  about  the  trial  function  rather  than  using  the 
simplest  trial  function  given  by  (1.20). 

To  take  care  of  the  expected  deviations  from 
Matthiessen ' s  rule  we  have  solved  the  Boltzmann 
equation  in  the  presence  of  other  mechanisms  of 
scattering,  but  in  order  that  our  method  applies  in 
its  present  form,  we  require  to  put  the  restriction 
that  the  scattering  due  to  these  other  mechanisms 
depends  only  on  the  relative  orientation  of  5  and  Icf  . 

Chapter  6  deals  with  the  solution  of  the 
Boltzmann  equation  and  Chapter  7  with  the  calculation 
of  the  resistivity.  The  resistivity  due  to  both 
parallel  and  random  dislocations  is  discussed.  For 
random  dislocations,  an  alternative  to  the  usual  method 
for  determining  the  resistivity  is  discussed. 

+ 

In  as  much  as  the  core  a  screw  dislocation  also  has 
a  dilatation  [ see  e.g.,  Nabarro  (1968)  p.  6±T]  ,  the  simple 
model  used  here  provides  an  estimate  for  the  resistivity 


due  to  a  screw  dislocation  also. 


■ 
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Fig.  1.  A  unit  edge  dislocation  in  a  simple  cubic 
lattice.  The  dislocation  is  along  the 
z-axis  and  the  slip  direction  along  the  . 


x-axis . 


Chapter  6 


THE  SOLUTION  OP  THE  BOLTZMANN  EQUATION 
IN  THE  PRESENCE  OP  A  SET  OF  PARALLEL  DISLOCATIONS 

In  this  chapter  we  obtain  a  solution  of  the 
Boltzmann  equation  for  the  distribution  function  in 
a  crystal  in  the  presence  of  a  set  of  parallel  dis¬ 
locations  treating  each  dislocation  to  be  a  row  of 
scatterers  having  a  spherically  symmetric  potential. 

Along  with  dislocations,  other  mechanisms  of  scattering 
may  also  be  present  with  the  restrict ion' that  the 
scattering  due  to  the  other  mechanisms  depends  only 
on  the  angle  between  £  and  £’ .  The  energy  surfaces 
in  k-space  will  be  assumed  to  spherical. 

One  approach  to  this  problem' is  to  try  to  solve 
the  integral  equation  (1.7)  for  the  anisotropic  relaxa¬ 
tion  time  tj.  For  the  simple  case  of  isotropic  scatterers 
(1.7)  was  solved  by  Bhatia  (19^9)  but  the  method  does 
not  appear  to  be  adaptable  to  an  easy  generalization. 

The  method  we  use  here  is  based  on  an  approach 
used  previously  by  Bhatia  and  Moore  (1961).  It  consists 
in  reducing  the  solution  of  the  Boltzmann  equation  to 
the  solution  of  a  difference  equation,  the  unknowns  of 
the  difference  equation  being  fLM(k),  the  coefficients 
expansion  of  the  distribution  function  fg  in  spherical 
harmonics.  In  general,  this  problem  may  be  as  difficult 


. 
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as  the  original  problem  itself.  But  in  our  case  it 
turns  out  that  the  solution  of  the  difference  equation 
can  be  reduced  to  the  solution  of  a  number  of  linear 
equations  for  the  unknowns  From  here 

f kM ’ s  can  be  obtained  if  certain  quantities  associated 
with  the  scattering  due  to  the  dislocations  and  the 
other  mechanisms  can  be  expanded  in  a  finite  number  of 
Legendre  polynomials.  If  £q  and  Lq  respectively  are 
the  degrees  of  the  highest-degree  polynomial  in  the 
expansion  of  the  relevant  quantities  in  the  dislocation 
scattering  and  the  scattering  due  to  the “other  mechan¬ 
isms,  Lq  is  required  to  be  finite  so  that  the  number  of 
equations  we  have  to  solve  is  finite  and  is  required 
to  be  finite  so  that  the  coefficients  of  the  unknowns 
in  these  equations  may  be  calculated.  Since  it  is  reason¬ 
able  to  assume  that  in  a  physical  problem  both  £q  and  Lq 
would  be  finite,  the  whole  problem  boils  down  to  solving 
a  finite  number  (usually  of  the  order  of  5)  of  linear 
equations . 

In  §6.1  we  shall  obtain  the  difference  equation 
for  the  f^’s;  the  solution  of  that  will  be  discussed 
in  §6.2. 

§6.1  Derivation  of  the  Difference  Equation  for  the  f Ljvi T  s 

The  Boltzmann  equation  for  the  distribution  function 
is  given  by  (1.1).  Substituting  for 


les  +- 


—  *P  n  o  1  H 


from  (1.2b) 


~  tP  "  ri  . 
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and 


9  t 


coll 


from  (1.5)  and  recalling  that  the  kf- 


integration  contains  a  factor  ft/8iTJ  which  had  been 
left  understood,  we  get 


e  ^  3  fk 
fi  £,k  9k 


ft 


8 


7T 


3 


/P(k,k' )(f£-fgt )dk ’ 


(6.1) 


We  shall  assume  that  the  probabilities  of  P^(k,k’)  and 
P^Ckjk')  of  scattering  due  to  dislocations  and  other 
mechanisms  respectively  are  additive.  Thus 


P(K,&»)  =  P1(J,$»)  +  P2(K,S»)  ,  (6.2) 

in  which,  as  mentioned  before,  P2(£,k’)  will  be  assumed 
to  depend  on  |k-k'|  only. 

Expanding  fj  in  terms  of  spherical  harmonics, 

we  have 


f£  fk  +  ELM  fLM  YLM  5 

where  ( k ,  3 , )  are  the  spherical  coordinates  of  k  with 
respect  to  the  system  of  coordinates  in  which  the 
polar  axis  is  along  the  direction  of  the  dislocation 
lines.  Here  the  sum  represents  the  deviation  of  f+ 
from  f°  on  the  application  of  the  electric  field.  The 
L-summation  is  over  all  positive  integers  and  M  takes 
all  integral  values  from  -L  to  +L.  The  prime  over  the 
summation  sign  stands  for  the  exclusion  of  the  term 
L=0  but  it  will  be  left  understood  in  future. 


Writing  a  similar  expansion  for  fj,  we  have 


f£'  =  fk  +  lM  fLM  YLM  (a’<t>)  ’ 


(6.4) 


where  [k f  ( =k)  ,a ,  (J)}  are  the  polar  coordinates  of  kT  . 
Note  that  the  scattering  being  elastic,  the  co¬ 
efficients  f^M  and  the  isotropic  part  f°  are  the  same 


for  f^  and  f^, . 
k  k 1 

Using  ( 6 . 2 )  -  ( 6 . 4  )  and  again  taking  the  energy  surfac 
to  be  spherical,  equation  (6.1)  may  be  written  as 


9  f 


o 


e  -*  r  '  k  _  t  ,  T 
~  fi  £-k  3k  Ii  +  I; 


(6.5) 


where 


h  ■  A  /k'^Ck.kUE  fLM[YLM(6^)-YLM(a)4>)]dk' 

OTT  LM 

sinadadcf)  ,  (6.6) 

and  Ip  is  given  by  a  similar  integral  but  with  P^ 
replaced  by  Pp. 

First  let  us  consider  1^.  We  have 
Pl(£,£’)  =  ^[M1(k,K1  )  ]26(Ek-  Ek,)  ,  (6.7) 

where  M1(k,£t)  is  the  scattering  matrix  element  due  to 
the  dislocations.  For  one  dislocation  M1(k,k’)  is 
given  by 
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so  that 


m1(SE,««)|2  =  ^  [w(q)] 


N 


where  n  is  the  number  of  scatterers  constituting  the 
row,  a  the  distance  between  the  scatterers,  w(q) 
their  form  factor  and  N  the  total  number  of  ions  in 
the  crystal.  (Note  that  since  the  potential  due  to 
individual  scatterers  has  been  assumed  to  be  spheri¬ 
cally  symmetric,  w  depends  only  on  the  magnitude  of  q.) 

If  there  are  Nq  dislocations  in  the  crystal, 
P^(£,k’)  will  be  given  by 


(6.8) 


Since  q  depends  upon  k  and  0  (q=2k  sin0/2),  where  0  is 
the  angle  between  S  and  5’,  we  can  write 


[w(q)]2  =  E  c  p 1  ^  ( k)  P  (cos0) 
£=  0  *  * 


(6.9) 


where  P^(cos0)  are  the  Legendre  polynomials  .  Substitu¬ 
ting  (6.8)  and  (6.9)  in  (6.6),  noting  that 


No  confusion  should  arise  from  the  use  of  symbol  P 


for  probabilities  as  well  as  for  Legendre  polynomials. 
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/ k’2  6(Ek-Ek,)dk'  =  k2(|§)  1  , 


(6.10) 


and  using 


P^(cosQ)  = 


2ITT  l  * 

m 


(6.11) 


we  get 

n  n  k2  c  p 1 ^ 

Z1  =  IffiN  ^dk^  E  2£+I  fLM  Y£m^*^  x 

£mLM 


2 

*  sin  [Jgnka  ( u)  ’ -go  )] 

x  - p - 

sin^  [h  ka(o)f  -w)] 


[YLM^^^"YT.M^a^^]da3’d(i)  , 


LM 


(6.12) 


where  a)  =  cos$  and  w1  =  cosa.  Here,  the  integration 
over  a) '  behaves  as  a  delta  function  and  we  get 


11  £LM  °'Z  fLM  {YlM  (e,4,)^o(u)J 


Y£M ^  5  ^  ^£M ^ w  ^  w  ^ 


=  E  cif^LM(3^){^Pn^)]2-PP.M^)J2>,  (6.13) 


£LM 


£o 


£M 


where  /  n 

N  nfi  kc 

i  o  o  £  /dEx-1 

C  £  dk 

1  (2£+l)Nfta 


(6.14) 


and  fry,  are  defined  by 


YLM(e’^  = 


/  4  7T 


^LM  (u)  e1”'1' 


(6.15) 


/ 
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or,  equivalently,  by 


.  M 

^LM  =  NLM  PL^  5 


(6.16) 


M,  N 

where  P^(oj)  are  the  associated  Legendre  polynomials  and 
Nlm  are  given  by 


N 


LM 


(2L  +,  1) 


(L-M) I 
(L+M) ! 


(6.17) 


o 


Noting  that  the  density  of  dislocations  is  equal  t 
(NQna/fl),  we  can  write  c^  in  terms  of  N^.  This  gives 
Nnfi2kci1)  , 

„  t  _  Do  Z  /  dE  ^  —  1  f  c  q  \ 

c* '  771^7  <®)  '  <6'181 


Now  let  us  evaluate  I^.  For  P^Ckyk’)  we  can  write 
down  an  expression  similar  to  the  expression' ( 6 . 7 )  for 


P^(k,Kf).  But  now  the  matrix' element' would' depend  only 
on  the  relative  orientations  of  k  and  k*.  Therefore  we 
can  write 


00 


— /  k f  2P9  (  |  Ic-S  ’  |  )dk f  =  k2  Z  c(2)(k)Po(cos0)  . 
8ttJ  ^  £=0  36  36 


(6.19) 


Substituting  (6.19)  in  (6.6)  and  using  (6.11),  we  get 


J2  =  E  c£fLMY£m^J^)/Y£m(a,(J))  [YTM(0,i|;)-YTM(a,(l))Jdmfd4) 
£mLM 


Em  (co  °L)  fLM  YLM  > 


(6.20) 


where 

"  _  (2)  4Trk2 

C  £  C  £  2  £+1 


(6.21) 


. 
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Next  we  consider  the  left  hand  side  of  the 
Boltzmann  equation  (6.5).  If  we  choose  the  z-axis 
along  the  dislocation  lines  (which  is  also  the 
direction  of  the  polar  axis  for  the  spherical  co- 
ordinates  we  are  using),  e.k  is  given  by 


e.k  =  sing  cosip  +  e^  sin3  sinip  +  e^  cos3* 


(6.22) 


Hence 


.o 


J-  £.k  =  Z  A,  Y-.  (3,iJ/)  , 

fi  9k  m  lm  lm  * Y  9 


(6.23) 


where 


e  9f°  /47Txi 


A 


10 


(^) 


A 


1,±1 


fi  9k  v3 


e  3  f  f  2  tt  r  -  a  \ 
fi  3  k  3  x  y 


(6.24) 


Using  (6.23),  (6.13)  and  (6.20),  the  Boltzmann 
equation  (6.5)  reduces  to 

(co  °L)  fLM  YLM 


+  ^  °£  YLM  YLM^ ^  j  ^ 


&LM 


ilo 


M 


E  A-,  Y-.  (3,40 

m  lm  lm  5  Y 


(6.25) 
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The  Difference  Equation 


We  have  thus  expressed  the  Boltzmann  equation 


(6.5)  in  a  form  where  instead  of  determining  f 


5 


directly,  we  can  determine  its  components  f 


LM* 


The 


fact  that  T  s  are  infinite  in  number  does  not 
necessarily  make  this  reformulation  more  difficult 
because  firstly,  many  of  them  would  be  identically 
zero;  and  secondly,  for  determining  a  physical  property 
of  interest,  we  may  not  require  all  of  the  nonzero 

j 

fl^’s,  e.g.,  as  will  be  seen  in  the  next  chapter,  for 
determining  the  three  components'  of  J,  hence  for 
determining  the  conductivity  (or  resistivity)  tensor, 
we  need  only  three  of  them,  namely' f^Q  and  f-^ 

Now  (6.25)  contains  all  the  f’s  mixed  together. 
The  f’s  we  want  may  be  disentangled  from' the  rest 
partly  by  using  the ■ orthogonality • of ' the  spherical 
harmonics,  which  gives"  a' difference' equation  for  f^, 
and  partly  by  solving  the  difference  equation  thus 
obtained . 

To  derive  the  difference  equation  for  f^  from 
(6,25)  we  proceed  as  follows:- 


1. 


In  order  to  be  able  to  take  advantage  of  the 


orthogonality  of  the  spherical  harmonics,  we  try  to 
rewrite  (6.25)  in  such  a  manner  that,  just  as  on  the 


. 
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RHS,  each  term  on  the  LHS  also  is  a  single  spherical 
harmonic.  This  may  be  done* in  two  steps:  First 
express  the  terms  involving  products- of  B )  and 

other  angle  dependent  quantities'  as ~ a" series  of  single 
spherical  harmonics;  then,  rearrange  the  terms  to 
get  the  desired  form. 

Thus  we  write 


(  Cn  °L  ^  YLM^  ^  +  ^  0  C  P.^T.M^  ^  W  ^  1?  p  M  ^  W  ^  ^ 


£  £  LM 


£o 


£M 


=  ZjC(j,L,M)  YL+jjM  (&,if>)  , 


(6.26) 


where  j  takes  positive  as  well  as  negative  values  and 
c(j,L,M)  can  be  determined* in ‘ terms' of  the  coefficients 
cf  and  c" .  Combining  (6.26)  and  (6.25)  gives 

j LM  >LJM)fLMYL+j  =  ^mAlmYlm^  ^  ^  9 

i 

which  may  be  rewritten  as 


£  o(jsL-j,M)f  Ylm(6^)  =  ImAlmYlm(e^)  . 

j  LM  u  5 


(6.27) 


which  is  the  desired" form . 

2.  Using  the  orthogonality  of  spherical  harmonics, 

now  gives 


Zc(j,L  j,M)  A1M  6l1 


Mm 


(6.28) 


so  that  the  RHS  of  (6.28)  is  zero  unless  L=1  and 
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M=m(=0,  ±1).  If  the  number  of  terms  in  the  above  sum 
is  finite,  (6.28)  is  the  required  difference  equation. 

Note  that,  in  general,  the  RHS  of  (6.26)  will 
be  of  the  form  c(  j  ,L,M,Mf  )YL+J-jM+M»  (  3,^) ,  in  which 

case  the  LHS  of  (6.28)  would  be  Z  , c ( j ,L-j ,M-M ’ ,M ’ )  x 

f^  .  jyj _ jyj |  .  Here  a  simpler  form  for(6.26)  results  from 

the  fact  that  Y^(3,^)  is  multiplied  by  terms  indepen¬ 
dent  of  ifj .  Therefore,  while  there  is  a  mixing  of  L  in 
the  difference  equation,  there  is  no  mixing  of  M.  In 
a  more  realistic  model  of  the  dislocation,  that  may  also 
be  expected. 

Now  we  proceed  to  evaluate  c(j,L,M)  for  the  three 
values  of  M(=0,  ±1)  for  which  the  RHS  of  (6.28)  is 
different  from  zero. 

For  M=0,  we  have  from  (6.26)  that 

Zj  c(j,L,0)YL+j  y3,iM  =  (c"  -  c^)Ylm(3,40  . 


c( j 5L,0)  =  (c”  -  c”)  6jQ  . 


Combining  this  with  (6.28)  gives 

r  -  Al°  *  _  A10  * 

L0  "  M  LI  "  M  LI 
c  —  c  T  c  —  c 


(6.29) 


As  regards  M=±l,  it  is  clear  from  (6.26)  that 
c(j,L,M)  will  be  the  same  for  both  of  these  values. 
So  let  us  rewrite  (6.26)  for  M=l.  This  gives 


Ill 


(co  ■  CL)  +  11  0iYLl(e’,('){Pjlo(to)n2_P)ll(“)]2} 


i 


Z  c(j,L,l)  Yl  (B,i|0  . 
J 


(6.30) 


To  determine  c(j,L,l),  we  note  two  points:  one, 

E  c^{p^0(w)J  -  0^£]_(w)j  } when  written  in  powers  of  w 

will  contain  only  even  powers  of  oj;  two,  co  n  (n  is 
integer)  when  multiplied  by  P^(w)  gives  a  series  of 
terms  P^+2r  w^ere  r  varies  from  -n  to  +n.  Thus  we  may 
write 


j 


l  ch^to(“)J2-l?)!.i(“)^2}=  ao+  al“2  + 


+ 


2£ 


o 


(6.31) 


and 

w2n  P2(w)  =  nk°(L)PY(o))  +  nk[2)(L)PL+2(“) 

+  nk22)(L)PL-2(“)  +  nki4)(L)Pi+4(“>  + 


nk^)(L)PP_4(o)) 


+  . +  nkl  ^  ('L-)PL+2n('“->+nk22n'l('L')PL-2n('“^  - 

( 6 . 3  2 ) + 

£q  in  (6.31)  is  the  highest  Legendre  polynomial  in  the 

2 

expansion  of[w(q)J  -  see  (6.9).  In  (6.32),  the  super¬ 
script  on  k  denotes  the  magnitude  of  the  change  in  L,  the 


t 


If  2n  >  L,  clearly 


k0 
n  2 


( 2r ) 


(L)  will  be  nonzero  only  for 


2r  <  L . 


» 
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pre-subscript  is  self-explanatory  and  the  post-subscript 
indicates  whether  L  has  decreased' or  increased  -  1  denotes 
increase,  2  denotes  decrease,  and  no  post-subscript  means 
no  change  in  L.  It  may  be  mentioned  that  the  coefficients 
a  in  (6.31)  and  the  coefficients~k  in  (6.32)  can  be 
determined,  but,  as  will  be  seen  later,  we  will  not  need 
to  determine  them. 

Substituting  (6.32)  and  (6. 31)  in  (6. 30),  we 
immediately  obtain 


c 2  2  ^  5 


L,  1)  =  [a 


(L) 

r  r  1 , 2 


a  k(2r)(L) 

ar+l  r+1  ± , 2 


o  Nt  .. 

,  2r  ,T  v-r  LI 

+  ....+  a  k  2(L)J  n 

o  o  5  L+2r,l 


+  (c"-c")6  ,  (6.33) 

o  L  ro 


where,  just  like  k^  and  k2 ,  and  c2  correspond  respec¬ 
tively  to  increase  and  decrease  in  L. 

Combining  this  with  (6.29),  we  obtain  the  difference 
equation 

-  M  M  . 

(C0  "  CL  fL  + 

£ao+  a1  k°(L)  +  a2  2k°(L)  +  ....+  a£  £  k°(L)]fL 

o  o 


(2) 


(2) 


+  [a-^  ^k|^'/(L-2)  +  a2  2k^'(L-2)  +  ... 


+  a*  l  k22) 
o  o 


X  (L-2)]f\_ 


L-2 


+  [a-,  -1k^2^(L+2)  +  a0  0k^2^(L+2)  +  ....  +  a  „  „  k  ^  2  ^ 


'1  12 


£  £  2 
o  o 


x  (L+2)]f 


L+2 
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+  [a2  2k^(|)(L-4)  +  -  +  a  k^4;(L-i))]fL_4 

*0  o 

+  [a.  2  2^2  ^  (^+4 )  +  ••••  +  a^  ^  ^  ( L+4  )J  f  +....+ 

o 

+  a*  l  ^2to>(L-240)f2_2t  +  a£  k<2V(L+2i  )f 

oo  o  o  o  o 


(4) 


=  A  6 


LI 


(6.34) 


where  f^  5  f^  NL1  and  A  =  A^  . 

The  difference  equation  for  M=-l  will  be  the  same 

except  that  on  the  RHS,  we  will  have  A,  ,  Instead  of 

-L  ) 

An* 


§6.2  Solution  of  the  Difference  Equation 


The  method  of  solving  (6.34)  is  suggested  by  the 
relation 

a(L-2r)  mk^2r)(L-2r)  =  a(L)  mk^2r)(L),  (6.35) 

where  a(L)  =  ( 2L+1  )/t2L(L+l  )]  5  which  is  proved  in  Appendix 
II. 


First  we  rearrange  (6.34)  to  obtain 


(cM  -  c j  +  a  )  fT 
o  L  o  L 


+  ax  [1k0(L)fL  +  1k[2)(L-2)fL_2  +  ’  (L+2  )fL+2] 


(2) 


+ 


+  a«.  U  k°(L)fL  +  1  k["4L-2)fL_2  +  %  k<2)(L+2)f 
o  o  o  o 
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+  ••••  +  £  ki2ib)(L-24  )fL_  +  ,  k<2V(L+2*o)f 

o  o  o  o 


3 


=  A  6 


LI 


(6.36) 


It  Is  instructive  next  to  consider  the  special  case 
where 


"  "  o 

Ct  C  6-rp, 

L  o  LO 


(6.37) 


With  this  (6.36)  becomes 


(co  +  ao^fL  +  al^lk  ^L^fL  +  lkl  ^  2^fL-2  +  lk2  ^L+2^fL+2^ 


+  - +  a0  [.  k° ( L ) f T  +  0  k^;:(L-2)fT  Q+  k^;(L+2)f 


(2) 


(2) 


'£  L£ 
o  o 


L  £  1 

o 


L-2  £  2 

o 


L+2 


+  t  ^2V(L-2*o)f  +  ,  k(2V(L+2*  )f 

o  00  o 


] 


=  A  6 


LI 


(6.38) 


The  solution 


^  _  2L+1  A 

L  2L(L+1)  " 


+  1 
/ 


2  4  2  £ 

cQ+ao  ^  1+y^o)  +y ^oo  +...+y^  w  o 

o 


(6.39) 


is  suggested  from  the  relation  (6.35).  It  may  easily  be 
verified  by  substitution.  This  is  done  in  Appendix  II  , 
where  we  also  prove  (6.35).  Yr  in  (6.39)  stand  for 

m 

a  /  ( c  +  a  )  . 
r  o  o 

Now  we  can  readily  see  the  modification  that  would 


be 


necessary  if  a  larger  number  of  c^  are  nonzero,  but  if 


. 
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for  L  greater  than,  say  L  ,  are  zero.  For  this 
purpose,  we  rewrite  (6.36)  such  that  the  LHS  is  the 
same  as  in  (6.38),  so  that  RHS  would  be 


(A  +  c1  f 1 ) 6L1  +  c  f  6L3  +  -  +  cL  f  6ll  (6.40) 

o  o  o  5 


where  for  obvious  reasons,  we  have  taken  L  to  be  odd 

5  o 

IT 

Thus  c-^  is  the  highest  nonzero  coefficient, 
o 

Comparison  of  (6.40)  with  the  RHS  of  (6.38) 


suggests  the  solution 


fL  =  (B  +  W  2LIE7IT  ' 


P^(o))P^(o))  dm 

2  2l~ 

1+Y-,  a)  +  .  .  .  .  +y  Q  a)  o 

36  o 


,  o  2L+1 
P3  2L( L+l )  J 


P^(u))P3(aj)  da) 

2 

1+Y-i  w  +....  +Y  0  w  0 
1  *0 


4.  ,  o  2L+1  , 

* ‘ •  PL  2L(L+1)  J 
o 


pRuOPl  (u>)  du) 

i-j 


2  2  £  5 

1+Y-i  w  +  .  .  .  .  +Y  n  w  o 

o 


(6.41) 


where 


p  -  A 

b  -  — n - 

c  +a 
o  o 


and  3r  =  — 


c  +a 
o  o 


The  verification  of  ( 6 . 4 1 )  follows  closely  that  of  (6.39) 
and  will  not  be  shown  explicitly. 

L  +1 

Using  (6.4l)  for  L=l,3*  ....,  Lq  we  get  — | — 

equations  in  as  many  f’s  from  where  f, ,f_,  ...,  fT  may 

d  Lo 

be  obtained  (and  once  these  f's  are  known,  f^  for  any 


.c  It  s  91*  t  a  « fieri:*  is:*£9‘i3  J  rrc  ,o 

. 


J  crs^ficf  sveri  ew  ,  anos  t©i  euoivdo  *10!  3  isriA' 


0.  ,  y+  *  •  •  • 


(  i  ,  i)  )  B  t  x,  >eoIa  8 wo.  o;  (I  .6)  ictde*  !  *t<  v  t:;,. 
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higher  L  may  be  evaluated  from  ( 6 . 4 1 )  directly.  (For 
the  purposes  of  electrical  resistivity,  of  course,  we 
need  only  f  . ) 

The  system  of  linear  equations  we  obtain  is 


(C1  "  l)fp  +  c3  1-^2  + 


+  CL  J1L  fL  +  A  Z11  0 

o  o  o 


C1  I31  fl  +  (c3  I33  1)f3+ 


It 


+  cl  IOT  fT  +  A  I0,  =0 
Lo  3Lo  Lo  31 


II  -r  „  .  M 

C1  ZLol  fl  c3  ZL  3  f3+ 


+  ^CL  IL  L  ~1’)fL  +  A  ZL  1  0  ’ 

o  o  o  o  o 


where 


'LL 


2L+1 


+  1 
/ 


PZ(m)  PZ  f (m)  dm 


2L(L+1)  \  "  .  2  ,  2 £ 

-1  c  +a  +an  a)  +....+a„  w  o 

001  £o 


(6.42) 


2L+1  +Z 
2L ( L+l )  { 


PZ(m)  PZf(m)dm 


-1  c"  +  E.c{ ■  {[*>,.„  (co)]2-p,,(“)]2} 


(6.43) 


(Note  that  1^,  and  1^,^  are  not  equal  but  are  related 

.  n  ^  2L+1  L '  ( L  ’  +1 )  \ 

by  a  simple  constant  ^  ) —  )  • 

For  f11  this  gives  (recall  f11=  f1/N11J  An=  A/N-^) 
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A11 

f  =  — — 
11  A 


'll 


(c^I-^^^l)  •  •  *o  1^-p 
J  o  J  o 


o  '*  T 

3  53 


•  • -CL  I5L 
0  0 


C3IL  3  ‘  *  *  (  CT.  ZT.  T ) 


'L  L  L 
o  o  o 


-I 


31 


c3Il3---°LoIlLo 


"  "  T 

353 ’ - •cLo15Lo 


C3IL„3‘ • • (eL  XL  L 
o  o  o  o 


+  (-lfLo  ±S)/2  I 


L  1 
o 


IT 

c3 

hs  •••• 

TT 

•••  CL 

1 — 1 

H 

0 

0 

tt 

TT 

°3 

T  -1,3 

• 

• 

• 

O 

tr1 

lL  , 

0  5 

0 

0  5 

(6.44) 


where 


/  M  \ 
C1I11~  Z' 

It 

C  ->  i  .  .  . 

....  cT  i 

L  1L 

0  0 

o"T 

131 

( C3I33— 1 ) 

TT 

. . . .  cL  I3L 

0  0 

TT 

C1ZL  1 

0 

TT 

c3ZLo3  “ 

1 — 1 

1 

0 

0 

H 

O 

=  ^ 

O 

• 

It  may  be  noted  that  f^  (and  so  also  f^i5  ^51  . )  Is 

proportional  to  A-q  so  that  the  deviation  of  f  from  f°  is 
just  proportional  to  certain  combination  of  components  of 
the  electric  field,  as  it  must  be. 

The  value  of  f^  1  (or  of  f^  f^  ^  etc.)  may  be 

obtained  from  f^1  (or  from  f^,  f^  etc.)  by  replacing  A^1 

by  Ai,-r 


X 


Chapter  7 


CALCULATION  OF  THE  RESISTIVITY  DUE  TO 

DISLOCATIONS 

In  the  last  chapter  we  solved  the  Boltzmann 

equation  for  a  crystal  containing  parallel  dislocations 

and  some  other  mechanisms  of  scattering.  Here  we  first 

use  these  results  to  calculate  the  resistivity  tensor 

p.  .  for  a  metal.  Then  we  discuss  the  resistivity  due 
^  J 

to  random  dislocations.  Finally,  we  present  some 
numerical  estimates  of  the • resist ivity  obtained  from 
the  model  of  the  dislocations  we  have  used. 

§7*1  Resistivity  Tensor  for  Parallel 

Dislocations 

As  discussed  in  Chapter  1,  the  evaluation  of  . 

is  based  on  determining  J.  in  terms  of  e . .  J  is  given 

*  J 

by  (1.9).  Assuming,  as  before,  the  Fermi  surface  to  be 
spherical,  we  get 

J  =  — f  k  k2dk  dft 

4tt  3-fx  k  9  k 

=  — I-  /  fr*  k  k2  dE  dft  . 

,,  3^  k 

4tt  ti 

Therefore  J  would  be 
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Substituting  for  f^  from  (6.4),  we  get 


Jz  = 


4'jT:5fi 


57  r  [f|  +  'z  fLMYLM(6,4-)]  [(^^^(e^^k^ndE 
'*S|  LM 


e  /  4tt  \h 


3*  3 


(^-)  2  /  f-,  n  k  dft  dE  . 


10 


(7.1) 


Similarly,  for  Jx  and  we  get 


J  =  — I —  (l2-)1^  /  (f, ,+  f,  .)k2dfldE,  and 
x  3  11  i’_i 


(7.2) 


Jv  =  —I —  (kL)k  /  (f,  ,  -  f, ,)  k2dft  dE  . 

y  4TT3fii  3  u-i  ii 


(7.3) 


For  the  model  of  the  dislocation  we  have  used, 
fios  f11  and  f1  are  proportional  to  A1Q,  A11  and 
Ai  ^  respectively  because  there  is  no  mixing  in  M. 
Therefore  we  can  write 


f  =  G  ( c f  c’M  A 
1M  1M  5  c  '  h1M 


(7.4) 


where  G1^(c’,cM)  stands  for  G1M(c|,  c^,  .  .  .  .  ,c^  ; 
c”,  ...,  c”  ).  G10  involves  only  c^  and  c^,  and  may 

be  immediately  obtained  from  (6.29);  however,  G^(=Gi  _^) 
in  general  depends  on  all  of  these  coefficients  and  has  to 
be  evaluated  from  (6.44)  numerically  [in  some  particular 
cases  (see  below)  the  dependence  turns  out  to  be  simple]. 

Substituting  (7.4)  in  (7.1  to  3),  we  get 


X.  !  as  1  i ol 

d  5ljj  9Vfld  9W  ftoldfiools  b  >r!3  To- <£  ^om  iritf  '10'? 

scf  Inv  n  5  3^  s^oleiedT 
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J  = 


J  = 
x 


J  = 


y 


2 

6  e  z  2 

2-2  f  G10k  ^ 
3tt  ti 

2 

e  ex  2 

— Hf  !  Gllk  ( 
3ir  ft  11 

2 

6  £  v  2 

,  2,2  S  Gllk  ( 
3  7T  Tl 


3  f  0  \ 

3  E  ' 

dE 

dk 

3  f°v 

3  E  ; 

dE 

dk 

3  f°x 

3  E  ' 

dE 

dk 

dE 


dE 


dE 


(7.5) 


Using  (1.11),  we  get 


2 
e  e 


J 


3tt  n 


z  _  /  t  tt\.2  /9  E  \ 

2.2  G10^cF>  cF5  kF  ^T^F 


(7.6) 


where  =  c(E„).  The  expressions  for  J  and  J  are 

similar  except  that  instead  of  G^q,  we  have  and 

instead  of  t  we  have  e  for  J  and  e  for  J  . 

z  xx  y  y 

In  all  the  three  cases,  depends  only  on  , 
which  shows  that  the  axes  we  have  used  are  principal 
axes.  The  three  principal  resistivities  are 


„  2.2 
3  7T  fl 


zz 


2  2 

e  kpG10 


/  t  tt  X 

( cp J  cp) 


(dE)-! 

vdk;F 


(7.7) 


P 


xx 


,  2.2 

3  TT  fl 


2,  2n  /  t 

e  kFGll(cF 


(dE)-l 

vdk;F 


(7.8) 


Note  that  since  G^q= 


the  value  of  p  is 

z  z 


the  same  as  in  the  absence  of  the  dislocations,  which 
is  to  be  expected  since  a  dislocation  does  not  affect 
the  periodicity  of  the  lattice  parallel  to  itself. 


'  • 
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Furthermore,  for  the  model  of  the  dislocation  we  are 
using,  when  referring  to  the  principal  axes,  there 
are  only  two  independent  components  of  p. .. 


For  ApD ,  (7.7)  and  (7.8)  give 


( APD ) 3  =  °> 


( ApD)i,2 


2  2  -1 
3tt  n  /  dEh  1 

2,  2  ldk;F  ~  9 
e 


D 


(7.9) 


where  D  =  D-/Gi:l(  c  '  ,c" )  -  (e"  -  c")}  .  (7.10) 

Here  the  subscript  F  is  left  understood  in  the  c's.  It 

2 

should  be  noticed  that  since  c’  «  k,  cM  «  k  ,  the 
resistivity  due  to  the  other  mechanisms  (which  depends 
only  on  cM)  depends  upon  energy  in  a  different  manner 
than  the  resistivity  of  dislocations  (which  involves  both 
c*  and  cM)o  Therefore,  the  effect  of  dislocations  on 
thermoelectric  power  which  is  proportional  to  d  in  p/dE, 
could  be  very  different  from  that  on  the  resistivity. 

In  general,  this  is  about  as  far  as  one  can  go 
analytically.  Now  we  discuss  a  few  particular  cases  of 
interest . 

1.  c!  =  6„  (i.e.  the  scatterers  constituting  the 

X/  Xj 

dislocation  scatter  isotropically):  In  this  case,  D 

t  t 

becomes  a  which  is  equal  to  c  when  other  cn  are  zero. 
This  is  the  same  result  as  obtained  by  Bhatia  (19^9) 
with  the  difference  that  to  show  this  result,  we  have 
not  made  any  assumption  regarding  the  coefficients  c" 

M  M 

while  Bhatia  assumed  cT  =  c  orn. 


rf}  ,  'A  »  "o  tA  10 
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I!  | 

2.  (co+aQ)  >>  for  £=1,2,... ,J lo:  In  this  case 

we  get 


D  -  a0+  |  f  rPiU)]2(aia)2+a2A - +a^ow2£o)da),  (7.11) 

where  the  correction  terms  are  of  order 


£ 

-tr+~a'  ~  =  1  s  2, . ,  £o)  and  depend  upon,  apart  from 

o  o 

c^  ,  on  c ”  and  c”  also.  These  correction  terms  repre¬ 
sent  the  deviations  from  Matthiessen ’ s  rule.  (In  the 
case  1  above,  there  are  no  deviations.)  Note  that  the 
deviations  referred  to  here  are  of  the  type  Ag  (see 
Chapter  4),  since  in  this  theory  we  have  taken  to  be 
zero  (see  Chapter  6). 

Using  (6.31),  we  get  from  here 

D  =  |  2  c l  f  [Pi(oJ)]2{[PAo(>**)]2  -  [Pj.hw)]2}  •  (7-12) 

£ 

This  can  be  easily  evaluated  and  it  gives 


D  =  -3 


E 

£ 


c^/(2£-l)(2£+3) 


32 

£ 


( 2  £-1 ) ( 2  £+3 ) 


( 7 . 13 )  + 


•4* 

The  expression  on  the  RHS  of  (7.13)  is  nonnegative 
it  must  be.  This  can  be  seen  by  using  the  definition 


c^  (eqn.  6.18)  and  by  noting  that 


£ 


P^tcos  ) 
(2  £  —  1 ) ( 2  £+3 ) 


1 


since  2  (2£_i)  (2£+3)  3  and  LP£(COS0)]<  1 

X/ 


1 


as 

of 
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where  prime  over  the  summation  sign  denotes  that  &=o 
is  to  be  excluded. 

This  approximation  essentially  corresponds  to 
that  used  by  Mackenzie  and  Sondheimer  to  calculate  the 
resistivity  for  their  model  of  the  dislocation  (see 
Chapter  5).  Also,  as  may  be  readily  verified,  one 
would  obtain  for  p  the  expression  (7.9)  with  D  given 
by  (7*13)  if  one  uses  the  variational  expression 
for  the  resistivity  with  the  simplest  trial  function. 

Numerical  values  of  resistivities  with  and 
without  using  the  approximation  (2)  above  are  presented 
in  §7.3. 

§7*2  Resistivity  of  Randomly  Oriented 

Dislocations 

The  problem  of  the  resistivity  of  random  disloca¬ 
tions  is  not  as  simple  as  it  may  seem  at  first  sight. 

There  appear  to  be  two  possible  ways  of  proceeding.  One 
is  to  take  the  ensemble  average* of  the" resistivity  such 
that  each  member  of  the  ensemble  contains  parallel 
dislocations  of  density  equal  to  the  density  of  random 
dislocations  in  the  crystal  the  resistivity  of  which  we 
wish  to  determine,  but  the  direction  of  the  dislocations 
in  the  various  members  of  the  ensemble  being  random.  The 
other  method  is  to  determine  fj,  hence  the  current  density 


. 
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and  the  resistivity,  in  the  given  crystal  by  taking 
the  dislocations  to  be  random  at  the  outset. 

All  the  theories  we  are  aware  of,  use  the 
first  method  -  without  even  mentioning  the  second. 

It  appears  to  us  that  both  the  methods ■ are  correct 
but  they  apply  under  different  experimental  situations. 
The  most  appropriate  criterion  for  determining  which 
of  the  methods  would  apply  in  a  given  situation,  seems 
to  be  the  relative  value  of  the  mean  free  path  A  of 
the  electrons  and  the  mean  distance  d(=N^  2)  between 
the  dislocations.  If  A  <<  d,  as  is  the  case  at  room 
temperature  when  the  dislocation  density  is  low  (N^ 

~  10°  cm  ;,  the  electrons  see  the  dislocations  of 
various  orientations  individually  so  that  one  should 
use  the  first  method.  However,  if  A >>  d,  as  is  the 
case  at  low  temperatures  for  a  relatively  pure  sample, 
the  second  method  would  apply.  Indeed  there  are  situa¬ 
tions  (A  ~  d)  intermediate  between  these  two  extremes 
when  neither  of  the  two  methods  would  be  strictly 
applicable.  For  the  intermediate  case,  it  is  not  clear 
what  one  should  do.  The  two  limiting  cases  will  be 
discussed  below. 

First  Method  (A  < <  d ) 

The  resistivity  in  the  direction  of  the  flow  of 
the  current  for  each  member  of  the  ensemble  is  J.e/J2. 
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Therefore  the  resistivity  Ap^  due  to  disordered  disloca¬ 
tions  is  given  by 


Ap 


o 

D 


i 


where  <  >  denotes  ensemble  average. 


J 

y 


is  given  by 


~  cos  y  +  — sin  y 

z  y 


where  y  is  the  angle  between  J  and  the  z-axis  (or  the 
dislocation  line).  This  gives 


A  o  _  i 
ApD  3 


( Ap 


D 


)3  +  3  (ApD)1 


5 


which  is  the  familiar  result  [see,e;g.,  Wilson  (1953) 
p.  1971  for  the  mean  resistivity  when  there  is  an  axis 
of  symmetry  for  the ' resist ivity .  (In  this  case  the  axis 
of  symmetry  is  the  direction  of  the  dislocation  for  each 
member  of  the  ensemble.)  Using  (7*9),  this  gives 


Second 


2TT2«2,dE>-l 
2,  2  ldkJF 

e  kp 

Method  (A  >> 


D  . 

d) 


(7.1*0 


Leu  P  (k,k')  be  the  probability  of  scattering  due 

to  random  dislocations  whose  total  number  in  the  crystal 

is  N  .  From  (6.8)  we  see  that  the  contribution  to 
o 

P°(£,£' )  from  dislocations  lying  between  6  and  6+d6,  and 


* 
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X  and  x+d-X,  where  6  and  X  are  the  polar  and  the  azimu¬ 
thal  angles  with  q  as  the  polar  axis,  we  get  from  (6.8) 


dP°(£,K ' )  = 


2MN  0  .2/3. 

or  (  \~i  2 r  / -r-,  „  \  sm  (Jgnqa  cos6) 

— 2C«(q)]  «(Ek-Ek'}  .  2,1 - — 

*RN  sin  (^qa  cos6) 


x  sin6d6dx 


(7.15) 


or  noting  that  k(cosa  -  cos$)  =  q  cos6.  This  gives 


2TTN 

p°(Jc,4«)  =  — | 

firr 


mr 

aq 


[w(q)f  5(Ek-Ek')  ' 


(7.16) 


The  result  for  the  resistivity  is  now  obtained  either 
from  (1.16)  or  directly  by  comparison  with  the  resis¬ 
tivity  given  by  P0(j£,$f)  which,  like  P°(k,kf),  depends 


only  on  the  magnitude  of  q.  We  get 

2  2 


3Tr2(m*)2fi  N  n  2 

O  OOrr/x-|  d.d  . 

f  \y{x)_\  x  dx 


PD 


,,  2^2,  3 
4e  fi  kp 


Na 


0 


3TT2(m*)2fi^  N  2 

- 3  -  -  ? -  /  [w(x)]2  x  dx  , 

4e  npk;;a^  0 


"F 


(7.17) 


1  . 


where  x  -  q/k  .  Note  that  because  of  the  factor  —  in 

F  q 


-r  -r  O 

the  expression  (7*16)  for  Pu(k,k’),  [w(x)J  in  the 
expression  for  is  weighted  with  x  rather  than  the 


usual  weight  factor  x  . 


o 


In  order  to  be  able  to  compare  it  with  given 


by  (7.1*0,  we  rewrite  this  as 


■ 


where 


N-p^  k-p  , _i  2 

\  =  §)p  /  w(x)  P  (!-§-)  dx  .  (7.19) 

2n  a 

Two  important  differences  between  p°  obtained  by 
the  two  methods  we  have  used,  may  be  noted 

1)  Since  in  the  first  method,  is  obtained  by 
averaging  the  corresponding  result  for  parallel  dislo¬ 
cations,  there  are  deviations  from  Matthiessen’s  rule 

for  random  dislocations  just  as  for  parallel  dislocations. 
However,  in  the  second  method  p^  is  obtained  from  P°(k,£’ ) 
which,  like  depends' only  on  |q|.  Therefore,  in 

this  there  are  no  deviations  from  Matthiessen ’ s  rule  (other 
assumptions  being  the  same) . 

2)  If  we  rewrite  the  expression  (6.18)  for  c^  as 
[using  (6.9)J 

^  (i^PTwU)]2  Pt(l  -  f)  d*  ,  (7.20) 

we  note  that  c^  and  A^,  occurring  in  the  first  and  the 
second  methods  respectively,  are  of  the  same  form  but 
are  weighted  differently  with  x.  Moreover,  while  the 
first  method  involves  all  the  (nonzero)  c^,  the  second 
method  involves  only  the  first  two  A0 . 


. 


b  rts 


128 


§7.3  Numerical  Results  and  their  Discussion 

The  object  of  the  numerical  calculations  was 
to  make  a  systematic  investigation  of  for  the 
model  used.  For  parallel  dislocations  we  have  inves¬ 
tigated  the  effect  on  p1[=  (fip^^/N^]  of  other 
mechanisms  -  the  manner  in  which  p-^  changes  with  the 
nature  of  these  other  mechanisms  and  with  the  relative 
value  of  their  resistivity  p^-  and  the  dislocation 
resistivity  (6p^)-^  -  and  thus  estimate  the  deviations 
from  Matthiessen ’ s  rule.  As  regards  random  dislocations, 
their  resistivity  p°  in  the  first  method  can  be  deter¬ 
mined  from  p^  itself  and  would,  of  course,  depend  on  the 
type  of  the  external  mechanism  present  and  on  p^/fip^. 

We  have  also  calculated  p°  by  the  second  method  in  order 
to  compare  the  results  of  the  two  methods'."  The  metal 
chosen  for  this  study  was  A1 .  This  choice  was  made,  on 

the  one  hand,  because  of  the  reliability  with  which 

+  3 

scattering  from  an  A1  ion  in  the  metal  can  be  deter¬ 
mined  from  the  effective  potential  theory  and,  on  the 
other,  because  of  the  extensive  experimental  study  which 
has  been  made  on  the  dislocations  in  this  metal.  Just 
as  done  by  Harrison  (1958),  the  scatterers  constituting 
the  dislocations  were  taken  to  be  vacancies.  The  other 
mechanisms  in  the  presence  of  which  the  restivitity  due 
to  dislocations  was  investigated,  were  taken  to  be 
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impurities.  For  simplicity  and  because  of  the  lack  of 
information  on  the  structure  factor  in  solids,  thermal 
scattering  was  not  included.  The  impurities  we  tried 
were, Mg,  Zn  and  Cd  whose  form  factors  in  A1  (host)  were 
calculated  in  Chapter  3* 

f 

For  carrying  out  the  numerical  calculations  for 
parallel  dislocations,  it  is  convenient  to  work  with 

»  »t 

integrals  1^  and  1^  defined  by 

2  2 
I ^  =  /  x  [w(x)J2  P£(l  -  |— )  dx 


I 


it 

£ 


f  x  [Aw ( x )]  2  P  (1  -  )  dx  , 

0  x  ^ 


(7.21) 


T  If 

which  occur  in  c^  and  c^  respectively  rather  than  working 
directly  with  the  coefficients.  Here  w  is  the  form 
factor  for  A1  and  Aw,  the  difference  in  the  form  factors 

f  " 

of  A1  and  the  impurity  ion  in  A1 .  From  1^  and  1^, 
can  be  determined  from  (7*9)  by  noting  with  the  help  of 


(7.10)  and  (6.44)  that 


D  ( c 


£ 


= 


a  D  ( 


£ 


a 


cl 

—  ) 

a 


(7.22) 


(  a  independent  of  l) ,  so  that  if  a  and  6  are  taken  to 

|  t  It  It 

be  equal  to  c^/I^  and  c^/I^  respectively  and  y  equals 
3/a,  we  get 


D(<q,  Op  =  a  D(ip  Y  Ip  . 


(7.23) 
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Note  that  y  determines  the  relative  value  of  the  resis¬ 


tivities  due  to  dislocations  and  impurities: 

~yl"/In  .  Thus  if  we  determine  D(I^  ,y  I0)  for  various 


£ 


values  of  y,  (7.23)  combined  with  (7*9)  gives  p^  for 
various  values  of  p^/6pD.  Writing  p^  explicitly,  we 
have 

.2 


pi  = 


_  3tt  -ti  /dE^  2  ^o  ^  f  -p f  T"N 
dk  F  2  D  (-t-JL*  Y  ZJ^ 


o  2 1 

2e  k. 


a 


3tt2  (m* ) 2 

o  2^3 
2e  n 


U‘ 


o 


kFa2 


D(l’,  Y  i") 


i 


(7.24) 


|  If 

D(I^,  y  I  )  for  different  values  of  y  was  calcu¬ 
lated  on  the  IBM  360/67.  The  number  of  integrals  1^  we 

,  tt 

used  was  0  and  that  of  1^  was  10.  The  results' for 

tt  .  * 

=  y  I  /I  for  the  three  impurities 


various  values  of  R  e 
are  shown  in  Table  IV.  The  expression  (7.13)  which,  we 


recall,  corresponds  to  the  limit' R  ->  <»,  gives  D=0.0552  . 

It  is  clear  from  Table  IV  that  D  increases 
progressively  with  R  so  that  p^  increases  as  more  and 
more  impurities  are  added.  The  two  limiting  values  of 
p^  are  0.72  x  lO-^  and  1.0  x  10-1^  ft  cm^ .  The  increase 
in  p^  with  R  has  a  simple  explanation  on  the  basis  of  the 
variational  principle.  The  effect  of  adding  the  impurities 


is  that  the  correct  variational  function  cj)£  becomes  a 

ir 
k 


mixture  of  (J)£1S  and  cf)imp  ,  where  the  superscripts  are 
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Table  IV.  Values  of  D(I^,  y  I^)  for  various  values 
of  R  for  dislocations  In  A1  with  Mg,  Zn 
and  Cd  as  Impurities.  Here  R=  I^/I T o~p^/6pD 
where  p-j-  is  resistivity  due  to 
impurities  and  6p^,that  due  to  dislocations. 


R  x  (l . 64  )  1 

D 

.2 

in  ryd 

Mg 

Zn 

Cd 

0.0 

0.0401 

0.0401 

0.0401 

1.0 

0.0468 

0.0462 

0.0472 

2.0 

0.0492 

0.0486 

0.0495 

3.0 

0.0504 

0.0499 

0.0507 

4.0 

0.0512 

0.0507 

0.0514 

6.0 

0.0522 

0.0518 

0.0523 

8.0 

0.0527 

0.0524 

0.0529 

10.0 

0.0531 

0.0528 

0.0532 

,  12.0 

0.0533 

0.0531 

0.0534 

16.0 

0.0537 

0.0535 

0.0538 

20.0 

0.0539 

0.0538 

0.0540 

self-explanatory.  This  gives  rise  to  deviations  from 

Matthiessen f s  rule,  which  are  positive  in  sign  (see 

§4.1),  and  hence  the  increase  in  p^.  (The  upper  limit 

of  p^,  of  course,  corresponds  to  the  use  of  the  trial 
imp 

function  .  ) 

As  regards  the  resistivity  p°  of  random  dis- 

2 

locations,  its  value  in  the  first  method  is  just  -j  p-^ 


. 
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so  that  In  the  two  limits  R  -*  0  and  R  ->  °°,  p°  is 

0.48  x  10~19  and  0.66  x  lO-"*"9  57  cm^ .  These  may  be 

compared  (see  Table  III)  with  the  most  recent  experi- 

o 

mental  value  1.8  57  cm  given  by  Rider  and  Foxon  (1966). 
Considering  the  crudeness  of  the  model  of  the  dislo¬ 
cation  used  in  the  theory,  the  agreement  with  experiment 
is  satisfactory. 

o 

We  have  also  obtained  by  the  second  method 
and  it  was  found  to  be  0.66  x  10  ^57  cm  .  This  happens 

to  be  the  same  as  p°  of  the  first  method  for  R+°°  .  This, 
of  course,  must  be  accidental  because  as  is  clear  from 
(7.14),  [with  D  given  by  (7«13)J  and  (7.17),  the  two 
expressions  for  p^  are  quite  different-even  for  R  ■+  00 . 

§8.4  Summary  and  Conclusions 

Using  a  simple  model  of  the  dislocation,  where  it 
is  treated  to  be  a  row  of  weak  scatterers  with  a  spheri¬ 
cally  symmetric  potential,  we  have  made  a  systematic 
investigation  of  the  resistivity  due  to  parallel  and 
randomly  oriented  dislocations  in  the  presence  of  other 
mechanisms  of  scattering.  By  obtaining  an  exact  solution 
of  the  Boltzmann  equation  for  parallel  dislocations,  we 
find  that  on  account  of  the  deviations  from  Matthiessen ’ s 
rule,  the  apparent  resistivity  of  dislocations  increases 
as  the  scattering  due  to  the  other  mechanicms  becomes 
larger  and  larger  compared  with  that  due  to  dislocations. 
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Numerical  estimates  made  on  dislocations  in  Al,  taking 
the  scatterers  to  be  vacancies,  show  that  the  two 
limiting  values  (R  -*  0  and  R-*-°°  ,  where  R=pI/6pD)  of 
the  component  of  the  resistivity  tensor  perpendicular  to 
the  dislocations ' differ  by  30%. 

As  regards  numerical  values,  the  calculated 
resistivity  p°  of  disordered  dislocations,  agrees  with 
experiment  sufficiently  well  to  encourage  further  in¬ 
vestigations  along  the  same  lines.  In  particular,  a 
repetition  of  the  calculations  presented' here  using 
better  estimates  of  the  scattering  due  to  the  core 
region  of  the  dislocation  would  be  desirable. 

On  the  experimental  side,  it  would  be  of  interest 
to  have  more  information  on  the  deviations  from  Matthies- 
sen's  rule  as  a  function  of  temperature  so  that  one  has 
a  better  idea  as  to  which  of  the  two  methods'  of  averaging 
(see  §7*2)  is  more  appropriate  in  a  given  temperature 
region.  This  will  enable  one  to  make  with  some  confi¬ 
dence  estimates  of  anisotropy  in  the  resistivity  of  a 
sample  where  dislocations  are  only  partially  aligned. 
There  is  some  experimental  evidence  for  anistropy  but  in 
the  present  state  of  the  theory  a  comparison  is  not  very 
meaningful . 


’ 
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Appendix  I 

Evaluation  of  the  Integral  I(q)  occurring  in 

Nonlocal  Screening 

From  (2.88),  I(q)  is  given  by 

I C q )  =  P  /  dhg+q|wR|g)  s  (i.i) 

k  - | k+q | 

where  P  denotes  the  principal  value  (henceforth,  we  shall 
leave  out  P,  but  it  will  be  understood  that  we  are  interes¬ 
ted  in' the  principal  value  of  the  integral). 

Two  methods  are  available  for  evaluating  I(q);  one 
is  due  to  Harrison  (1966)  and  the  other  due  to  Animalu 
(1965a).  Below  we  shall  discuss  both  of  them. 

Harrison’s  Method 

Let  us  choose  a  cylinderical  set  of  coordinates 
with  z-axis  along  q.  Let  the  coordinates  of  k  be 
(kp  p  ,  <f> ,  kpZ).  (p  and  z  have  been  defined  in  this 
manner  to  make  them  dimensionless.)  Now  the  various 
quantities  k,|k+q|  and  cos0  occurring  in  (k+q|w  |k) 
can  be  expressed  in  terms  of  q,  z  and  p  as  below: - 

1  1  /~27  2“ 

k  =  kp  /z  +  p 

| k+q |  =  k FJ  (z+q/kF)2  +  p2  (1.2) 

,  2,  ,->i2  2  k  +kFqZ 

k  + 1 k+q 1  -q  _  _ 

2k | k+q |  k | k+q 


COS0 


. 
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Then 


(k+q | I k) 


R,  n 

=  w  ( q  j  z  j  p ) 


(1.3) 


Substituting  (1.3)  and  the  first  two  of  (1.2)  in  (1.1) 


we  get 

kF  +1  dz  1_z2  R 
I(q)  "  f_±  q/k~p+~2"z~  ;Q  w  * 


(1.4) 


where 


y  =  P 


(1.5) 


The  integration  over  y  can  be  done  straightforwardly.  But 

in  the  z-integrat ion  we  come  across  a  singularity  at 

z=-q/2kp  if  q  <  2kp.  Harrison’s  method  consists  in  making 

a  polynomial  expansion  in  (z-q/2kp)  of  the  remaining 
1  — z^  R 

factor  ,  w  (q,z5y)  in  the  region  of  the  singularity 

/  ay 

0 

and  integrating  each  term  through  the  singularity  exactly. 


In  detail  the  various  steps  are  as  follows 

2 


1  —  z 1 


1. 


R 


For  a  given  q,  evaluate  Y  =  j  dy  w  (q,z,y)  for 

o 

various  values  of  z  in  the  interval  -1  to  +1. 


2. 


Split  the  integral  I(q)  into  three  parts 


J(q)  = 


-q/2kw-e 

/ 

-1 


dz  Y( z ) 
•q/2kp+z 


+ 


1 

/ 

•q/2kp+e 


dz  Y( z ) 
■  q/2kp+z 


+  / 


q/2kp+£  dz  y(z) 
-q/2kp-e  -q/2kF+Z 


(1.6) 
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3.  If  we  choose  e  to  be  sufficiently  large  that  1/e 

can  be  handled  conveniently  by  the  computer,  the  first 
two  terms  in  (1.6)  can  be  easily  evaluated.  To 
evaluate  the  third,  we  express  Y  as 

Y  =  a  +  bz?  +  cz'2  +  dz’^  +  .  ,  (1.7) 

where 

z’  =  z  -  q/2kp  .  (1.8) 

Then  the  last  term  in  (1.6)  becomes 

f  — -  (a  +  bz’  +  dz’  +  dz’  +  ....)  , 

z 

-e 

2 

which  is  equal  to  2e(b  +  de  /3  +  ....). 

Animalu’s  Method 

In  Harrison's  method  we  circumvented  the  difficulty 
of  integrating  through  the  singularity  by  expressing  the 
integrand  as  a  polynomial  around  the  singularity  and  then 
integrating  analytically.  In  Animalu’s  method  one 
evaluates  the  integral  I(q)  by  the  Gaussian  6-point 
formula  and  tries  to  pick  up  the  points  symmetrically 
about  the  singularity.  The  result  would  be  the  required 
principal  value. 

Let  (k,  a,  6)  be  the  spherical  polar  coordinates 

of  k  with  respect  to  the  system  of  coordinates  whose  z- 

|  -+■  -+■  | 

axis  is  along  q.  Now  since  | k+q |  and  cos0  occurring  in 


f§  ■ 


(k+q|w  |k)  can  be  expressed  in  terms  of  k,  q  and  m 
where  m  stands  for  cosa,  from  (1.1)  we  get 


I(q)  =  - 


(1.9) 


The  denominator  in  (1.9)  vanishes  on  the  hyperbola 


2km  +  q  =  0 


(I. 10) 


so  that  in  the  Gaussian  formula  the  points  cannot  be 
picked  up  symmetrically  about  the  singularity  with  k 
and  w  as  the  variables  of  integrations.  If  we  make  the 
substitution 

y  =  k 

(I. 11) 

x  =  2kco  +  q  , 

The  points  A,  B,  C,  D,  E,  F  of  Fig.  2(a)  go  over  to  A’, 

B',  C',  D’j  E’,  F’  of  Fig.  2(b).  Thus  the  new  line  of 
singularity  is  a  straight  line  and  this  makes  it  possible 
to  pick  up  the  points  symmetrically  about  the  singularity. 


The  Jacobian  of  the  transformation  (I. 11)  is 


T  =  |3  (a>,k)  |  =  1_ 
( x  j  y  )  2y 


(1.12) 


so  that  (1.9)  becomes 


I(q)  =  -  ^  ff  dx  dy 


x 


(1.13) 


X 
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k 


The  plot  of  (I. 10)  for  q  <  k^  . 


Fig.  2(b) 

Remapping  of  the  points  A}  B,  C,  D,  E,  F  of  Fig. 2(a) 
with  the  transformation  (I. 11).  The  line  E'Q  is 
drawn  such  that  A'F'  =  F'Q.  The  points  P  and  R  on 
E’Q  and  C’B’  respectively  are  such  that  E T C * =E ’ P  and 
PQ=RB  f . 
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To  pick  up  the  points  symmetrically  about  the  line 
singularity  E'F',  we  draw  a  line  E’Q.  If  in  the 
region  A’E’Q  we  perform  the  x-int egrat ion  first,  we 
shall  have  achieved  our  objective.  To  take  care  of 
the  point  singularity  E'  we  take  E’C’P  as  another  region 
of  integration  but  here  we  perform  the  y-integrat ion 
first.  The  remaining  domain  is  free  from  singularities 
but  for  convenience  of  choosing  the  limits  we  divide 
it  into  two  parts  by  a  line  PR.  In  C1  PR  and  RB’QP  the 
order  of  integration  is  immaterial.  Supposing  we  do 
the  x-integrat ion  first,  and  denote  the  regions  A'E'Q, 
E’C'P,  C 1  PR  and  RB’QP  by  1,2,3  and  4  respectively,  (1.13) 
may  be  written  as 


I(q)  =  1^  +  I9  +  1.3  +  1^, 


2  3 


where 


h  -  - 


h  =  - 


k 


TT 


P 


2y-q 


%  f  ydy  / 

q  q/2 


dx  R  /  s 

w  (q,x,y) 


x 


IT  S  dx  Q+X)/2  R 


±  s 


f 


q  0  X  ( q-x )/2 


-(2y-q) 

ywlv(q,x,y)  dy 


h  = 


TT 


q 


q+2y 


£  /  ydy  f 

q  0  q 


dx  R ,  \ 

go  ( q ,  x  ,  y  ) 


x 


h  “ 


ir  >  .  2y+q 

-  /  ydy  / 
q  q  2y-q 


dx  R ,  n 

—  m  ( q ,  x  ,  y )  . 


x 


. 


.  1 1  ±  :  3  3  sniflLi  3  no  J  1 6*1 3 9  J ■ n  1  - 1 
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Note  that  in  Pig.  2  we  have  taken  q  <  . 

The  discussion  of  the  case  k^  <  q  4  2kp  is  exactly 
similar.  In  fact  as  long  as  correct  order  of  inte¬ 
grations  is  left,  the  computer  deals  with  both  of  these 
cases  in  one  go.  The  case  of  q  greater  than  2kp  does 
not  present  any  trouble  because  then,  the  singularity 
lies  outside  the  region  of  integrations  and  the  whole 
domain  can  be  included  in  one  step.  However,  it  may 
be  pointed  out  that  for  q  >  2k^,  though  the  singularity 
lies  outside  the  domain  of  integrations,  it  affects  the 
accuracy  somewhat  because  the  points  in  the  Gaussian 
formula  are  not  symmetrical  about  the  singularity. 

The  effect  on  the  accuracy  is  greater  the  closer  is 
the  singularity  to  the  region  of  integrations.  Thus, 
for  larger  values  of  (q-2kp),  the  accuracy  is  better. 

Comparison  of  the  Two  Methods 

We  have  tried  both  the  methods  and  find  that 
Animalu’s  method  is  much  easier  to  use.  However,  it 
may  be  added  that  the  accuracy  of  the  two  methods  is 


of  the  same  order. 


. 


■ 


X 
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Appendix  II 


Proof  of  the  Identity  (6.33) 

(6.35)  is  readily  proved  by  multiplying  (6.32) 
1  1 

by  pL+2r(w)  and  by  3  and  uslnS 


+  1 

/ 

-1 


pi(U) 


1 

) 

L 


(u>) 


a(TT  6LLf 


3 


where,  as  in  the  text,  a ( L )  =  ( 2L+1 )/ [2L( L+l )J . 
gives 


nk^2r)(L)  =  a ( L+2r ) 


+  1 

/  U) 
-1 


2n 


p  i(“) 


p3"  ,  0  (oj) 

L+2r  ' 


dm 


This 


(ii. i) 


and 

nk^2r)(L)  =  a(L-2r)  /  io2n  p£((o)  P2_2r,(o)  do  . 


Changing  L  to  (L-2r)  in  (II. 1)  gives 

kn(2r)(L-2r)  =  a ( L )  /  w2n  ?}  0  (a))Pd(oo)  dm  . 
n  1  ,  L-2r  L 


(II. 2) 


(II. 3) 


Substituting  for  the  integral  in  (II. 3)  from  (II.2), 
we  get 

a(L-2r)nk[2r)(L-2r)  =  a(L)nk^2r } (L) ,  (II. 4) 


which  is  the  identity  (6.35),  we  wanted  to  prove. 


X 
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To  Verify  that  (6.39)  is  a  Solution  of  (6.38) 


We  recall  that  (6.38)  and  (6.39)  give 


(co+ao)fL+  al  ^lk  (L)fT,+iki  2  )fT,_p  +  i  ^p^L+2  )fT.+p] 


L  1  1 


L-2  12 


L+2- 


+ - +a  [  k°(L)f  +  k(2)(L-2)f  +  k^2)(L+2)f 


A  Hi.  L  i  1 

00  o 


L-2  £  “2 

o 


L+2 


+ - +  „  kn(2V(L-2«,JfT  +  „  k^2ilo)(L+2<>„)fTJ.0„  ] 


l  1 
o 


o'  L-2£  £  2 

o  o 


o '  L+2  £ 


o 


=  A  6 


LI 


and 


(II. 5) 


+1 


fL= 


2L+1  A  f 

2L(L+1)  (c£+a„)  ' 


P2(co)P2(oo)  du) 


2  4  2f 

o  -l  1+y^o)  +y2^  +...+Y^  w  o 


(II. 6) 


o 


We  want  to  show  that  (II. 6)  statisfies  (II. 5).  For  this 
purpose  we  note  that  from  (II. 4),  (II. 6)  and  (6.32),  it 
follows  that 


1k°(L)fL  +  1k[2)(L-2)fL_2  +  1k^2)(L+2)fL+2 


2L+1 


A  go2P^  ( 03 )  ( 03 )  doo 

A  f  -  L  1 


2L ( L+l )  c"+a  '  ,  2,  ,  22£  5 

v  00  l+Y-,w  +....+y«  w  00 

^o 


2k°(L)fL  +  2k^2)(L-2)fL_2  +  2k^2)(L+2)fL+2 


(II. 7) 


+  i|k|4)(L-4)fL_4  +  2kC4)(L+4)fL+1| 


2  L+l  A 
2L(L+l )  cM  +a 


I 


034P^  ( 03 )  P^  ( 03 )  do) 


2  £ 


o  o  l+Y-^w  + .  .  .  .  +y ^  03  o 


o 


(II. 8) 


u.b  (U>)x'  (*>)j  *M  __  A  if.,  : 
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and  so  on.  On  substituting  (II. 7),  (II. 8 
(II. 5)  it  is  immediately  seen  that  (II. 5) 
fied  since 

2L(L+1)  ^  PL ^ w ^  =  6L1 

Thus,  with  the  help  of  (II. 4)  and  (6.32), 
shown  that  (II. 6)  satisfies  (II. 5). 


)  etc.  in 
is  satis- 


we  have 
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Expl  anation  of  the  Tables 


The  following  tables  give  the  form  factors  for 
eleven  alloys  calculated  from  the  Harrison  pseudopoten- 
tial^.  For  the  impurity  ion  a  proper  account  has  been 
taken  of  the  change  in  the  surroundings^  of  the  ion 
compared  to  the  case  when  it  was  in  the  pure  metal. 

Since  most  applications  depend  on  the  relative  values 
of  the  host  and  the  impurity  form  factors  and  since  var¬ 
ious  pseudopotential  calculations  differ  considerably 
in  detail  (e.g.  in  the  calculation  of  the  screening, 
inclusion  of  the  effects  of  exchange  and  correlation, 
inclusion  of  the  influence  of  the  effective  mass,  and 
so  on)  we  also  give  the  values  of  the  form  factors  of 
the  host  obtained  in  a  manner  similar  to  the  results 
for  the  impurity. 

The  eleven  alloys  for  which  the  calculations  were 

made  are  (Mg,  Zn,  Ca)  in  Al;  (Al,  Li)  in  Mg;  (Na,  K)  in 

K;  (Li,  K)  in  Na.  For  each  host  we  first  give  the  form 

factors  of  the  impurities  and  then  for  the  host.  For 

each  ion  the  first  column  is  the  reduced  unit  q/k^  where 

•  * 

k.-,  is  for  the  host.  The  next  two  columns  give  the 

r 

screened  form  factors  in  rydbergs.  The  first  set  of  form 
factors  was  obtained  by  setting  m*/m  =  lj  in  the  other 
the  correction  for  the  influence  of  effective  mass  has 
been,  applied . 
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